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Abstract 


Techniques are described for separating gas from glacier ice and analyzing for oxygen, 
nitrogen, and carbon dioxide. A systematic survey of a temperate glacier for enclosed gas 
showed 1) no general loss of gas from the glacier, but a specific loss of oxygen from the gas 
from the upper to the lower part of the glacier, and 2) the gas pressure and ice density 
increased from the upper to the lower part of the glacier. Present evidence indicates that a 
chemical reduction is responsible for the oxygen loss. The various data are discussed in relation 

to glacier mechanics. 


Introduction 


A prominent characteristic of glacier ice is 
the inclusion in the ice of gas bubbles or 
cavities. These gas enclosures may originate in 
two different ways (cp. SCHOLANDER, et al., 
1956). In high-arctic glaciers where melting 
does not occur in the ice formation area, 
glacier ice is formed by snow compaction. The 
gas enclosed in such ice is simply atmospheric 
air trapped between the fused snow crystals, 
and hence of complete uniformity in composi- 
tion. The atmosphere is 20.94 % O, and 0.03 % 
CO, and the rest nitrogen including noble 
gases (cp. Hock, et al., 1952). 

However, if the freezing of water plays a 
role in the ice formation, one would expect 
the gas to have a composition other than 
atmospheric air. Water fully saturated with 
air at 0° C and at 760 mm. Hg contains 2.9 % 


1 These studies were aided by a contract between 
ONR, Department of the Navy, and the Arctic Institute 
‚of North America. 

1 Contribution No. 873 from the Woods Hole 
Oceanographic Institution. 
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by volume of dissolved gases of which 34.9 % 
is O, and 1.75 % is CO, (Handbook of Chem- 
istry and Physics, 1952 ed.). As the solubility 
of gases in ice is at least 1000 times less than in 
water, probably zero, these gases will separate 
from the ice when the water freezes (SCHOLAN- 
DER, et. al., 1956). Hence when air-equilibrated 
meltwater percolates down through snow on 
a glacier and freezes to ice it will enrich the 
enclosed gas with oxygen and carbon dioxide, 
and one would expect considerable local as 
well as seasonal variations. 

Laboratory experiments have indicated that 
CO, diffusion through pure ice is some 40,000 
to 70,000 times slower than through water and 
O, is even slower (SCHOLANDER, et al., 1953). 
Considering the enormous distances in a 
glacier and the relatively large amounts of 
trapped gas it seems likely that diffusion losses 
would be extremely small even over long 
periods of time. If there are no gas changes 
through oxidation, it is possible that very 
ancient atmosphere will be preserved in the 


depths of the high-arctic inland ice. 
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Several workers have determined the com- 
position of this gas in isolated pieces of glacier 
ice (STEENSTRUP 1883, HAMBURG 1895, BERNARD 
1922, BARNES 1927). Though in general the 
techniques employed for separating the gas 
from the ice do not yield accurate data, these 
indicated that the composition of the gas was 
mote or less close to air. 

Mire recently BADER (1950) investigated the 
pressure and amount of gas in the bubbles of 
the terminus of the Malaspina Glacier, Alaska. 
The pressure varied between 11/,—21/, atmos- 
pheres and the gas quantity between 3.8 - 
4.8 % of bulk volume. 

SCHOLANDER, KANWISHER, and NUTT (1956) 
measured the gas composition (O, and N,), 
quantity, and pressure in pieces of glacier ice 
taken from icebergs along the Labrador coast. 
They concluded that the gas in four out of six 
icebergs examined was indistinguishable from 
present day atmosphere, but that the other two 
had a composition significantly lower in oxy- 
gen (20.0%). They suggested that this could 
possibly signify a low atmospheric O, at the 
time the snow compacted to form the ice, and 
he sp eculation was made that this might be 
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due to a lowered photosynthetic activity during 
a period of climatic decline and advanced 
glaciation. In these determinations the ice was 
melted under mercury, the gas phase and the 
dissolved gases were analyzed separately and 
the results totalled. This method had an accuracy 
of +0.2 %. The pressure was found to vary 
between 1.0 —6.2 atmospheres. 

Many questions arise with regard to the gas 
enclosures, such as: to what extent is the atmos- 
phere trapped in the ice preserved since the 
formation of the ice? to what extent do the 
bubbles reflect the pressures in the ice? to what 
extent is the ancient precipitation forming the 
ice preserved in its original composition? Only 
by systematic studies can one hope to answer 
these questions. 

In the present study we have undertaken such 
a systematic investigation on a temperate gla- 
cier, Storbreen, in the Jotunheimen district of 
Norway. This was selected for study because 
it is readily accessible from a log cabin well 
enough constructed to serve as a field labora- 
tory, and because it is being examined in other 
respects by the Norsk Polarinstitutt. The field 
investigation took place in the month of March. 


Fig. 1. Map of Storbreen, from 

surveys by the Norsk Polarinsti- 

tutt. The dark circles are sample 

locations, and the elevations are 
in meters, 
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Sampling of Ice 


Samples were taken one-half meter below 
the surface of Storbreen ice in a line from the 
terminus to well above the firn line (Fig. 1). 
The glacier topography determined the posi- 
tion of the line, but except possibly for the 
upper stations (8 — 9) the samples are from more 
or less the same line of flow. Samples TA and 
8A were taken 30—4o meters laterally distant 
from 1 and 8 respectively for duplication of 
the upper and lower ice. 

The samples were brought to the field 
laboratory for analysis. Duplicate samples were 
brought to the cold room at the Institute of 
Zoophysiology, Oslo, where studies were made 
on single bubbles and various chemical tests 
conducted. 


Methods 
1. Bulk determination of gas bubbles in ice 


In this method the gas is separated from the 
ice by shaving the ice in high vacuum and in 
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cold. The gas thus released by the opening of 
the bubbles is collected over mercury and 
analyzed. Data obtained give gas quantity, 
pressure, and composition, and ice density. 

Apparatus. A piece of ice is rough-hewn, 
placed in a small specially constructed lathe, 
and turned to a cylinder. The ice is transferred 
to a jig and the ends are carefully sheared flat 
and perpendicular with a wide-bladed knife. 
Dimensions and weight of the piece are accu- 
rately taken. The ice cylinder is inserted in the 
barrel of an ice shaving box (Fig. 2) where it 
occupies 75 % of the volume to allow room 
for the shavings. The top and bottom of the 
shaving box are vacuum sealed with standard 
O-rings. 

Inserted in the box is a stainless steel rotary 
knife designed like a wood drill of the German 
“forstner” type which cuts a hole with a flat 
bottom. The drill shaft extends through a 
bearing in the top of the box fitted with four 
O-ring vacuum seals, and is turned by an 
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Fig. 2. Sketches of the apparatus. A is a section through the ice shaving box, B a sketch of the ice lathe, 
and C a schematic diagram of the vacuum system. 
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ordinary hand brace. The rotary shaving can 
proceed until the shavings fill the space above 
the knife, leaving an intact piece of ice at the 
bottom. 

The shaving box is connected to a high 
vacuum pump through a glass manifold. Also 
connected on the manifold are 1) a U-tube 
mercury manometer, 2) a measuring cylinder 
and leveling bulb filled with ethylene glycol, 
and 3) a 500 cc. gas extraction flask and leveling 
bulb filled with mercury (Fig. 2). 

Procedure. The shaving box containing the 
prepared piece of ice is connected to the mani- 
fold and the system evacuated until the vapor 
pressure of ice is obtained (1.7 mm Hg at 
-10°C.). Air at ambient pressure is intro- 
duced from the measuring cylinder, giving a 
measure of the space evacuated before the ice 
is shaved!. After shaving the gas pressure in the 
system is rcad on the manometer. Part of the 
liberated gas is transferred to the gas extractor 
and from there to a § cc. gas pipette where it is 
stored under overpressure, and later analyzed. 
The last step is measuring by caliper the un- 
shaved piece of ice. All these proceedings are 
conducted at temperatures below freezing. 

Gas analysis. The gas was analyzed in a 1/, 
cc. volumetric analyzer for O,, CO,, and N,,” 
with an accuracy of + 0.015 % (SCHOLANDER, 
1947). Accurate analyses can only be made in a 
well thermo-regulated room and not below 
20° C., which was accomplished by means of 
accessory stoves. Checks on atmospheric air 
were run ahead of all ice gas determinations to 
insure that full accuracy was attained. 

Calculations. The volume of the cylindrical 
ice sample is calculated from the measure- 
ments, then divided by the weight to obtain 
the density (r). The amount of gas space (s) 
in the ice is 


(.917 -r) 
Ay 
+917 


s= 


where .917 is the density of pure ice (Handbook 
of Chemistry and Physics, 1952 ed.), and v the 
volume of ice shaved. 

The quantity of gas (Q) at STP is 


1 This can only be used in ice where the gas cavities 
do not interconnect (cp. check 3). 

2 N, as used herein contains the other non-absorable 
gases. 
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Pzyp 
-F 
( 760 ) 
where V, is the volume of the evacuated 
system, F the volume of the shaft introduced 
on shaving, P the pressure of the gas, vp the 
vapor pressure of ice and 760 atmospheric 


pressure in mm Hg. The pressure (P,) exerted, 
by the gas in the ice is then 


pa 


S 


Checks 1) Diffusion tightness of the system. 

Introduction of gas mixtures other than air 
into the system revealed a very slight diffusion 
leak of air through the vacuum tubing and 
joints. By reducing the amount of vacuum 
tubing to the least possible this error was 
reduced to an insignificant amount. 
2) Shaving ice in the presence of gas. 

It was thought that adsorption of gas on 
the enormous surfaces of the shavings might 
significantly affect the gas composition. To 
test this 5 cc. of air were introduced into the 
evacuated system containing nearly bubble-free 
lake ice, the ice was shaved, and the gas re- 
moved and analyzed. The O, remained un- 
changed (-+0.015 %) but there was a slight 
increase of CO, (to 0.09 %) which might 
have come from the few tiny bubbles visible 
in the ice, or perhaps was present dissolved in 
intercrystalline brine films (SCHOLANDER, et 
al., 1953). The same procedure was followed 
using a gas mixture of 15 % O, and 5 % CO», 
and there was a slight loss of CO, and a slight 
gain of O, corresponding to the above men- 
tioned diffusion leaks of air. It may be con- 
cluded that ice shavings do not measurably 
alter the gas composition. 

3) Gas tightness of the ice cylinder. 

Determination of the bulk gas content of 
the ice requires that the gas bubbles are dis- 
crete and that the enclosing ice structure is 
vacuum tight. To test this pieces of the samples 
were evacuated for five minutes in kerosene 
dyed with Sudan III. The volumes of the 
melted samples and the oil were then taken to 
determine the amount of oil that had entered 
the cylinder. In most cases the oil entered only 
those bubbles that were cut open at the sur- 
face, and in these cases the gas pressure estima- 
tion had a maximum error of —9 %. In samp- 
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les 1, 7, 8A, and 9 due to many channels inter- 
connecting the bubbles, the oil penetrated into 
the pieces, and the pressure measurements on 
these samples are considered invalid. 


2. Micro determination of composition of single 


gas bubbles 


In this method the gas bubbles are opened 
anaerobically, the gas removed in gas pipettes 
and transferred directly to composition analy- 
zers. 

A well is carefully shaved over a gas bubble 
with a spoon-shaped chisel. Concentrated acid 
citrate and then mercury are placed in the well. 
A 0.8 mm drill, cut off and resharpened so 
that only 2 mm of spiral grove remain, is 
inserted through the mercury and citrate and 
a hole drilled into the bubble. The gas is 
removed in a mercury micropipette inserted 
through the fluids, and either 1) transferred 
directly to a syringe analyzer, and analyzed 
for O,, CO,, and N, with +0.2 % accuracy 
in 40 mm? samples (SCHOLANDER, et al., 1955), 
or 2) stored behind mercury in a cup for ana- 
lysis by a micrometer micro-analyzer with 

+0.2 % accuracy in I mm? samples (ScHo- 
LANDER and EVANS, 1947). 


Results 


The results of the various measurements on 
Storbreen ice are presented in Figures 3 and 4, 
all plotted against the sample locations. 

Pressure (Fig. 3). The gas pressure in the ice 
may be seen to increase more or less regularly 
from atmospheric pressure in the upper part 
of the glacier to over 2 atmospheres near the 
terminus. 

Density (Fig. 3). The ice density fluctuates 
widely throughout the glacier (.854—.904 
‘gm/cc) and in samples from the same location. 
However there appears to be a trend of in- 
crease of ice density from the upper part to 
the lower part of the glacier. 

Gas Quantity (Fig. 3). The amount of gas 
varies widely throughout the glacier (1.8—7.8 
volume-percent) with no relation to location, 
i.e. there is no general loss of gas going down 
the glacier. 

Oxygen (Fig. 4). The O, composition of the 
gas in the upper part of the glacier (8, 8A, 9) 
is approximately equivalent to air, and the 
values vary from bubble to bubble about 2 % 
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Fig. 3. Plots of the gas pressure, ice density, and gas 
quantity against sample location. The open circles are 
samples 1 À and 8 A. 
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(19.5—21.9 %). Going toward the terminus 
there is a general decrease of O, to values as 
low as 4.8 %, and the bubble to bubble 
variations appear to increase. The terminal 
values are in two groups, the average of the 
lower group (5—12 %) continuing the curve 
of O, decrease down the glacier. The higher 
group of values (16—21 %) suggests that per- 
haps recently the thin terminal ice may have 
cracked and allowed air or meltwater to enter 
and contaminate the original ice. 

Carbon Dioxide (Fig. 4). The CO, in the 
samples near the firn line and at the terminus 
varies from air to values ten times greater 
than air (0.05—0.29 %). All the remaining 
samples have atmospheric CO, in the bubbles. 
There is no inverse correlation with the O, 
values. 


Chemical Tests of Ice Purity 


The amount of dry matter, including 
particulate matter, in the ice was determined 
by evaporating samples 1 A, 4, and 8 A to 
dryness. The residue of these samples weighed 
0.6, 0.7, and 0.5 mg. per 100 cc. ice respec- 
tively. 

The glacier ice is essentially frozen precipi- 
tation, and it is of interest to compare it 
with present day precipitation. The salt con- 
tent would also give an idea of the amount 
of brine that may be present as intercrystalline 
films (cp. Dorsey, 1940, SCHOLANDER, et. al., 
1953). The specific conductivity and pH of 
the melted samples were measured by standard 
methods. The conductivity was found to be 
very uniform between 7.0—11.5 ohms-lem-1 x 
x 1076, and the pH ranged 6.08—6.44. 


Discussion 


A quantitative study of the gas enclosures in 
Storbreen has given the following main results: 
that there is a loss of oxygen from near atmos- 
pheric in the upper part of the glacier to well 
below atmospheric in the lower; that the O, 
varies from bubble to bubble, and this varia- 
tion increases from the upper to the lower part 
of the glacier; that the O, decreases without 
an increase in COg,; and that the gas pressure 
increases from the upper part of the glacier 
towards the terminus. We shall below discuss 
these findings in relation to each other and to 
known data on glacier mechanics. 
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Oxygen loss. There are two processes that 
present themselves as possible causes for the 
O, loss: 1) diffusion loss through the ice to 
the atmosphere, and 2) oxygen loss by chemical 
action. 

If there is a diffusion loss it would most likely 
take place through intercrystalline brine films, 
and would then be a function of salt content! 
and temperature. If there is a positive total 
pressure gradient from the bubble to the air, 
oxygen would be preferentially lost because of 
its higher diffusion rate relative to nitrogen. 
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The O, tension in the bubbles, calculated from 
content and pressure, are given in Figure 5. 
It may be seen that in most cases there exists 
a positive pressure gradient which could be 
responsible for a preferential diffusion loss of 
oxygen. However in sample 1A the O, loss 
has proceeded to a tension below that of the 
ambient atmosphere, and this clearly is not 
due to diffusion loss to the atmosphere. We 
must therefore look for other causes for the 
oxygen loss. 

The increasing spread of individual bubble 
O, values going down the glacier suggests that 
the reduction is not uniform from bubble to 
bubble and is another indication ofthe consider- 
able diffusion tightness of ice. This agrees 
with the finding of gas bubbles containing 
40% CO,and 60% N, 3 cm below the surface 
of 6 months old pond ice (SCHOLANDER, et al., 
1953). 

If the O, loss were caused by bacteriul 
decomposition of organic matter one woald 
expect a concurrent increase in CO,. This 


! 100 cc. meltwater with a conductivity of 10 ohms-1 
cm! x 10-8 would leave about 10 mm? unfrozen brine 
at —10° C and so mm at —2° C. 
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intriguing possibility is unlikely, as most of the 
samples with reduced O, had CO, values 
equivalent to air, in fact the highest CO, 
encountered in any sample amounted to only 
a few tenths of one percent. 

In dissolved salts and acidity Storbreen ice 
compares well with present day precipitation. 
The range of conductivity and pH of precipi- 
tation taken from measurements published for 
four inland Swedish stations (EGNER and ErIKS- 
SON, 1955) are presented in Table 1. 


Table 1. Range of Conductivity and pH in Preci- 
pitation. Four Inland Swedish Stations. 


Conductivity 


Station pH 
ohms "cm 108 
Riksgränsen .. 6—54 5-6—6.3 
SET pee 9—19 5.2—6.2 
Mora Set. 19—90 6.8—7.3 
Ter eee Soe ré 6703 
Storbreen Ice 7—12 6.I—6.4 


It may be concluded that the evidence so far 
points to a chemical oxidation of dust taking 
place within the individual bubbles, and that 
the precipitation preserved in Storbreen is es- 
sentially the same as today’s. 

Carbon Dioxide. Slightly elevated CO, values 
occurred in the samples taken near the firn line 
(8, 8A, 9) and near the terminus (1, 1A). All 
samples in between had a CO, content equi- 
valent to air. Although the analytical accuracy 
is poor for such low values it seems clear that 
the CO, tension in most of these intermediate 
samples is higher than in air, and that diffusion 
equilibrium has not been attained. 

Firnification Process. A fresh fallen snow mass 
has low density and is permeable to air and 
water. Through various processes the snow 
crystals become polygonal granules and the 
snow density may increase to .4—.5 in one 
year, and in succeeding years these granules 
compact to ice by the addition of load, the 
spaces between the granules becoming sealed 
from each other and from the atmosphere at 
about density .8 (cp. BADER 1939, 1950, 
MATTHES 1942). 

Examination of the firn at three meters’ 
depth at location 9 showed alternating layers 
Tellus VIII (1956), 4 
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of coarse-grained snow and bubbly ice. Sample 
9 was of this ice and although the oxygen was 
not elevated there was a slight elevation of 
carbon dioxide indicating gas separated from 
water on freezing. However most other samples 
had CO, values corresponding to air, and it 
appears therefore that simple trapping of air 
by compaction of ice granules outweighs by 
far gas formed by freezing of meltwater. 
Relation to Glacier Mechanics. In Figure 6 is 
given a schematic sketch of the lines of flow 
of a temperate glacier (KOECHLIN, 1944). In 
this concept the ice moves through the glacier 
along specific levels and is thereby subjected 
to different loads, the ice nearest the terminus 
having been subjected to the greatest load. It 
would appear that the gradient found in gas 
pressure agrees with this concept, and the 
surface samples reflect the compaction in the 
glacier interior. However the pressures are not 
equal to the theoretical ice load, and as Bader 
concluded, some of the pressure appears to be 
lost as the ice surfaces from inside the glacier. 
It is also evident from this concept that sur- 
face sampling down the glacier would sample 
ice of increasing age, and the gradient in 


PRN 
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Fig. 6. Schematic diagram of a longitudinal section 
through a temperate glacier, after KOECHLIN, 1944. 
The arrows indicate the lines of flow. 


oxygen content found in Storbreen also agrees 
with this. If the oxygen reduction proceeds at 
a uniform rate and the rate were known, the 
ages of the samples could be estimated. For 
example, the O, in sample 9, estimated by 
other criteria as some 15 years old (Liestol, 
personal communication), was about 0.3 % 
below air, corresponding to a reduction of 
0.02 % per year. In the terminal samples the 
oxygen averages about 12 % below air and 
hence the age of this ice might possibly be some 
600 years. 


Summary 

A systematic investigation of the gas enclo- 
sures in a temperate glacier has been under- 
taken. The apparatus and method are described 
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for separating the gas from the ice in bulk and 
analyzing it for oxygen, nitrogen, and carbon 
dioxide with +0.015 % accuracy. This meth- 
od gives data on gas composition, quantity, 
and pressure, and ice density. A technique is 
also described for extracting and analyzing 
the gas in individual bubbles with +0.2 % 
accuracy. 

Samples were taken from the surface of the 
glacier in a longitudinal line from above the 
firn line to the terminus. It was found that: 
1) the gas pressure increased from atmospheric 
near the firn line to three times atmospheric 
near the terminus, 2) there was a trend of 
increasing ice density from the upper to the 
lower part of the glacier, 3) there was no 
general loss of gas from the glacier, 4) the 
oxygen content was equivalent to air in 
the upper part of the glacier, but there was a 
loss of oxygen going toward the terminus, 
5) the oxygen content varied from bubble to 
bubble, and this variation increased going down 
the glacier, and 6) the oxygen was lost without 
an increase in carbon dioxide. 

The possible causes for the oxygen loss are 
discussed. Loss by diffusion through the ice 
appears unlikely, and the present evidence 
points to dust oxidation within the individual 
bubbles. The glacier ice was tested chemically, 
and compares well with present day precipita- 
tion. 
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The various gas ‘parameters are discussed in 
relation to glacier mechanics. It appears that 
simple trapping of atmospheric air is responsible 
for most of the gas in the glacier, and that gas 
frozen out of meltwater is of minor occurrence. 
The various data support the current concept 
of glacier movement. 
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Shower Formation in Large Cumulus 


By F. H. LUDLAM and P. M. SAUNDERS. Imperial College, London! 


(Manuscript received July 14, 1956) 


Abstract 


A distinction is drawn between the appearance of evanescent fibrous streaks in decaying 
cumulus towers (fibrillation), and the glaciation of cumulus. A technique of measuring cumulus 
is described, and used to show that fibrillation, which marks the onset of shower formation, 
occurs when cumulus towers reach a fibrillation zone or a fibrillation level; towers which rise 
appreciably above this zone glaciate if supercooled. The fibrillation level varies from occasion 
to occasion, apparently according to the strength of the cloud updraughts, and in a manner 
which indicates the coalescence process to be the general shower-forming agency in clouds 
observed in Sweden, although the temperatures at the summit-level of incipient showers 


may be as low as — 20° C. 


I. Introduction 


It is well known that the size which cumulus 
clouds attain before producing showers varies 
from place to place and occasion to occasion. 
Sometimes, especially over or near the ocean, 
showers fall from cumulus only about two 
kilometres deep (see, for example, ALPERT, 
1955), while at other times, especially in the 
interior of continents, no showers form in 
clouds having a depth of several kilometres. 
Attempts have been made to relate shower 
formation to the temperature at the level of 
the cloud summits, but this, just as the critical 
cloud depth, seems to vary between wide 
limits. For example, in New Mexico Bra- 
HAM, REYNOLDS and HARRELL (1951) detected 
radar echoes from only one fifth of summer 
cumulus whose summit-level temperatures 
lay between —12 and — 24° C; PEPPLER (1940) 
concluded from observations made on routine 
aircraft soundings in Germany that showers 
usually occur when the summit-level tem- 
peratures fall below —12° C; ScHWERDTEE- 


! This work was done while the writers were on 
attachment to the International Institute of Meteorology, 
Stockholm. 


GER (1948) found —3°C to be the most fre- 
quent value of this temperature in over 250 
aircraft observations of shower clouds near 
60°N 13°W, while it is now generally recog- 
nised that in the tropics showers commonly 
form in clouds whose summit-level tempera- 
tures are several degrees above the melting- 
point. Such observations cannot be recon- 
ciled with the idea that showers form only 
when the cloud tops become cold enough 
for ice nuclei to promote the appearance of 
significant concentrations of ice crystals 
amongst supercooled water droplets, as en- 
visaged in the Bergeron-Findeisen theory of 
shower formation, and it is now accepted that 
processes of droplet coalescence are important, 
not only in the small shower clouds of trop- 
ical regions, but probably also in those 
clouds large enough to have strongly super- 
cooled tops which become “glaciated” in 
the course of shower production. LUDLAM 
(1951, 1952) emphasised that processes of 
cloud-particle growth require a certain period 
to accomplish the production of particles of 
precipitation-element size, and regarded the 
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period available in a particular cumulus as 
determined by the time taken for the cloud 
updraughts to lift the small cloud particles 
from near the base, where they are formed 
(or some higher level in the case of ice crys- 
tals), to the level of the summits, where they 
are rapidly evaporated. Thus the vertical ex- 
tent, and hence also the summit-level tem- 
perature, of those cumulus which become 
shower clouds should depend very much on 
the strength of the updraughts in their interiors, 
a parameter not previously considered. The 
observations discussed in the following para- 
graphs were made to examine more closely 
the size of cumulus in which showers form, 
and the relation to the speed of the cloud 
updraughts. 


2. The significance of the fibrillation and 
glaciation of cumulus tops 


In an ordinary population of cumulus the 
clouds reach a variety of heights before sub- 
siding and shrinking or dissolving. Usually 
the number of summits attaining a particular 
level falls rapidly as this level is raised, and this 
is true also of the individual towers, or bubbles 
(LupLAM and Scorer, 1953), of which the 
larger clouds are composed. 

During the diurnal development of cumulus 
overland, showers do not form unless the 
largest clouds reach a certain size, or summit- 
level. At first the showers are slight and widely 
scattered; it is uncommon for the develop- 
ment to cease just at this stage: usually they 
very soon become more intense and numerous 
as the clouds become still larger and greater 
numbers exceed the critical size. In temper- 
ate latitudes the shower clouds are nearly 
always more than 2 km deep, and so the tem- 
perature in their tops is usually at least a few 
degrees below the melting-point. The forma- 
tion of a shower in one of the larger clouds is 
characteristically accompanied or preceded 
(often several minutes before precipitation 1s 
seen below the cloud base) by a transformation 
in the appearance of its upper parts. At first 
the cumulus tops are well-defined, with sharp 
outlines and a wealth of details provided by 
many minor domes or nodules and their 
shadowed interstices; if the cloud tower 
continues to rise its upper cap retains a sharp 
and bulging outline, but soon on its sides 
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and lower parts the details dissolve into a silky 
sheen, and the edges become more diffuse 
(as in Figs. s—7). Eventually the whole of 
the tower acquires the lustre and the fibrous 
texture which is typical of ice clouds, and 
the largest clouds often develop persistent 
fibrous anvils which from optical phenomena 
and flight observations are known to be 
composed largely, if not wholly, of ice par- 
ticles. Consequently this striking transforma- 
tion is called glaciation. 

Sometimes, in a cumulus tower which has 
reached the peak of its ascent, traces only of 
this transformation can be seen (as in Figs. 
1—3). It may be difficult or even impossible to 
detect until the evaporation of the tower 
begins; usually this is accompanied and 
hastened by a subsidence, and within a minute 
or two the tower becomes emaciated and 
ragged; amongst the last shreds of cloud, 
which have well-defined edges, diffuse or 
distinctly fibrous upright streaks become 
visible. If the whole cloud is subsiding, or if a 
wind shear has displaced the dissolving tower 
to one side of it, then the streaks are left in 
clear air, and may themselves survive only 
another minute or two before evaporating. 
No shower develops, although a brief, faint 
radar echo may have been detected: evidently 
the streaks consist of particles much larger 
than the cloud droplets, and for this reason 
sink rather quickly and are noticeably more 
reluctant to evaporate. It is remarkable that 
cloud towers do not usually shrink steadily as 
they dissolve; rather, after they have become 
emaciated, a stage is suddenly reached at 
which only a few ragged fragments exist 
here and there throughout the volume pre- 
viously occupied by the tower. Sometimes 
these frayed residues may have a rather fibrous 
arrangement, but in a clear atmosphere an 
experienced observer rarely has difficulty in 
distinguishing this appearance from the more 
diffuse and persistent streakiness which is the 
only trace of the transformation called glacia- 
tion. 

The fibrous texture characteristic of ice 
clouds is attributed to the large size of their 
particles, which fall in trails from localised 
condensation-regions. The crystals can be 
regarded as precipitation-elements, which can 
fall hundreds of metres or even some kilo- 
metres from their birth-place before evapo- 
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rating (LUDLAM, 1948; 1956). Their large size 
(linear dimension about 200 microns) is 
attained because of their very small initial 
concentration (about 1/1000 of the concen- 
tration of droplets characteristic of liquid 
clouds), and because, forming only at approxi- 
mate saturation with respect to liquid water, 
their fall commonly takes them through a 
layer a kilometre deep or more, which although 
cloudless is supersaturated with respect to ice, 
so that an abundant store of vapour is available 
for condensation. Consequently the particles 
of small ice clouds can by condensation alone 
reach the size of precipitation elements, 
whereas droplets can reach this size only by 
aggregation processes inside rather thick 
clouds. Nevertheless, the trails of raindrops 
which fall from the localised regions of 
copious condensation inside cumulus also 
have a fibrous structure, although their 
nearness to a ground observer and poor 
illumination may make this coarser than seen 
in the high ice clouds, and it is clear that the 
fibrous structure does not signal the nature of 
the particles, but only shows that they have 
the size of precipitation elements (dimen- 
sion 200 microns or more). The transformation 
of cumulus tops into fibrous trails often occurs? 
in those tropical shower clouds which are 
wholly below the level of the 0° C isotherm, 
and is therefore not appropriately called a 
“glaciation”. It seems equally inadvisable to 
use this term to describe the traces of streak 
formation which occur in cumulus tops which 
are supercooled, but apparently insufficiently 
for the development of an undoubted glacia- 
tion. We therefore propose that thisname should 
be reserved for the transformation of cumulus 
summits into a large and persistent ice cloud, 
while the more evanescent streakiness which 
occurs in decaying towers should be called a 
‘fibrillation’. 

One of the most common and striking 
features of glaciation in a large cloud is the 
cessation of sinking motions on the sides of 
the cloud summits. The nodules which con- 
stantly emerge on the upper cap of a growing 
cumulus tower effect the rapid mixing of 
the cloud air with the clear environment; 
this mixing is accompanied by a chilling as 
the cloud particles evaporate, and thus by a 


1 As seen in films of J. S. Markus and V. J. SCHAEFER. 
See also, for example, SCHAEFER, 1949, p. 7. 
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loss of buoyancy.” The nodules can be seen 
to sink aside from the cap of the cloud, which 
grows upwards by the protrusion of its in- 
terior. The sinking motions on the cloud 
sides are particularly conspicuous on time- 
lapse films. At temperatures considerably 
below 0°C the humidity in the surrounding 
clear air may be well below 100% with 
respect to liquid water, but may correspond 
to near-saturation with respect to ice. Eva- 
poration during mixing processes at the 
boundaries of a glaciated cloud then becomes 
less intense, the cloud edges become less 
sharply defined, and the chilling becomes 
insufficient to cause noticeable sinking motions 
and to promote evaporation. Instead, the 
residues of successive glaciating towers slowly 
accumulate and produce the characteristic 
anvil cloud, so that the typical shape of 
cumulus, which taper upward, gives place to 
the typical shape of mature cumulonimbus, 
which taper towards the middle levels. On 
the other hand, the streaks which are observed 
in fibrillated cumulus summits are essentially 
evanescent; they are observed to sink rapidly, 
dissolving into clear air, disappearing behind 
fresh cloud towers or settling into the bulk 
of the cloud. 

Evidently fibrillation marks a critical stage 
of shower formation, occurring in clouds 
which have been barely capable of producing 
particles of precipitation-element size. This 
transformation can be related to the maximum 
height, or dimensions, attained by the parent 
cloud tower. In contrast, it is hardly possible 
to define the exact stage at which the glaciation 
of a large shower cloud occurs, because this is a 
transformation which occurs more slowly, 
and which spreads sideways to later towers, 
and downwards, so that eventually the whole 
cloud, and the precipitation from it, may 
acquire the fibrous texture. In the observa- 
tions to be described, particular attention was 
paid to determining the stage of cumulus 
development in which fibrillation occurred. 


3. Observational techniques 


The observations were made during the 
summer months of 1954 and 1955 near Oster- 
sund, in the province of Jamtland, central 
Sweden (Fig. 8). The field-station is a few 
kilometres north of Östersund, on a hill 
near the east shore of the lake Storsjün. On 
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Figs. 1—7. Fibrillation and glaciation in cumulus towers. On the left and extreme right of Fig. 1 are disintegrating 

towers. That on the left has fibrillated and has left diffuse streaks near its peak level, above shreds of droplet cloud 

which evaporate more quickly. On figs. 2 and 3 traces of diffuse streaks can be seen beneath the subsiding 

remnants of the right-hand tower (‘marginal’ fibrillation). The central cloud towers build rapidly to much greater 

heights, and in figs. 5 to 7 acquire the soft texture associated with glaciation. The interval between successive pictures 
is about 2 min. 
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Oviksfjällen 


l Anarisfjällen 


River 


Fig. 8. Location and surroundings of observing station at As, Sweden. The theodolite base-line, just over 2 km 
long, extends between Torsta and Landsom. The treeline (about 600 m above sea level) is marked, and peak 
heights are given in metres. 


days of scattered cumulus the air over the cool 
waters of this lake is usually cloud-free, so 
that we have unobstructed views of the clouds 
which develop over the mountainous areas 
which lie farther west. The nearest of these 
are the Ovik and Anaris, which have a number 
of peaks above 1,000 m. Over the whole 
eastern hemisphere the ground undulates 
between 300 and 400 m above sea-level; it is 
mostly covered by forest, with numerous 
small rivers and lakes. 

The air is very clean, and on bright days it is 
possible from the ground to survey and meas- 
ure cumulus clouds at distances of 100 km 


or more. The measurements have been made 
by observers equipped with pilot-balloon 
theodolites and in telephonic communication, 
at either end of a N-S base-line just over 
2 km long. The technique of observing has 
already been described in a preliminary report 
(LupLAM, 1955). The observations made from 
the ground are supplemented by aircraft 
soundings, under the direction of ground 
observers, and by routine and special radio- 
sonde ascents from the nearby station in Oster- 
sund. Occasional use has also been made of 
3 cm radar, but simply for the detection of 
showers when the clouds are scattered the 
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eye is superior. In very cloudy weather, and 
when the visibility is comparatively poor, 
this is of course no longer true. On one or 
two occasions when the visibility decreased 
to about so km the observational work was 


hampered. 


4. Accuracy of the observations 


After some practice the observers are able 
to identify and train their theodolites upon a 
particular detail of a cloud tower some 20 to 
30 seconds after one of them has chosen it for 
measurement. This detail is usually the top 
of one of the nodules which continually 
emerge from the cap of the tower, or some 
other detail of size less than 10 m, which at a 
distance of about so km subtends an angle of 
less than o.o1°. Its azimuth and elevation are 
measured simultaneously with each theodolite 
on a signal from one observer. The tower is 
kept under observation; the measurement 
may be repeated, the rate of rise of the tower 
may be measured, and the peak height reached 
by the tower before its dissolution begins is 
recorded. Each theodolite! is orientated by 
sighting the other, and can be read to 0.01°; 
after each few observations the theodolites 
are levelled and freshly orientated: the bearing 
of the other instrument is rarely found to 
have wandered by more than 0.01° (usually 
in the same sense at each theodolite). 

The observations are given to a computer 
who uses a 20-inch slide rule to evaluate the 
distance (d) and the height (h) of the cloud 
detail with respect to the control station, from 
the equations 


d= csin À /sin (¢ —A) (1) 


where c is the length of the base-line (2,154 m), 
and ¢ and A are the azimuths from the control 
and second station, respectively; and 


h=d tan © +6.8 x 1078 d? (2) 


where © is the elevation from the control 
station; h and d are expressed in metres. The 
last term in equation (2) is a correction-term 
to take into account the curvature of the earth 
and the refraction of light; it was used after 
it had been found to give correspondence 


1 The theodolites used were British Meteorological 
Office balloon theodolites, made by E. R. Watts and 
Sons. 
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within a few metres with the paths of rays 
calculated to have an average refraction. 
The principal errors which arise in the cal- 
culated distance of a cumulus detail are due to 
errors in the angular measurements of the 
azimuths. If these amount to d¢, 6A, it may 
readily be shown that the error in the distance 


d is given by 
öd= -d (66-64) /sin (6-A) (3) 


The algebraic difference of the azimuth errors 
is not likely to exceed 0.02°; using this value 
in equation (3), the percentage error in the 
calculated distance of the cloud can be derived 
as a function of the true distance and azimuth, 
and is shown in Fig. 9. Two-theodolite meas- 
urements are not usually attempted on clouds 
lying within 20° of the direction of the base- 
line. Fig. 9 also shows the resultant error in 
measuring the height of a cloud detail 5,000 
m above the control station, assuming an 
additive error of 0.02° in observing its angular 
elevation. It can be seen from the diagram 
that, since most of the observations are made 
on clouds at ranges between 25 and 75 km, 
the errors in the height determinations are 
not likely to exceed 100 m. 

It is not always possible for both observers 
to identify a detail in the cloud bases, and 
their height is therefore estimated by noting 
the elevation of the cloud base each time a 
tower is measured, and assuming it to have 
the same range. The heights deduced in this 
way are probably correct within so m, and 
agree well with measurements made by air- 
craft. There is usually a marked diurnal varia- 
tion in the height of cloud base, and if this is 
taken into account it appears that on most 
days the cloud base at a particular time is 
uniform within about 100 m, although a 
consistent difference exceeding this may occur 
over certain sectors of the sky, or even in the 
same locality when the wind is light or when 
there are showers. When therefore the gener- 
al level of the cloud base is known, a useful 
estimate of the height of a cloud top can be 
made with one theodolite simply by measuring 
the angular elevations of the top and base 
of the cloud. The error of such a determination 
is not likely to exceed about 10 %, i.e., $00 m 
at 5 km. In this way, and more accurately by 
using radar to determine the range of clouds 
containing precipitation, or aircraft observa- 
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Fig. 9. Maximum likely errors in distance and height of cloud tops derived from theodolite observations, as a 

function of distance and angular bearing away from the direction of the base-line. An error of 0.02° is assumed 

in the difference of the azimuths observed from the ends of the base-line, and an additive equal error in the 
angular elevation measured at the control station. 


tions, measurements may be made on clouds 
obscured from one observer and on clouds 
lying in directions close to that of the base-line. 

The theodolite observations are supple- 
mented by 16 mm time-lapse colour films 
and sequences of photographs on 6x9 cm 
films, which form an invaluable record. Ex- 
perience of making measurements from these 
pictures leads us to suppose that by using a 
large negative size and taking a number of 
precautions, the analysis of photographs taken 
from either end of a similar base-line could 
provide height determinations with errors 
not more than three times as great as those 
arising in the 2-theodolite measurements. 
Although this photographic method has the 
advantage of providing a permanent record, 
and possibly more information at less personal 
trouble, its reduction has to be postponed until 
the film has been processed and is more te- 
dious: the practised theodolite observers and 
their computer obtain a more accurate range 
and height within a minute of selecting a 
cloud. When this rapidity has been achieved, 
between so and 100 clouds can be measured 


in a day, and the number of observations is 
more often limited by the number of visible, 
interesting clouds than by the stamina of the 
observers. Their measurements, supported by 
the film records, are probably more satis- 
factory, and certainly cheaper, than those 
obtained primarily by the photographic 
method. 


5. The fibrillation level 


Observations on a number of days of scatter- 
ed cumulus are illustrated in Figs. 10 to 14. 
These diagrams show the peak heights reached 
by individual cumulus towers, and whether 
or not fibrillation or glaciation of the towers 
occurred. When only traces of fibrillation 
were observed the occurrence was noted 
(before the height of the tower was computed) 
as a “marginal” transformation. The diagrams 
include estimates of the height of the cloud 
bases, from which the considerable diurnal 
variation is apparent. 

These days are typical of those on which 
the depth of the cumulus exceeds about 2,500m: 
it is then possible to recognise a layer within 
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Fig. 10. Observations of cumulus bases and tops on several days in 1955. 


a) peak heights reached by ascending towers: an open circle indicates that no 
fibrillation was observed; a black dot indicates that fibrillation or glaciation 
occurred. If only traces of fibrillation were detected the letter ‘m’ is added; 
if a fibrillated tower was observed only after it had subsided an unknown 
but small distance from its peak height, an upward-pointing arrow is added. 


b) cloud bases: a large open circle indicates the convective condensation 
Ss 1% 
level of air about 2 m above the ground. Shower symbols are placed 
beside points representing bases measured near showers, which are often 
below the general base-level. Distinction is sometimes made between the 
general level of the cumulus bases in the mountainous (SW) sector and 
in the flat terrain near and east of the field station (NE). 
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which some of the cumulus tops become fi- 
brillated; usually this layer, the fibrillation zone, 
is not more than a few hundred metres thick, 
and contains most of the observations of 
Tellus VIII (1956), 4 
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marginal transformations, and 
all clouds which rise above the 
upper limit of the layer show 
marked fibrillation or glacia- 
tion. Quite commonly the zone 
is apparently only one or two 
hundred metres deep, and it 
is permissible to speak of the 
fibrillation level. On most days 
the cumulus have not become 
large enough for their tops to 
become fibrillated until about 
noon, after which time the var- 
iation in the height of cloud 
bases is rather small, so that, 
although we shall relate the 
fibrillation to the height of 
the towers above the bases, we 
have not usually tried to find 
a diurnal trend in the fibrilla- 
tion level. Sometimes a varia- 
tion has seemed to be present 
and it has been desirable to 
distinguish between conditions 
in the morning and afternoon, 
but mostly the observations are 
hardly sufficient to justify this. 

That it should be possible to 
recognise a fibrillation level, 


soc and that it should be found at 


practically the same height on 
successive similar days (for ex- 
ample, 10 to 12 July, 1955), 
gives us confidence in the ac- 
curacy of our observing meth- 
ods. It is still more encouraging 
to discover such evidence of 
orderly behaviour in popula- 
tions of cumulus whose indi- 
viduals superficially have a very 
varied character, growing over 
a variety of terrain and attaining 
a range of sizes and durations. 

The height and temperature 
of the fibrillation level, and its 
height above the cloud bases, 
vary considerably according to 
the properties of the air mass, 
as shown in Table 1, which 


includes all days during two summer periods 
when a fibrillation level could be determined. 
The height of the fibrillation level varies over a 
range from 4,000 to 6,900 m, and its tempera- 
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ture from —9 to — 23°C, while 
its height above the cloud base 
varies from 2,500 to 4,600 m. 

On days when the depth of 
cumulus is restricted to about 
2,500 m or less it is not usually 
possible to observe fibrillation 
in the cloud summits. Obser- 
vation is frequently hindered 
because the individual towers 
do not rise far out of the gen- 
eral cloud mass, having a small- 
er size and upward speed than 
is usual in larger clouds. Often 
they have entered a very stable 
layer or inversion which is sup- 
pressing the convection, have 
“overshot” their equilibrium 
level and acquired a negative 
buoyancy, which causes them 
to subside and disappear again 
from sight before disintegrat- 
ing, unlike towers below the 
extreme summit-levels on days 
of much larger clouds, which 
are brought to rest by erosive 
mixing after emerging rapidly 
from the parent clouds. Fre- 
quently, too, the towers of the 
small clouds spread out at the 
base of the stable layer, pro- 
ducing expanding shelves of 
stratocumulus cumulogenitus 
which obscure succeeding tow- 
ers. 

When, therefore, showers 
develop in such small clouds, 
it is not possible to relate their 
formation to a fibrillation level, 
but only to refer to the level 
which the summits attain. Usu- 
ally this is only a few hundred 


metres above the peak heights HEIGHT, km. © 


reached by the towers of still 
smaller clouds which fail to 
produce showers. Table 2 then, 
which summarises our observa- 
tions on days of small shower 
clouds, can be regarded as an 
extension of Table 1, and shows 
that showers have occurred in 
clouds having a depth as little 
as 900 m, and summit-level 
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Fig. 14. See Fig. ro. 
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Table 1. Observed values of height and clear-air 

temperature of the fibrillation level, or of the 

average height of the fibrillation zone, and its 
height above the cloud base. 


: Height 
Height | Temper- = 
Date above S.L. are above 
°C cloud base 
m m 

1954 | 
Sr AUB Sere 4,800 — 13 3,500 
TRNAS Pe. rn 4,000 — 9 2,500 
DA AUS tens Goes 4,000 — II 2,500 
TS SED reger 5,200 — 23 4,000 

1955 

28 Jun (3)! 5,000 — 19 2,600 
OMI A. 5,300 — 16 3,100 
ie jC Gy) Somes 5,000 — 13 3,500 
TOUT (0). ee. 6,800 — 19 4,300 
RING 6,800 — 19 4,300 
2 Tule(8)e. Je 6,900 — 21 4,400 
Pal (Qe) nas, tc 6,000 — 13 4,000 
Ee pl (TO)... 6,600 — 19 4,600 
8 Aug (Ir) 4,300 — 9 3,300 
TOMA ER ((T2)G 3: 4,800 — 10 3,300 
PS AUS (T3)... 5,600 — 15 3,300 
116 Aug (14). 5,400 — 12 3,300 


1 The numbers in brackets identify the occasion of 
the curves of Fig. 17 and the points of Fig. 19. Two 
of these occasions are listed in Table 2. 


Table 2. Observed values of height and clear-air 

temperature at the tops of cumulus which be- 

came small shower clouds, and their height above 
cloud base. 


Height | Temper- A 
Date aboveS.L.| ature SS 
cloud base 
m BK m 
1954 
19 Aug a.m. (I). 3,900 — 10 2,400 
p-m. (2). 3,800 — 9 2,000 
2OMAUSIA MN... 3,200 — 4 1,700 
Pen. 3,300 — 4 1,400 
BIN ee 2,300 — 3 900 
SY Site Sooo nen 2,400 — 2 1,500 
LAISSES MONO 2,900 — 7 1,600 
1955 
PO | (DORE cho oe 3,900 — 9 2,300 
2 Aug. 2,700 — 4 I,000 
2 INT eG 0 ao 2,300 + 3 1,200 


(On some of these occasions the cumulus tops 
penetrated a few hundred m above an inversion and 
the temperatures within the tops may conceivably 
have been 3 or 4° C lower than the clear-air tem- 
perature at the summit level) 
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temperatures as high as +3°C. Together, 
these tables illustrate the great variation in 
the sizes and summit temperatures at which 
showers form, even in one place during the 
same season, and doubtless the range could 
be increased further by including observations 
from other regions. 

In 1954 it was noticed, as previously re- 
ported (LupLAM, 1955), that the height of 
the fibrillation level above the cloud base 
seemed to be related to the vigour of the 
convection, being great when the cloud 
towers ascended rapidly, and small in the little 
clouds of shallow polar air masses, which 
probably contain only weak updraughts. This 
circumstance, and the remarkably high tem- 
peratures at the lower of the fibrillation levels, 
suggest that it is the coalescence process, 
rather than the growth of ice particles, which 
produces the precipitation elements in fi- 
brillating cloud tops, and hence which is the 
predominant shower-forming agency. It is 
therefore necessary to discuss the course of 
the coalescence process in cumulus clouds. 


6. The coalescence process in cumulus 


The coalescence of cloud droplets in cumulus 
has been considered in a previous paper (Lup- 
LAM, 1951). Here the discussion was based 
upon two premises: 

1. In the lowest few hundred metres of the 
cloud there are some droplets of radius 20 to 
40 4, in concentrations of order 100/m%, 
amongst a droplet population having uniform 
radii of several microns. 

2. The droplets are lifted through the cloud by 
convective bubbles, in which the concentration 
of condensed water is that corresponding to 
adiabatic rise from the cloud base, until the 
bubbles arrive at the cloud summits and emerge 
as individual towers. There the mixing with 
clear air causes a rapid evaporation of all 
small droplets. 

It was shown that if, during the rise of a 
bubble to the summit-level, the large droplets 
have attained radii of about 150 u, then their 
further growth into raindrops and the for- 
mation of a shower is practically inevitable. 
This is because droplets of such a size have a 
fall speed exceeding 1 m/sec, so that they 
tend to have settled some hundreds of metres 
through the bubble and remain poised in the 
general cloud tops with a good chance of 


434 


escaping the evaporation to which the towers 
are exposed; because even if involved in the 
evaporation they are able to settle through a 
few hundred metres of unsaturated air to 
re-enter the cloud mass with only slight 
shrinkage, and because in a moderate concen- 
tration of cloud water further growth by 
coalescence beyond this radius to raindrop 
size is very rapid, occupying only a few min- 
utes. 

Equally it may be considered that the pres- 
ence of droplets of about this size and concen- 
tration is necessary for the occurence of fi- 
brillation in a dissolving cloud tower, so that 
the summit-levels calculated in the previous 
paper (as a function of cloud base temperature, 
which governs liquid water content within 
the cloud, and of rate of rise of bubbles) to be 
necessary for the formation of showers, may 
be regarded as theoretical fibrillation levels. 

However, these premises need some recon- 
sideration, and are discussed in the following 
paragraphs. 


7. The size spectrum of droplets in cumulus 


In the previous paper it was shown that 
droplets of the size assumed to occur in 
rather small concentrations near the cloud 
base and efficiently initiate the coalescence 
process, are readily introduced into the clouds 
of maritime air masses by the giant hygro- 
scopic nuclei, which there have concentrations 
of order 10°/m?. These nuclei may be absent, 
or present in much smaller concentrations, 
in continental air masses, and even in maritime 
air masses if the winds are too light to produce 
sea-spray. It was therefore suggested that in 
these circumstances the coalescence process 
may be hindered. 

Recently observations have been made of 
the size spectra of droplets in the lower parts 
of cumulus in Russia (ZAYTSEV, 1950), in 
New Jersey and Florida (WEICKMANN and 
AUFM Kamps, 1953), and in England (Mur- 
GATROYD, 1955). These writers apparently 
agree that small cumulus always contain 
abundant concentrations of droplets having 
radii of about 30 u. 

According to Zaytsev droplets of radii up 
to 25 u are encountered 200 m above the 
base of small clouds, and “some droplets of 
the order of 100 u diameter are found, on 
the average, from 350 m above the cloud 
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base”. The diagrams of Weickmann and aufm 
Kampe show that most of their observations 
in cumulus humilis were made at soo or 900 
m above the cloud bases, and their summary 
of the droplet size spectrum typical of these 
clouds shows the presence of droplets of radius 
about 30 u in concentrations of order I0-1/cm?. 
Finally, Murgatroyd states that well inside 
cumulus at heights of more than 300 m above 
the level of the bases the mean droplet radius is 
about 10 u, and that there are a few large 
droplets of radius 40—50 u above 300 m, 
and of radius exceeding 100 u above 1,000 m., 

Remarkably enough, none of these authors 
specifies the volumes of the samples of cloud 
air from which their droplet spectra were 
obtained, but since all of the measurements 
were made by microphotography of small 
slides exposed during flight, it can be assumed 
that their spectra refer to volumes of up to 
100 cm?. Naturally if much smaller volumes 
had been sampled, there would have been no 
mention of the large droplets, and equally 
it may be speculated that if much larger 
volumes are examined, smaller concentrations 
of even bigger droplets will be discovered. 
Clearly it is important to relate a given 
spectrum to the volume which it represents. 

However, on the basis of the measurements 
quoted, it appears that at heights about 300 m 
above the bases of most small cumulus, droplets 
of radius about 30 wu exist in concentrations 
between 10°? and 1071/cm®. It is most unlikely 
that these large droplets could arise by conden- 
sation alone on giant hygroscopic nuclei. 
According to Woopcock (1953), the typical 
concentrations of giant nuclei near cloud base 
in air over the oceans with fresh surface winds 
(Beaufort force 4—5) are those shown in the 
second column of Table 3. The last columns 
show the radii which are attained by these 


Table 3. Growth of giant sea-salt nuclei by 
condensation. 


Radius in 
cloud after 


Typical | Radius at 
Mass of salt] concentra-| cloud base 


in nucleus tion (Keith and 
(Woodcock)| Arons) |300 sec} ı hr 
g cms u u u 
toi) Ra oe Io 14 22 
10-9 10m 18 2 35 
307" 10> 35 4 56 
4 
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nuclei at the level of the cloud base if they 
are lifted from low levels at speeds of about 
I m/sec, according to KEITH and Arons 
(1954). The fourth column shows the radii 
attained at a height of 300 m above the cloud 
base if rise continues at this rate, calculated as 
described previously (LUDLAM, 1951, p. 404), 
but taking into account the heat economy of 
the droplets. 

A large cloud may commonly have existed 
for about an hour before reaching a size at 
which it is likely to produce a shower, and 
it is possible that at the end of such a period a 
bubble rising through the lowest parts of the 
cloud becomes infected with droplets which 
have survived the whole of this time. The last 
column of the table therefore shows the sizes 
of the giant nuclei after they have been 
suspended an hour inside the cloud. It can be 
seen that even if they infected fresh bubbles in 
undiminished concentration, and even if this 
concentration is that found in oceanic air, 
then the large droplets grown upon the 
nuclei would still be about 10 times fewer 
than those actually encountered. 

It is therefore evident that the observed 
large droplets are grown predominantly by 
the coalescence process, which is presumably 
active in small clouds even though it does not 
succeed in producing precipitation. Recent 
developments in the bubble theory of con- 
vection (Scorer and LUDLAM, unpublished) 
recognise that during the rise of a buoyant 
bubble there is an intense mixing with the 
environment through the cap of the bubble, 
while there are strong stirring motions within 
the bubble itself. It must therefore be supposed 
that a bubble which rises through a cloud 
which has been in existence for some time 
incorporates numbers of the particles of 
preceding bubbles, in ever-decreasing concen- 
tration the greater their age. In this way a 
small proportion of the cloud particles are 
able by prolonged condensation and coales- 
cence to attain an abnormal size, and the size 
spectra of the droplets become very broad: 
it is impossible otherwise to account for the 
breadth of the observed spectra, for the rise of 
air unmixed from the level of the cloud base 
produces very narrow spectra containing 
droplets of practically uniform size (see, for 
example, HOwELL, 1949). According to this 
view, even in small clouds the size spectra 
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may include radii much greater than 30 u if 
the sampling volume is sufficiently increased. 

In this connection we recall previous obser- 
vations in England (Browne, Day and Lup- 
LAM, 1955), in which droplets of radius about 
150 u were detected in concentrations of order 
1/m$ inside cumulus of depth 1,500 m: they 
were found on three successive days at heights 
between about 600 m and 1,200 m above 
the cloudbases, in clouds which did not 
produce showers. We believe that droplets of 
this radius and concentration commonly 
exist in cumulus having this size, and that 
this accounts for the readiness with which 
such clouds produce echoes from suitable 
radar (see, e.g., PLANK, ATLAS and PAULSEN, 
1955). The size spectra which have been 
detailed in recent literature must be regarded 
as extending continuously to radii of at least 
150 u in small cumulus, and of course to 
values appropriate to raindrops in larger 
clouds which produce showers. 


8. The concentration of liquid water in 
cumulus 


It has previously been emphasised (LUDLAM 
and Scorer, 1953) that only within actively 
rising bubbles, which occupy only a fraction 
of the total volume of a large cloud, can it be 
expected that the concentration of condensed 
water approximates to the value corresponding 
to adiabatic rise from the cloud base. The 
rate of dilution of a rising bubble is very 
considerable, and the departure from the 
saturated adiabatic process within one bubble 
must greatly depend upon how soon it 
follows in the path of another, that is, upon 
the rate of bubble production. Only if this 
rate is increased beyond a certain value can 
a whole cloud grow, with successive bubbles 
retaining a greater buoyancy during their 
rise and hence reaching greater heights. Prob- 
ably only during the building phase of a 
large cloud do the conditions within rising 
bubbles closely approach values appropriate 
to the rise of air unmixed from the cloud base. 
If the rate of bubble production falls, then the 
cloud ceases to grow or even shrinks, and 
although bubbles may still emerge as towers 
in the summits, there must be a stronger de- 
parture from adiabatic conditions at all levels 
in the cloud. 

Perhaps for this reason, and because of the 
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difficulty of selecting for examination clouds 
in the building phase, and then making trav- 
erses through the internal bubbles, only 
rarely have flight measurements with im- 
proved instruments indicated liquid water 
concentration close to the adiabatic values 
(WARNER, 1955, MURGATROYD, 1955). 


9. Calculation of the course of the coales- 
cense process in cumulus 


In view of the conclusions of the two pre- 
vious sections, it must be questioned whether 
the considerations mentioned in section 6 form 
a satisfactory basis on which to estimate the 
course of the coalescence process. Every cloud 
which approaches the size associated with 
shower formation has been in existence for 
some time and probably already contains very 
small concentrations of drops of nearly 
the size of raindrops, and almost certainly 
larger than the radius of 150 u previously 
specified as a critical one for shower formation. 

Nevertheless, it is observed that when the 
cloud is large (depth about 2,500 m or more), 
during a building phase successive bubbles 
form separate towers in which the production 
of fibrillation, of a detectable radar echo, 
or even of a shower, occurs rather suddenly, 
suggesting that within a particular bubble the 
coalescence process has progressed beyond a 
well-defined threshold. Evidently at this stage 
the concentrations of drops of radius greater 
than about 150 w have suddenly increased 
from insignificant values of 1/m% or less to 
values exceeding about roo/m?, and it will 
be shown below that this is readily inter- 
pretable as due to the growth to this size, 
within the individual bubble, of those droplets 
of radius about 30 u with which it became 
infected at the beginning of its ascent through 
the cloud mass. In such a large cloud, in a 
building phase which brings it to the point 
of shower formation, it may still be reasonable 
to expect the liquid water content within 
rising bubbles to approximate to the adiabatic 
values, and altogether the previous simple 
model still seems a reasonable basis for a cal- 
culation of the course of the coalescence 
process. 

On the other hand our observations suggest 
that in small clouds, of depth less than about 
2,500 m, the showers develop more slowly, 
for it is often difficult to define a stage at 


F. H. LUDLAM AND P. M SAUNDERS 


which a shower has: formed. Usually no 
fibrillation can be seen in the summits, and 
often the first visible signs of precipitation 
are faint fallstreaks beneath spreading shelves 
of stratocumulus cumulogenitus. Gradually, 
over a period of perhaps half an hour or more, 
these become more extensive until most of 
the cloud acquires a glaciated structure, so 
that it appears as a miniature anvil cloud. 
The showers which fall from these clouds 
are very slight, and they may even fail to 
reach the ground. In such clouds the evapora- 
tion near the summits is very slow and the 
growth of the cloud particles into precipitation 
elements proceeds over a long period within 
a succession of convective bubbles. The simple 
model in which the cloud droplets within a 
bubble are rapidly evaporated when it reaches 
the summit-level of the cloud, and in which 
the period of droplet growth is regarded as 
limited to the period taken for the bubble to 
traverse the cloud from base to summits, is 
therefore quite inappropriate, for it considers 
only the history of the bubble and not the 
history of the entire cloud. 

Only during the building phase of a large, 
vigorous cloud is the model likely to be at 
all satisfactory, and even then it is necessary 
to circumvent consideration of the cloud 
history (i.e., of the early stages of the coa- 
lescense process), by appealing to the flight 
observations which show the existence of 
significant concentrations of 30 u-radius 
droplets near the bases of practically all cu- 
mulus, and which therefore can be supposed 
to be present at the beginning of a bubble 
history. 

With this reservation in mind, and remem- 
bering also that the assumption of adiabatic 
concentrations of liquid water is not then 
unreasonable but is still not supported by 
direct measurements, we can estimate the 
course of the coalescence process as a bubble 
rises from the cloud base towards the summits. 

As before (LUDLAM, 1951) we restrict atten- 
tion to the growth of an individual droplet of 
radius R and fall speed V, settling through a 
cloud of smaller droplets of uniform radius r 


and fall speed v, as described by the equations 
40dR=E,E, (V—v) wdt (4) 
where @ is the density of water, E, is the 


collision efficiency, E, the coalescence effi- 
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ciency, and w is the concentration of water 
in the form of the smaller droplets; and 


dz = (U1-V) dt (5) 


where U! is the updraught speed and z is 
height (conveniently measured above cloud 
base). 

Certain simplifications of these equations 
are reasonable. In the absence of any informa- 
tion to the contrary, the coalescence efficiency 
E, is assumed equal to 1. As before, the collec- 
tion efficiency E, is taken as the value E, 
calculated by Langmuir, adjusted to take 
some account of the finite size of the collected 
droplets. Under average conditions of air 
density and temperature, and assuming the 
collected droplets to have a uniform radius of 
10 u, the value of E, varies with R as shown in 
Table 4; these values cannot be regarded as 
very secure, particularly at the lower radii. 
Reported attempts at experimental verification 
give efficiencies both higher and lower than 
Langmuir’s values (see, for example, KINZER 
and Coss, 1956). 


Table 4. Fall-speed V and efficiency of catch E, = 

(E; 2: + r/R)?, where E, is Langmuir’s efficiency 

of catch, of a droplet of radius R falling through 

a cloud of droplets of radius r (= Io u). Air pres- 
sure 700 mb, temperature —7° C. 


R y = E, 
17 cm/sec 

21 5.5 0.31 1.08 
30 212 0.44 0.98 
50 26.7 0.55 0.89 
100 ER 0.67 0.85 
150 127 0.72 0.85 


The liquid water concentration w is assumed 
to be that corresponding to adiabatic ascent 
from the level ofthe cloud base, m. In equation 
(5), which has been used to determine the 
height of the growing droplet and hence 
the liquid water concentration in its environ- 
ment, it is preferable to replace (U! - V), 
where U! is the speed of the updraught in 
the cloud, by U, the speed of ascent of the 
bubble containing the droplet. This is per- 
missible in view of the circulatory and stirring 
motions within a bubble, and because we are 
not concerned with values of V greater than 
about 1 m/sec, which may allow the growing 
droplet to fall out of a bubble. 
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The equations can then be written 


40dR = E, (V-v)mdt (6) 
and de = Orde (7) 
whence 4dR/E, (V-v) = mdz/Uo (8) 


I 


If now E, and V are regarded solely as 
functions of R, and as m and U are functions 
of z, we obtain the following relation for the 
growth of a droplet from an initial radius Ro, 
at the cloud base, to a final radius R, at a 
height z, above the cloud base: 


21 


J sarin, (V-v)= f (m/Ue) dz (9) 


O 


The value of the integral on the left can be 
computed for various values of R, and R,, 
while that on the right can be computed 
if the height and temperature of cloud base, 
and the variation of U with z, are given. 
The hypothesis that fibrillation is due to the 
growth of droplets from radii of about 30 u 
to radii of about 150 u can be tested by in- 
serting these limits to the integral on the left 
and comparing its value with that on the right 
when the height above the cloud base of the 
fibrillation level is taken as z,. To make 
this comparison it is necessary to determine 
the variation of U with height. 


10. The speed of cloud updraughts 
Probably the best method of observing 


the speed of cumulus updraughts is to measure 
the rising speed of a glider soaring inside 
the clouds, since the pilot is able to manoeuvre 
his aircraft into the updraught cores. Estimates 
of vertical air speed can also be made during 
horizontal traverses of clouds by powered 
aircraft, but are difficult to interpret, for the 
aircraft may not have passed through the 
core of an updraught, which occupies only a 
portion of the cloud. On the other hand, 
the existence of a circulatory motion within a 
cloud bubble results in the vertical air speed 
in its interior exceeding, by a considerable 
fraction, the speed of rise of the bubble (as 
defined by the motion of its cap), which is 
the appropriate value to consider when examin- 
ing the growth of droplets within it. In view 
of these practical difficulties in exploring 
cumulus with aircraft we have tried to assess 
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the rising speeds of bubbles by observing 
from the ground the rate of ascent of the caps 
of individual cloud towers. 

The towers become visible when ascending 
bubbles emerge from the bulk of the cloud. 
In our experience, provided they are not 
glaciated, they are nearly always within one 
minute observed to be decelerating, and come 
to rest after rising a distance about equal to 
their diameter at the time of emergence 
(about 400 m). This is attributed to the strong 
effect of mixing with clear air, as opposed 
to the saturated air through which the bubble 
rises within the cloud. The evaporation which 
occurs during mixing with clear air is accom- 
panied by a buoyancy-destroying chilling, 
and the volume of air which can continue to 
rise quickly diminishes, so that the bubble is 
said to suffer a rapid “erosion” (as discussed by 
Malkus and Scorer, 1955). According to the 
bubble theory, a bubble cap rises at a speed 
approximating to a limiting speed U given 
by a relation of the form 


U? = kRgB (10) 


Here R is approximately the half-width of 
the bubble, B is the buoyancy (usually ex- 
pressed as AT/T, where T is the virtual 
temperature and /\T the excess of virtual 
temperature within the bubble, adjusted in 
the case of cloudy air to include the effect of 
the condensed water upon the air density), 
and k is a constant. The deceleration of cloud 
towers in clear air is associated with a decrease 
in R, and perhaps also in the buoyancy B, 
during mixing and erosion. 

Evidently only the rate at which a cloud 
tower emerges from a cumulus is at all likely 
to be representative of the rate of rise of 
bubbles within the cloud mass. Our attention 
was therefore concentrated on observing the 
rate of rise of cloud towers as soon as possible 
after they became identifiable. Measurements 
were made directly by theodolite and subse- 
quently also from time-lapse film records. 
Both kinds of observations provided an 
angular rate of rise, which was converted 
into a true rate of rise by assessing the distance 
and height of the towers, and, equally neces- 
sary, their horizontal motion. Usually the 
distance and height were known from the 
double-theodolite measurements; on other 
occasions they were deduced from the angular 
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height of the cloud base (whose true height 
throughout the day was observed) or ob- 
tained from the range of associated radar 
echoes. The horizontal motion of the cloud 
towers was assumed to be that of the wind at 
their level, which was known from obser- 
vations of the radio-sonde and special pilot 
balloons, and from theodolite observations 
of the displacement of measured clouds. On 
occasions of‘ considerable wind shear the 
horizontal motion of a cloud tower may differ 
from the wind at its level by several m/sec 
(Scorer and LUDLAM, 1953); however, 
the wind shear within a layer occupied by 
vigorous convective clouds is usually rather 
small, and the effect of horizontal motion was 
kept as little as possible by observing clouds 
mainly on azimuths at right angles to the 
wind direction. Consequently we think that 
the errors made in deducing the true rates of 
rise of the cloud towers are not likely to 
exceed 10 %. 

When the rates of rise of emergent towers 
are plotted as a function of height, as in 
Fig. 1$—16, a cloud of points is obtained. 
We have drawn an envelope to this cloud 
to indicate the maximum speeds observed, 
and assume that this curve represents the 
maximum rates of rise of bubbles within 
building clouds of all sizes, as a function of 
height. Where the observed values are smaller 
we suppose this to imply that the observations 
were made after some deceleration had already 
occurred, or that the bubbles had a size or 
buoyancy less than the maximum appropriate 
to their level. Strictly, it would have been 
better to take into account the varying size 
of the observed towers. However, in this 
exploratory study it was supposed that the 
smaller bubbles had already been partially 
eroded or had originated at levels well above 
the cloud base; in either circumstance they 
were likely to contain significantly less con- 
densed water and therefore be less likely to 
produce precipitation elements than the bub- 
bles having the maximum upward speeds. 
From a few observations, also, it seemed that 
the introduction of a tower width was not 
certain to reduce greatly the scatter of the 
points on the height-speed diagrams. 

Fig. 17 shows the maximum rate of rise 
of cloud towers, as a function of height, 
from near the cloud base to the fibrillation 
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Figs. 15, 16. Rates of rise of emergent cumulus towers, as a function of height, on two days in 1955. Vertical lines 

show mean upward speeds between levels corresponding to the ends of the lines, as measured from time-lapse 

films. Circles represent emergent speeds over intervals of about 200 m, as measured by theodolite. The full lines 

show speeds deduced from soundings, on the basis of bubble theory. The pecked curves, drawn as an envelope 

to the points representing individual towers, are assumed to represent the maximum rates of rise of cloud bubbles 
inside building clouds of all sizes, as functions of height. 


level, observed on one day in 1954 and on 
12 days in 1955. We assume these curves 
define the rates of rise of bubbles inside 
building clouds, and use them to evaluate the 
integral on the right of equation (9). 


11. The coalescence process and the produc- 
tion of a fibrillation level 


Two of these curves are reproduced on 
figure 18, on which the ordinates now re- 
present height above cloud base. One curve 
represents conditions on a day of vigorous 
convection (12 July, 1955) when rising speeds 
were as great as 6 m/sec; the other is for the 
afternoon of 19 August 1954, when the speeds 
barely exceeded 2 m/sec. We have not meas- 
ured smaller rising speeds, these being char- 
acteristic of days of small clouds on which 
the difficulties of observation are greater. 
Tellus VIII (1956), 4 


The diagram also shows the value of the 
expression S =f (m/ Uo) dz, as a function of 


O 
the height z above the cloud base, on these 
two occasions. In so far as adiabatic values of 
liquid water concentration are actually attained 
within rising bubbles, and in so far as the inte- 
gral on the left of equation (0) has been 
correctly obtained, the value of the integral 
S at any level corresponds to a particular 
radius of a growing droplet lifted from near 
the base of the cloud by the bubbles. The 
value of this radius depends upon the initial 
radius, and scales of radii appropriate to initial 
radii of 20 u and of 30 are entered along 
the upper margin of the diagram. Thus at the 
level where the value of the expression S 
reaches 19.6 x 1074 sec, a droplet of initial 
radius 20 u will have attained a radius of 
150 u. It can be seen that in the lowest few 
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Fig. 17. Maximum rates of rise of cloud towers on the 

14 occasions numbered in Table 1, as a function of height 

above sea level, from near cloud base to the observed 
fibrillation level. 


hundred metres of the clouds, where the 
liquid water content is small, the contribution 
to the value of the integral S is small, partic- 
ularly when the updraughts are strong. 
Thus the exact level near the base at which 
the growing droplets are introduced into a 
rising bubble hardly affects the level at which 
they attain radii of 150 u, or any larger size. 

We assume that inside a rising bubble at 
any level the drop-size corresponding to con- 
centrations significant for shower production 
is determined by the size attained by the 
droplets of radius about 30 u which were 
present in these concentrations within a few 
hundred metres of the cloud base and have 
been lifted in the bubble. This does not 
preclude the existence already in the cloud 
at this level of smaller concentrations of 
drops of similar size, which have had a longer 
history. A very important feature is the rapid 
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Fig. 18. Maximum rates of rise of cloud towers as a func- 
tion of height above cloud base on 19 August 1954 (Us) 
and on 12 July 1955 (Us), and the variation with height 
of the values of the integral S = J (m/Uo)dz on these 
two days. The values at levels 300 m below the fibrillation 
levels are marked by dots, and the curves are extended 
as pecked lines to the fibrillation levels. The other dots 
show corresponding values on other observation occa- 
sions. The scales on the upper margin show the sizes 
of a droplet growing from an initial radius of 204 or 
30 u, which correspond to the values of the integral S. 


droplet growth which occurs in the bubble 
as it approaches the cloud summits. Supposing, 
for example, initial radii of only 20 w, then 
on 12 July 1955 the significant droplet radius 
at a level 3,400 m above the base is only so u; 
at 3,900 m it is 100 u and less than 200 m 
above it attains the critical value of 150 u, 
while in a further 400 m of ascent the droplets 
have achieved the size of raindrops (radius 
soo 4). This implies that in a cloud with 
summits about 3,900 m above the base there 
are no droplets in concentrations sufficient 
for shower production of greater radius than 
100 #, which is too small to allow them to 
escape evaporation in the summit towers. 
In a cloud only 200 m taller these droplets 
attain radii of 150 4, probably cause fibrillation, 
and are likely to settle into the cloud bulk 
during the dissolution of the towers and grow 
into raindrops at lower levels. If the cloud 
tops rise a further 400 m the towers them- 
selves become charged with raindrops. It is 
therefore understandable that there should be a 
well-defined fibrillation level, and that quite 
intense showers should form rather suddenly 
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in clouds which grow only slightly beyond a 
Critical size. 

In Fig. 18 we have marked the curves 
which show the trend of the integral S with 
a dot at a level 300 m below the fibrillation 
level, and extended the curve upwards as a 
pecked line. This is because during the last 
three hundred metres or so of the rise of a tower 
its erosion is very pronounced, so that in all 
probability its water content diminishes rap- 
idly and there can be practically no growth 
of the large droplets in this stage. As a tower 
dissolves the fibrillation is usually discerned 
about this distance below the peak height 
reached by the tower, and so we regard the 
value of the integral S at a level of 300 m 
below the fibrillation level as the value appro- 
priate to the growth of drops of 150 w radius. 

The remaining 12 dots show the values 
attained at this distance below the fibrillation 
level on other days for which the calculation 
could be made. It will be observed that on 
the days when the fibrillation level is well 
above the cloud base the integral Shas attained 
values of about 20 x 10”? sec, corresponding 
to the growth by coalescence of droplets of 
initial radius about 20 w. There is a noticeable 
tendency for the value of the integral corre- 
sponding to fibrillation to decrease as the 
height of the fibrillation level above the cloud 
base falls, and in fact for the afternoon of 
19 August 1954, on which the integral has 
its lowest value, the fibrillation level implied 
is actually the minimum height of the peaks 
of clouds which produced showers, fibrillation 
not being observed. 

In view of the impossibility of defining an 
initial droplet radius, and the uncertainty in 
the evaluation of the integrals in equation (9), 
the diagram can only indicate that there is a 
relation between the size of a cloud which is 
necessary for shower formation and the speed 
of its updraughts, and that the relation is 
consistent with what is known of the coales- 
cence process. It is evident that the process 1s 
more efficient in small clouds having weak 
updraughts, probably because the individual 
clouds have long durations and are often 
protected from rapid evaporation at the 
summits by the formation of stratocumulus 
cumulogenitus shelves beneath stable layers 
limiting the cloud growth. No satisfactory 
models for estimating the course of the coa- 
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lescence process in these small clouds and of 
its initial stages in large clouds have yet been 
evolved. 


12. Summary. The role of the coalescence 
process in shower formation 


We have shown that on a number of days 
when cumulus clouds were carefully obser- 
ved, showers formed in clouds whose towers 
for the first time rose above a rather well- 
defined level, the fibrillation level. During 
the decay of towers not far surpassing this 
level, an evanescent fibrous structure, due to 
the presence of precipitation elements, could 
be seen. The height of the fibrillation level 
varied from occasion to occasion. 

Previous flight observations of the spectra 
of droplet size in cumulus cloud indicate the 
activity of the coalescence process of droplet 
aggregation, even in small clouds. It appears 
that with some reservations the course of the 
coalescence process towards shower formation 
can be represented by a simple model of 
cumulus structure, in which droplets are 
lifted in convective bubbles from near the 
cloud base to the summits, where they are 
subject to rapid evaporation unless they have 
attained the size of precipitation elements 
(minimum radius 150 w). The model is 
applicable to large clouds containing vigorous 
updraughts, and not to small, weak shower 
clouds, in which the history of the entire 
cloud appears more important. 

The speed of bubbles rising inside building 
clouds can be estimated on particular days as a 
function of height, by observing the speeds 
at which bubbles rise out of clouds to produce 
towers. Assuming also that inside such clouds 
processes are not far from adiabatic, the course 
of the coalescence process within an ascending 
bubble can be estimated, and can be shown 
to have progressed to the stage of producing 
precipitation elements by the time the bubble 
reaches the observed fibrillation level. 

It is concluded that the coalescence process 
has been responsible for the formation of 
the observed showers, even though the cloud 
tops sometimes have been strongly super- 
cooled. 

Other observations, to be examined in a 
second paper, show that the ice phase may 
develop in the tops of clouds when they are 
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supercooled, even only slightly, but apparently 
not usually before the coalescence process 
causes the formation of a shower. It seems 
that the coalescence, having produced large 
drops of which some freeze, introduces crystal 
multiplication processes which after a time 
lead to a glaciation, or more extensive devel- 
opment of the ice phase than can occur only 
by the action of freezing nuclei in small cloud 
droplets. Glaciation in large clouds frequently 
causes what CRADDOCK (1949) has aptly 
called a “second phase” in the growth of 
cumulonimbus, in which new towers rise 
beyond the fibrillation level with high speeds 
and produce impressive anvils. These upper 
appendages to cumulonimbus may, however, 
have little significance for shower formation. 
Although in general there is a close relation 
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between glaciation and shower formation, it 
appears to us that these are not necessarily 
cause and effect, but are often both conse- 
quences of an advanced coalescence process. 
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The Micro-structure of Cumuli in Maritime and Continental Air 


By P. SQUIRES, Division of Radiophysics, C.S.LR.O., Sydney, Australia 


(Manuscript received June 6, 1956) 


Since 1950 a number of measurements have 
been made of the droplet spectra in clouds, 
chiefly cumuli, using the sampler described by 
SQUIRES and GILLESPIE (1952). This device takes 
ten sampling rods at each loading, and these 
can be exposed at intervals of about three 
seconds, taking some thirty seconds in all. In 
sampling cumuli, the general practice has been 
to make horizontal runs at various levels, to 
obtain representative data for the whole cloud. 

It has been found that the characteristics of 
the droplet spectrum change rather erratically 
from one sample to the next. It is evident that 
the “grain” of a cumulus is smaller than the 
distance flown in three seconds, that is, about 
200 m. Nevertheless, certain types of clouds 
have distinctive characteristics, and the erratic 
variation within the one cloud is not so great as 
to mask these characteristic differences. 

One of the most striking differences is that 
between cumuli in maritime and in continental 
air masses. In maritime cumuli the droplets are 
bigger, and the droplet concentrations smaller, 
than in continental cumuli. At corresponding 
levels, there seems to be no notable difference 
in liquid water content between the two types 
of cloud. The difference in micro-structure is 
most simply discussed in terms of droplet 
concentrations. 

Only those cumuli which were forming in 
an unmodified maritime air mass and which 
were sampled over the sea have been included 
in the class of “maritime cumuli”. These have 
been measured off the eastern and southern 
coasts of Australia at various seasons of the 
year, and off the island of Hawaii during 
Project Shower, October to December 1954. 
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The class of “continental cumuli” includes 
only cumuli observed over the land in an air 
mass which was estimated to have spent at 
least two days over the land. These have been 
sampled in south-eastern Australia only. There 
is some evidence of a seasonal variation in the 
micro-structure of continental cumuli in this 
region, the lower concentrations occurring 
during the rainier autumn months. This vari- 
ation is not however so great as to confuse the 
rather marked contrast between maritime and 
continental cumuli, and it will not be discussed 
here. 

This contrast is illustrated in Fig. 1, which 
gives histograms of the percentages of all 
samples, irrespective of position in the cloud, 
for which the droplet concentration fell in the 
specified ranges. Table 1 shows for each air 
mass type the median, first and third quar- 
tile and maximum droplet concentration. 


Table ı 
; es Continen- 
Air mass type Maritime ai 
ıst quartile droplet concen- 

EPA ON ee en ee 19 150 
Median droplet concentration 42 290 
3rd quartile droplet concen- 

ONE m 66 470 
Maximum droplet concen- 

tration observed ........ 400 2 810 


It will be noted that the continental values 
entered in this table are all about seven times 
the maritime values. This implies that, if the 
droplet concentrations found in the continental 
clouds were divided by seven, the resulting 
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60 


u 
fe) 


MARITIME CUMULI 


O per cm? 


> 
° 


Total number of somples 481 


CONTINENTAL CUMULI 


Percentage of samples with concentrations in each ronge of 5 


500 
Droplet concentration per cm? 
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4 per cent of samples had concentrations over 200 per cm’ 


Total number ot samples 233 


over IOOO per cm? 


. Histograms of the percentages of all samples taken in maritime and continental cumuli in which the 


droplet concentration fell in the specified ranges. 


histogram would resemble that shown for the 
maritime cumuli in Fig. 1. Thus, if the conti- 
nental values are grouped in ranges of 350, 
instead of 50 as in Fig. 1, the first three ranges 
contain respectively 58, 26, and 5 per cent of all 
continental samples. The first three ranges of 
the maritime droplet concentration shown in 
Fig. I contain respectively 60, 28 and 5 per 
cent of all maritime samples. As regards this 
series of observations, one may conclude that 
the continental cumuli had, in general, droplet 
concentrations about seven times those found 
in maritime cumuli. It may be remarked that 
minimum values of the concentration have no 
significance, since samples taken at the edge of 


a cloud will sometimes show very low values. 
quite unrepresentative of the cloud as a whole, 

It is well known that maritime cumuli are 
more liable to yield rain without the presence 
of ice than continental cumuli. It seems possible 
that the difference in microstructure described 
is partly responsible for this, and it is hoped to 
enlarge on this and other aspects of the micro- 
structure of clouds in a later paper. 
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The Electrification of Individual Cloud Droplets 


By S. TWOMEY 


Division of Radiophysics, C. S. I. R. O., Sydney, Australia 


(Manuscript received May 16, 1956) 


Abstract 


The results of some observations of the electric charges carried by individual cloud elements 
are described. Measurable charges were found on many droplets, and a reversal of polarity was 
observed which appeared to be associated with the presence of ice particles. The different 
mechanisms of charging are discussed in the light of the experimental results, and reference is 
made to the possible significance of droplet electrification in the precipitation process. 


I. Introduction 


The physical properties of many colloidal 
suspensions are known to depend on the 
electric charges carried by the suspended 
particles. By analogy it would be expected that 
the charges carried by the droplets in a natural 
cloud could be an important parameter where 
cloud stability and coagulation are concerned. 
Comparatively little is known, however, about 
this aspect of cloud physics. A few measure- 
ments have been made (GUNN, 1952; WEBB 
and GUNN, 1955; WIGAND, 1926) but in all 
cases only the net charge per unit volume of 
cloud could be measured by the methods 
used; no information was therefore available 
concerning the charges carried by individual 
cloud elements and the relation (if any) 
between charge and size. 

In experiments in 1954 a technique was 
developed which allowed the charge on indi- 
vidual cloud elements to be determined in 
sign and magnitude. Some preliminary results 
were obtained by allowing the droplets under 
observation to fall through a liquid (kerosene). 
These results indicated that many positively 
charged droplets were present in water clouds, 
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while on some occasions negative charges 
were measured which appeared to be associated 
with the presence of the ice phase in the cloud. 
The method was then modified to allow the 
measurements to be made with droplets 
falling in still air. The present paper is mainly 
concerned with reporting results obtained 
with the modified technique; reference is also 
made to the possible implications of the re- 
sults obtained. 


2. Experimental technique 


In order to measure their electric charge, 
droplets freely falling between a pair of parallel 
lates are subjected to a strong horizontal 
electric field and the deflection recorded photo- 
graphically. The applied voltage (6-12 kV) 
is switched in a manner which results in a 
charged droplet tracing out an asymmetrical 
“zig-zag” path from which terminal velocity, 
charge-to-mass ratio, and the polarity of the 
charge are directly obtainable. (As terminal 
velocity is a known function of drop size, diam- 
eter and charge can then be readily calcu- 
lated.) Falling patterns of charged droplets 
in the applied field can be seen in Figs ı and 
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Fig. 1. Photograph of the traces produced by a cloud of 
highly charged small droplets in the laboratory. 


AE a mr J 


é 


Fig. 2. Photograph of a trace produced by a highly 
charged natural cloud droplet; d=+3x10-? e.s.u., 


q=20 4). 


2. The measurements herein reported were 
more sensitive than the earlier measurements, 
as a higher voltage was available. As used, the 
method was suitable for droplets in the dia- 
meter range 5-100 although it was best 
suited for observations of droplets in the 
range 10-40 u. The minimum charge measur- 
able was dependent on droplet size and was 
as follows: 


Droplet diameter Minimum charge 


(1) (8: satt.) 
(u) (e.s.u.) 
Io 107% 
so 10? 

100 10-8 
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By further increasing the applied field it 
would be theoretically possible to reduce the 
limits to one-eighth of the above values but 
this would involve the use of an applied field 
very close to the breakdown field in air, and 
this seemed most undesirable, both for instru- 
mental reasons and because of the possibility 
that the release of ions by corona discharges 
might give rise to spurious charges on the 
droplets under observation. 

In the laboratory it was possible to produce 
uncharged droplets which remained unde- 
flected in the highest field available. This indi- 
cated that the droplets did not acquire a charge 
during their passage into the observation 
chamber and through the applied electric field. 
There was in fact no obvious way in which 
the droplets sampled in natural cloud could 
become artificially charged-they were per- 
mitted to fall freely through a large opening 
in the roof of the trailer in which the apparatus 
was housed, and entered the observation cham- 
ber through an opening (5 cm x $ cm) in 
the top of the chamber. However, as an 
added precaution, when an observation was 
being made, the high tension voltage was 
applied and the camera shutter opened simul- 
taneously. 


3. Experimentel results 
(a) Observations in Winter 1955 

During July 1955 the apparatus was set up 
at the summit of Mt. Wellington, a 4,160 ft 
peak near Hobart, Tasmania. Measuremen s 
were made in stratocumulus and stratus clouds, 
winds being mainly in the south-westerly 
sector. With winds from this direction the 
sampling location was exposed to clean mari- 
time air which had travelled several thousands 
of miles over the Southern Ocean. The possi- 
bility of the results having been affected b 
man-made contamination could therefore be 
disregarded. 

The observations showed that a large pro- 
portion (around so %) of the cloud droplets 
sampled carried detectable charges. The samples 
were taken during the afternoon and early 
evening and included both precipitating and 
non-precipitating clouds. Cloud bases were 
on the average about 500 ft below the summit. 
The temperature at the sampling location 
varied from a few degrees below freezing to a 
few degrees above freezing. 
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Positive charges almost exclusively were 
found on the cloud droplets when the liquid 
phase only was present in the cloud. In such 
cases, 80 % or more of the charged droplets 
carried positive charges. In Fig. 3a charge 
has been plotted against diameter for all 
measurements of positively charged droplets. 
On occasions when negative charges were 
found it was noted that the ice phase was pres- 
ent in the cloud, ice crystals being visible in 
the intense beam of the illuminating lamp. 
The percentage of negative charges among 
the charged droplets varied from 40% to 
100 % when the ice phase was present. The re- 
sults of the observations of negative charges 
are plotted in Fig. 4. Throughout the obser- 
vations a proportion of the traces recorded 
represented droplets with appreciable charge, 
but were distorted (e.g. by turbulence within 
the observation tank) so that polarity could 
not be determined. All such observations 
have been omitted, as also have the many 
observations in which parts of the traces of 
large charged drops were recorded, but where 
enough of the “‘zig-zag” patterns to determine 
polarity was not available. 

These experiments confirmed qualitatively, 
at least, conclusions previously drawn from 
experiments carried out on Mt. Wellington 
in winter 1954, utilizing the deflection of 
droplets in kerosene. The conclusions referred 
to were as follows: (i) the majority of cloud 
droplets sampled at a mountain summit carried 
appreciable electric charges; (ii) on any one 
occasion one sign usually predominated, al- 
though mixed charges were sometimes ob- 
served; (iii) positive droplets were found in 
clouds which were thought to be composed of 
liquid water droplets only-negative charges 
were found only when the ice phase was 
present, but it could not be ascertained whether 
these negative charges were carried by water 
droplets or ice crystals or both. 


(b) Observations in Summer 1956 


In order to ascertain whether similar charges 
are found in warm clouds, a second series of 


1 A few isolated measurements of cloud droplet elec- 
trification made in trade-wind cloud on the slopes of 
Mt. Haleakala, Maui, Territory of Hawaii, in January 
1955, with the same equipment, led to signs and magni- 
tudes similar to those reported and discussed here. Drop 
diameters were from 74 to 264 and charges ranged 
from +10-7 e.s.u. to 4 X I0-® e.s.u. 
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Fig. 3 a. Relation of observed charge to droplet di- 
ameter (positive charges)—July 1955 results. 
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Fig. 3 b. Relation of observed charge to droplet diameter 
(positive charges)—February 1956 results. 


observations was carried out during February 
1956, at the same location. Measurements were 
again made in stratocumulus clouds, with 
maritime air from the east, cloud base some 
1,000 — 2,000 ft below the mountain summit 
and cloud tops rarely more than 2.000 ft above 
the summit. The temperature at the sampling 
location ranged from 8 °C to 17°C. It is cer- 
tain therefore that the clouds rarely, if ever, 


reached freezing level and it can safely be 


assumed. that the ice phase was absent. 
The results, which are plotted in Fig. 3b, 
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showed that about 80% of the cloud droplets 
carried positive charges, negatively charged 
droplets being absent. The relation between 
charge and diameter was very similar to that 
found for the positive charges measured in 
winter time. 

The above results thus confirmed the pre- 
viously deduced association between liquid 
cloud droplets and positive charge. The clouds 
sampled were typical (usually non-precipitat- 
ing) stratocumuli; the diameters of the cloud 
droplets lay mainly in the range 10-40 u, 
with a median drop diameter between 25 w 
and 30 u. 


4. Discussition 
(a) Positive charges 


From Fig. 1a, in which the observed values 
for droplet charge have been plotted against 
droplet diameter for the positively charged 
droplets sampled in July 1955, it is seen that 
there appeared to be a relationship between 
charge and diameter of the form q = Ad", 
the exponent n lying between 2 and 3. The in- 
dividual values of course showed large devia- 
tions, but there certainly appeared to be a 
significant relationship: the correlation cocfli- 
cient drawn from the original data for Fig. 3 a 
was 0.85, while the regression equation was 


q = 5.4 x 10710 2-34 (where q was in e.s.u. 
and din microns). 


For comparison the charge which would give 


the breakdown field strength (for air) at the 


20F 


Diameter (u) 


L À. 


si. 1 
10° 107? 10° 10° 104 10? 
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Fig. 4. Relation of observed charge to droplet diameter 
(negative charges)—July 1955 results. 
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surface of the droplet has also been plotted 
in Fig. 3a-this probably represents the 
maximum charge which a droplet of given 
diameter can carry. It is seen that the charges 
measured were always considerably less than 
a tenth of this limiting value. 

The data plotted in Fig. 3b yielded a corre- 
lation coefficient of 0.80, and the following 
regression equation: 

g =9.9 x 10710 d2- 08, 


The latter is seen to be very similar to the re- 
gression equation drawn from the earlier data. 

In Fig. 4 charge has been plotted against 
droplet diameter for the negatively charged 
droplets sampled. It is seen that the scatter 
was greater than was the case with positive 
charges—in fact there is very little ground for 
drawing a curve through the points plotted. 
The correlation coefficient for these points 
was 0.48, considerably less than that obtained 
in the case of positively charged droplets. 

The relationship existing between charge 
and diameter for positively charged droplets 
suggested that a straightforward charging 
mechanism was responsible. Two processes 
which would give rise to such an ordered 
distribution of droplet charge are (a) diffusion 
of ions into a droplet (GUNN, 1955a) and (b) 
the electrification observed by BLANCHARD 
(1955) when minute droplets are formed by 
bursting bubbles at an air-sea water interface. 
However, neither of these mechanisms provides 
a complete explanation of the observed effects: 
Gunn’s theory and also his laboratory experi- 
ments indicate that the charges would be 
distributed statistically with a slight excess of 
positive charges, due to the superimposition 
of a very small systematic positive charge. 
The magnitude of the charges given by Gunn’s 
theory is also considerably less than the charges 
measured in these experiments! — for example, 


! However it has been suggested by Telford (private 
communication) that a field much higher than the 
breakdown field can exist at the droplet surface without 
the occurrence of corona discharge; this arises from the 
fact that the field falls away rapidly with distance from 
the surface and hence the average field over a few diam- 
eters may be less than the breakdown field although 
the field at the surface exceeds this value. 

!In a recent paper, GUNN (1955b) shows that mean 
charges of the same order of magnitude as those discussed 
here, can arise as a result of the collisions of droplets 
charged by ionic diffusion; equal numbers of positive 
and negative charges, however, are expected. 
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the mean charge for droplets with diameter 
20 u would be only 16e (7-7 x 107! e.s.u.) with 
a fraction less than 0.005 of the total population 
possessing charges in excess of twice this mean 
charge. The corresponding observed values 
were of order 10% e.s.u. 

Blanchard’s results were obtained for drop- 
lets of sea water, newly formed from bursting 
bubbles at an air-sea-water interface. Although 
the charges found by Blanchard were positive 
and of the same order of magnitude as those 
reported here, the suggestion that this process 
could account for the charges found on natural 
cloud droplets is subject to criticism on two 
main grounds. Firstly, a large proportion of 
the total population of cloud droplets observed 
carried measurable charges; the concentration 
of saline nuclei in the free air is almost cer- 
tainly less by several powers of ten than the 
concentration of droplets in a cloud, so that 
it is not at all likely that each charged droplet 
observed contained a sea salt nucleus. Secondly, 
the work of GUNN (1954, 1955a) suggests that a 
highly charged droplet would be rapidly 
discharged (in a time of the order of ten min- 
utes) to an equilibrium value determined by 
ionic diffusion. For these reasons therefore it 
seems reasonable to disregard the bursting 
bubble mechanism as the main source of 
cloud electrification. 

At the present juncture a diffusion process of 
the kind described by Gunn seems most 
applicable to the explanation of the observed 
effects. The order of magnitude of the charges 
computed by Gunn would no doubt be greater 
if the adsorption of ions was selective. It is 
known, mainly from balloelectric phenomena 
(NATANSON, 1950; PANNETIER, 1952), that a 
regular arrangement of ions exists at a water 
surface, negative ions congregating at the 
surface in the case of pure water, with a 
balancing excess of protons in the interior of 
the liquid. This dipole layer at the surface 
does not affect the electric field external to a 
droplet but it is possible that forces may exist 
in the transition layer which allow positive ions 
to be adsorbed more readily. PHILLIPS and 
Gunn (1954), however, while measuring the 
electrification of a metal sphere by moving 
ionized air, observed that surface contamina- 
tions did not have any appreciable effect on the 
charge collected by the sphere. That selective 
adsorption would be sufficient to account 
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for the charges measured follows from Gunn’s 
computations and also from the work of 
ArENDT and KALLMAN (1925). The latter 
measured the charge on droplets in an environ- 
ment in which ions of one sign only were 
present and found charges which would lead to 
charge-to-mass ratios and surface fields greater 
than those reported here. 

The contention that the observed charges 
were produced by a diffusion process is sup- 
ported by the relationship between charge 
and droplet diameter. It has been shown 
that the charges measured tended to vary 
according to the 2.34 power of the diameter, 
in the case of the observations in July 1955, 
or the 2.05 power in the case of the observations 
in February 1956; computations of the charge 
produced by diffusion always lead to an 


approximately square-law relationship. 
(b) Negative charges 


The negative charges found seemed to have 
a connection with freezing or melting of the 
cloud elements, charges of this sign being 
apparently associated with the presence of the 
ice phase in the cloud. As mentioned earlier 
it was not possible to ascertain whether the 
negative charges were carried by liquid 
droplets or ice particles or both. However, it is 
fairly certain that some at least of the negative 
charges were carried by ice crystals because 
it was sometimes possible to identify the traces 
of falling ice crystals by the discontinuous, 
“dashed” appearance of the traces (corre- 
sponding to the scintillations often ‘observed 
visually from ice crystals). Some of the nega- 
tively charged elements gave traces of this 
kind and it was therefore concluded that at 
least some of these elements were in fact ice 
particles. However this did not serve to indi- 
cate whether ice particles alone, or both ice 
particles and liquid droplets, carried negative 
charges. However, it is relevant to note that 
LUEDER (1951) found that a body collecting 
supercooled droplets by being rotated through 
a natural cloud below o °C acquired in the 
process a large negative charge; the maximum 
charge density was given as 0.2 coulombs 
per cm? of collected water and it was found 
that the charging process ceased when the 
potential gradient at the surface reached a 
limiting value of 700 volts/cm. The charge 
which would produce this potential gradient 
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at the surface of a droplet or spherical ice 
particle has been plotted on Fig. 4. It may 
of course have been fortuitous that this line 
fell where it did, but it is definitely thought- 
provoking that the points plotted did appear 
to be bounded by this line. If in fact we 
suppose that the observed negative charges 
were produced by a mechanism of the type 
found by Lueder, it is seen that the sign and 
magnitude of the charges are readily explained. 
The scatter of the charge values would also 
be readily explainable if the charging was 
caused through the collection by a frozen 
droplet or ice crystal of supercooled liquid 
droplets by collision: according to Lueder’s 
findings the charge produced on a cloud 
particle would be a linear function of the vo- 
lume of supercooled water collected, which 
would be expected to vary considerably from 
particle to particle. It is of some interest to 
compute, for a few values of droplet diameter, 
the volume of supercooled water the collection 
of which would give rise to the maximum 
observed charge values on the assumption 
that each cm? of collected supercooled water 
causes a charge of 0-2 coulombs (6 x 10? e.s.u.). 
Reference to Table 1 shows that on this hypo- 
thesis the freezing of a droplet, the volume of 
which represents a very small fraction of the 
total volume, would be sufficient to account 
for the largest charges observed on cloud 
elements. It is possible to account for the 
largest individual charge observed (2.5 x 1078 
e.s.u.) by the collection of one supercooled 
droplet with a diameter of only 2. The 
cessation of charging at a potential gradient 
of 700 volts/cm as found by Lueder would 
explain the absence of extremely high charges. 

From the foregoing it is seen that the mech- 
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anism of charge generation by collection 
of supercooled water furnishes an attractive 
hypothesis for the interpretation of the obser- 
vations concerning negatively charged cloud 


droplets. 


(c) The possible role of electric charges in cloud 
physics 

The electrification of cloud droplets may 
affect cloud stability in two ways: the coales- 
cence of cloud droplets may be inhibited if 
charges of the same sign and of sufficient 
magnitude are carried by the droplets; simi- 
larly the coalescence rate may be increased if 
the droplets carry charges of adequate magni- 
tude and mixed polarity. Appreciable effects 
produced by electrification have been observed 
in laboratory experiments on the coalescence 
of droplets by TELFORD, THORNDIKE and Bow- 
EN (1955). Both the sizes of the drops and 
the charges involved, however, were consid- 
erably larger than those discussed here. 

The possibility that electric charges on 
cloud droplets may be sufficient to enhance or 
diminish the stability of a cloud is particularly 
interesting when considered in conjunction 
with recent theoretical work by TELFORD (1955) 
who examined some statistical aspects of the 
mathematics of droplet collisions. Previous 
computations (BOWEN, 1950; LANGMUIR, 1948) 
indicated that the formation by coalescence of 
drops of raindrop size would require times 
of the order of one halfhour or more. However, 
Telford’s computations suggested that rain- 
drops can be formed within a few minutes, 
assuming that the collection efficiencies - of 
cloud drops are determined by the aerodynam- 
ics of the process alone. This finding, howev- 
er, is contrary to observations, which show 
that large clouds often persist for times of 
the order of an hour or more without precipi- 
tating. These apparently conflicting factors can 
be reconciled if one assumes that collection 
efficiency is reduced greatly because the drop- 
lets in a cloud carry charges of the same 
polarity and of sufficient magnitude to inhibit 
the coalescence process. 

The charges measured were probably large 
enough to affect coalescence, particularly that 
of the smaller droplets (for which the charge-to- 
mass ratio was greatest). This was suggested by 
consideration of the magnitude of the elec- 
trostatic forces between charged droplets in 
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close proximity; for example, the force 
between a pair of 1o'w droplets carrying equal 
charges and just in contact would equal the 
weight of either droplet if the charge carried 
was 7X 10-7 e.s.u. The actual values of the 
observed charges were not much below this 
figure. TELFORD et al. (1955) found that charges 
of 3x1074 e.s.u. were sufficient to increase 
the coalescence rate in the case of 150 drops; 
on the assumption that the charge-to-mass 
ratio is the appropriate one to consider where 
coalescence problems are concerned, it follows 
that the observed charges would indeed have a 
significant effect on the collection efficiencies 
of cloud elements. 

The above considerations suggest that there 
may be an aspect of the Bergeron precipitation 


. . to) 
mechanism which has not been given much 


attention: if a cloud composed entirely of liquid 
water droplets possesses a degree of stability 
due to the positive charging of a considerable 
proportion of the droplets, then the glaciation 
of part of the cloud may drastically affect the 
cloud stability, not only as a result of the va- 
pour-pressure difference between water and 
ice, but also because the growing ice particles 
may acquire a negative charge sufficient to 
increase their collection efficiency when they 
fall to the lower and warmer parts of the cloud 
where positively charged droplets are present. 

It must be pointed out that the above sug- 
gestions are largely conjectural. They are based 
on measurements made in cloud on mountain 
tops, which do not necessarily provide results 
which can be applied generally to clouds in 
the free air. However, measurements of this 
kind are the only source of information avai- 
lable at the present juncture. If in fact the 
charges measured were a result of the 
mechanisms suggested earlier, it seems very 
likely that the charges in a free cloud would 
follow a pattern similar to that described. 


5. Conclusions 


Measurements made on a mountain summit 


showed that a high proportion of cloud ele- 
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ments sampled carried appreciable electric 
charges. Drops in liquid water clouds were 
consistently positively charged, negative or 
mixed charges occurring in clouds containing 
the ice phase. The magnitudes of the positive 
charges were clearly correlated with diameter, 
ranging from around 10°? e.s.u. for droplets 
10 microns in diameter to above 105 e.s.u. 
for 100 micron droplets. The negative charges 
were of a similar order of magnitude but were 
not clearly correlated with droplet diameter. 

The charges measured should be sufficient to 
reduce significantly the coalescence rate in the 
case of cloud elements of the same polarity or 
to increase it in the case of encounters between 
elements of opposite polarity. Hence it seems 
probable that electric charges play an impor- 
tant part in the physics of rain formation. 

It is suggested that positive charges were 
accumulated by diffusion of ions into the 
droplets. It seems necessary to postulate that 
positive ions are preferentially adsorbed, 
possibly as a result of the arrangement of 
molecules in the liquid surface layer. The 
production of negatively charged cloud ele- 
ments is attributed to a process of the kind 
described by Lueder, in which the collection 
of supercooled water resulted in the collecting 
body acquiring a negative charge. 

Further observations, both in freezing and 
non-freezing cloud, are desirable. Valuable 
information would result if a technique could 
be developed for airborne measurements of 
the charge carried by individual droplets. 
Laboratory investigations into the charging of 
droplets formed by condensation and the ice 
crystals formed by the freezing of such droplets 
might also make a useful contribution to our 
knowledge of cloud electrification and cloud 
physics generally. 
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Abstract 


The results of seeding thirty-five supercooled clouds, mostly cumuliform, are described. 
The silver iodide was released as a smoke from aircraft flying through the upper levels of the 
cloud or just below cloud base. It is concluded that silver iodide is an effective agent in inducing 
precipitation in clouds whose upper levels are at temperatures colder than —5° C, and that 
the precipitation is usually released 20 to 25 minutes after seeding in the case of cumulus clouds 


but later in the case of layer clouds. 


1. Introduction 


Although the effectiveness of silver iodide 
as a nucleus for ice crystal formation, and 
consequently its value for cloud seeding, was 
first reported by Vonnegut in 1947 there is 
still quite inadequate documentation of its 
effectiveness in releasing precipitation from 
individual clouds. Indeed most of the literature 
on the use of silver iodide for cloud seeding 
relates to experiments in which the material 
has been released as a smoke from the ground 
in an attempt to influence the rainfall over a 
comparatively large target area downwind 
from the ground generator. Decisions about 
the effectiveness of the seeding have been 
based largely on a statistical analysis of the 
rainfall distribution, both in time and space, 
and considerable doubt exists amongst meteo- 
rologists about the validity of many of the 
conclusions reached. Comparatively few air- 
craft experiments have been reported, and 
most of these have also been aimed primarily 
at seeding a target area; little attention has 
been paid apparently to observations of the 


Tellus VIII (1956), 4 


growth and development of the seeded clouds, 
results being assessed as before by a statistical 
analysis of rainfall. 

This paper presents the results of some thirty- 
five aircraft experiments in seeding individual 
clouds in which the cloud was kept under 
observation for at least 30 minutes, and often 
for up to 60 minutes after being seeded, and its 
development compared with that of its un- 
seeded neighbours. The data presented de- 
monstrate fairly conclusively the effectiveness 
of silver iodide as a seeding agent in the field, 
and confirm much laboratory experience about 
its action. For the most part the experiments 
were made in relatively isolated cumulus 
clouds located in southeastern New South 
Wales about 100 miles inland from the coast; 
a lesser number of experiments was made in 
inland South Australia and North Queens- 


1 At a conference on weather modification studies 
held at the University of Arizona in April 1956 aufm. 
Kampe, Kelly and Weickmann (private communication) 
reported the results of seeding supercooled stratus clonds. 
These experiments were aimed, however, at clearing 
thin cloud layers rather than releasing precipitation. 


454 


land, and only very few experiments have 
been made in maritime conditions near the 
coast. 


2. Method of feeding 
(a) Flight procedure 


Three different aircraft were used in these 
seeding experiments: an Auster, an Anson 
and a DC3. The first two aircraft were not 
fitted for instrument flight in cloud and 
seedings with them were made just below 
cloud base. Thereafter the aircraft stayed in 
the vicinity for half to one hour while the 
cloud was kept under observation. The DC3 
aircraft was fully instrumented and nearly 
all seedings made with it were at a level, 
usually near the cloud tops, at the temperature 
where silver iodide was considered to be 
immediately active as an ice nucleus. Since 
the aircraft was also fitted with radar the 
presence of precipitation in the cloud could be 
observed well before it had reached the base. 
Observations of the results of seeding were 
made both at the seeding level inside the cloud 
and from below cloud base. 

For the most part seedings were made only 
when favourable experimental conditions were 
known to exist: that is, when cloud top tem- 
perature was known or expected to be — 4° C 
or colder and the cloud was not already pre- 
cipitating, surrounding clouds were not al- 
ready raining, and the chosen cloud appeared 
to be durable. Wherever possible a control 
cloud or clouds in the vicinity of the seeded 
one, and of equal or greater initial develop- 
ment, was also kept under observation for the 
full period of the experiment. 


(b) Generation of silver iodide nuclei 


In the case of the Auster aircraft, nuclei 
were generated by injecting an acetone solution 
of silver iodide into the engine exhaust. 
A small hand-operated pump injected roughly 
5 cc of solution at a time, and usually a total 
of 50 to 100 cc, containing about 10 gm of 
Agl was used for seeding each cloud. Special 
burners were developed for use with the 
Anson and DC3 aircraft. These consisted of a 
long flame tube in which the acetone solution 
of Agl was injected as a fine spray where it 
was ignited by a spark from an engine type 
spark plug. Ram pressure delivered sufficient 
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air to produce a hot flame which was quenched 
by excess air just before it left the rear of the 
flame tube. Burners were fitted at the wing 
tips of the Anson and a few feet below the 
wing of the DC3 and were supplied with 
the silver iodide solution under pressure from 
tanks in the main cabin. The burners were 
designed primarily for continuous seeding of a 
large area and consumed 1 to 2 gallons of 
solution per hour of operation, each gallon 
containing soo gm of Agl. The burner was 
kept on for the complete passage under or 
through the cloud, and at the aircraft speeds 
used, about 4 gm of Agl would be distrib- 
uted per mile. Thus approximately 10 gm 
of Agl was injected into each cloud. These 
burners have been tested for effectiveness 
on the ground, and of the order of 1018 nuclei, 
active at —10° C, are produced for each gram 
of Agl burnt. The threshold temperature at 
which any of the smoke particles become 
active as ice nuclei is about -4 to -6°C. 
It is known from airborne cold box tests 
that active nuclei are produced in flight, and 
that the burner efficiency under these condi- 
tions is similar to that on the ground. 


3. Results 


The results of all the experiments are sum- 
marized in Table 1, which lists the most 
significant data and observations obtained. 

The cloud base temperature was measured 
on a calibrated aircraft thermometer; the 
temperature at the cloud top was not always 
measured however, particularly when the 
seeding was done at the base level. The figure 
quoted is then estimated from available radio- 
sonde and synoptic information, taking account 
of the position of any fronts and the degree 
of coincidence between the radiosonde and 
measured values of temperature at cloud base 
or at the highest level reached by the aircraft. 
This figure is thus sometimes only an estimate, 
but it is considered unlikely to be in error by 
more than 2° C, except in the case of ver 
deep clouds seeded at base level, where the 
error may be 5° C. 

The time interval given in column 9 is 
that between the commencement of seeding 
and the first appearance of definite precipita- 
tion under the cloud base. Intensity of precip- 
itation was estimated both visually by flying 
the aircraft through the rain and where pos- 
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sible by radar, and the descriptions given 
may be taken roughly to imply the following: 
Heavy rain 15 mm/hr Moderate rain 5 mm/hr 
Light rain 1.5 mm/hr Drizzle 0.2 mm/hr. 

It is clear from the summarized information 
in the table that the degree of confidence of the 
cause and effect relationship between seeding 
and precipitation varies from case to case. In 
the interests of brevity it is proposed to discuss 
in detail only a few of the cases given and only 
those in which fairly clear evidence can be 
produced for the effectiveness of seeding. 

Case A, 15 October 1953. Nhill, Victoria. 

General synoptic situation : High pressure sys- 
tem centred 100 miles west of Tasmania, 
giving rise to light easterly stream over most of 
south-east Australia. 

Summary of flight: The aircraft took off 
from Nhill at 1545 hours, when three isolated 
cumuli were observed to the south in an 
otherwise clear sky. These clouds arose as a 
result of three scrub fires which were burning 
in the area—an isothermal layer at 3,200 ft 
apparently prevented the formation of con- 
vective clouds elsewhere. Cloud base was at 
7,000 ft where the temperature was + 4.5° C, 
and cloud tops at 13,000 ft. The three clouds 
were almost identical in both size and shape, 
with base dimensions approximately equal 
to the cloud depth (6,ooo ft). The clouds were 
dense and compact. One cloud was selected 
and seeded underneath its base at 1628 hours, 
prior to which no changes had occurred in 
any of the three clouds. The seeded cloud 
was then kept under observation and periodic 
traverses made underneath its base. At 1650 
hours a few droplets appeared on the wind- 
screen during such a traverse; at 1655 hours 
rain was falling to the ground. The precipi- 
tation was light to moderate and occurred in a 
column which was traversed in I minute 
(i. e. about 1! miles wide), with a small 
region of moderate to heavy rain which 
occupied about one quarter of the entire width 
of the rain column. The rain ceased about 
1705 hours, by which time the seeded cloud 
had largely dissipated; shortly afterwards 
only a small and ragged remnant of the original 
cumulus congestus remained. The control 
clouds showed no visible change in shape or 
appearance throughout the period of obser- 
vation. 

Case B, 4 May 1955. Uranquinty, N.S.W. 
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General weather: Major low pressure sys- 
tem centred soo miles east of N.S.W. coast; 
secondary low pressure centred about 100 
miles NNE of operational area, with associated 
trough line extending to the NW; southerly 
maritime polar air over most of south-east 
Australia. 

Summary of flight: In the vicinity of Wagga 
an extensive field of cumuli was observed 
during the morning. An individual cumulus, 
together with two other slightly larger cumuli 
for use as control clouds, was selected at 1145 
hours. Cloud base was at 4,000 ft where the 
temperature was at+6° C and the tops were 
estimated to be at 12,000 ft. 

The aircraft entered the selected cloud at 
1154 hours at 10,500 ft, the clear air tempera- 
ture being —10° C. A well marked updraft 
was encountered for 45 seconds near the 
middle of the cloud. Silver iodide continued 
to be released until the aircraft emerged from 
the cloud at 1156 hours. A steep 90° turn was 
then executed and the cloud perimeter circled 
as the aircraft descended slowly. The cloud 
was approximately circular in cross-section 
with a diameter of 3 to 5 miles. At 1206 
hours a faint rain echo appeared on the radar 
screen at a range of 1,500 yards. This echo 
grew rapidly and was seen to extend to the 
ground at 1211 hours. At 1220 hours the 
aircraft was below cloud base and a traverse 
was made through the rain belt. Precipitation 
was heavy and concentrated in a column 
which was traversed in 40 seconds. The seeded 
cloud was still raining when last observed at 
1245 hours. During the period of observation 
no change in appearance or precipitation was 
observed in the two control clouds. The 
seeded cloud only changed appearance well after 
precipitation had reached the ground when 
the typical anvil head of a cumulo-nimbus 
cloud appeared. 

Case C. 26 July 1955, Hobart, Tasmania 

General weather: Low pressure system 
centred 500 miles to west producing NW 
winds. A few thin Cu over land and extensive 
isolated layers of Ac over land and sea. No 
high cloud apart from a trace of Cs to NW. 

Summary of flight: A temperature sounding 
showed that freezing level was at 6,500 ft 
and the Ac was based at 10,000 ft, where the 
temperature was -8°C, and had its tops 
between 12,000 and 12,500 ft where the 
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temperature was -12.5°C. A sheet of Ac 
was chosen for seeding which was roughly 
80 miles by 50 miles in extent and fairly uni- 
form in appearance except that it was somewhat 
deeper at its southern limit where the base 
descended to 9,000 ft. A seeding run was made 
running from north to south entering cloud 
at 1145 hours at 12,000 ft but descending to 
11,500 ft at 1148 hours to keep wholly within 
the cloud. At 11.57 hours the aircraft com- 
menced to descend through the cloud while 
turning on a reciprocal course and breaking 
through the base at 9,300 ft at 1205 hours. 
The aircraft descended to 8,000 ft while 
heading north. Seeding was terminated at 
1200 hours. 

During the seeding run no sign of snow or 
ice crystals was observed in the cloud which 
appeared to consist wholly of supercooled 
water droplets. From the rate in which ice 
built up on small exposed surfaces and from 
general experience it was estimated that the 
cloud water content was between o.ı and 
0.2 gm/m?. The aircraft headed north till 
1214 hours when it turned back southwards 
again. No precipitation was visible anywhere 
till 1217 hours when the first streaks were 
seen at the southern limit of the seeded track. 
The aircraft entered this at 1225 hours and 
it was found to consist of small snow flakes 
2 to 3 mm across, the precipitation rate being 
estimated at between 1 and 1 mm/hour. 
The aircraft turned back northwards at 1229 
hours and precipitation was observed from 
the northern limit of the seeded track at 1235 
hours. Thereafter the aircraft kept flying up 
and down beneath the seeded track. By 1300 
hours the precipitation area had extended 
over the whole length of the seeded track 
and had broadened in the east-west direction 
to cover roughly 3 to $ miles. However, the 
precipitation did not reach the surface, mostly 
evaporating before it had reached 4,000 ft. 
The aircraft then returned to base, landing 
at 1340 hours. The cloud was still precipitating 
and showed no sign of evaporation along 
the seeded track. There was no precipitation 
anywhere in the cloud sheet except along 
this track and no other Ac in the area pro- 
duced precipitation. It is estimated that the to- 
tal amount of precipitation released at cloud 
base in the 2-hour observation period was 
I mm over 100 square miles. 
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4. Discussion of results 


(a) Effect of cloud top temperature on results of 
seeding 


In Figure 1 the results are grouped according 
to cloud top temperature. It can be seen that 
of the 29 clouds colder than — 5° C at their 
tops 72 per cent precipitated after seeding, 
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Fig. 1. The effect of cloud top temperature on the results 
of seeding. 


21 per cent evaporated and 7 per cent remained 
unchanged. 

Of the six clouds warmer than - 5° C only 
one cloud precipitated and one evaporated; 
the other four, or 67 per cent, showed no 
change. 


(b) Time to release precipitation 


The results are grouped in s-minute time 
intervals in Figure 2. Most clouds are seen to 
precipitate between 20 and 25 minutes after 
seeding. There is slight evidence that layer 
clouds tend to take somewhat longer to re- 
lease precipitation after seeding than do 
cumuliform clouds. The cases where the time 
interval was markedly shorter than 20 minutes 
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Fig. 2. The time interval between seeding and the re- 
lease of precipitation. 
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may be real but it is suspected that in some at 
least of these cases the rain was due to natural 
causes. There is also some evidence that when 
clouds are seeded at the base, the time taken 
to release precipitation is longer than when 
they are seeded at higher levels where the Agl 
smoke particles can act as ice erystal nuclei 
immediately, but as shown in Figure 3, there 
is little sign that the depth of the cloud has 
much effect on the time interval between 
seeding and the release of rain. When layer 
clouds are excluded it would seem that deeper 
clouds seeded at their bases take longer to 
rain than do shallower clouds seeded in the 
same way, but the correlation is not very 
great. Presumably the more violent updrafts 
in the deeper clouds compensate for their 


greater depth. 
(c) Quantity of silver iodide required for seeding 


No experiments have been made with 
very small quantities of silver iodide. Indeed 
in many cases the number of effective nuclei 
that have been distributed must have been 
many tens or hundreds of times the number 
of raindrops that eventually precipitated. This 
certainly occurred in case C reported above. 


20000 > 
® 
€ Cu seeded at base 
15000 - O Cu seeded near top 
À Layer clouds 
J | °° 
ee e 
£ 10000- ose 
a 
Q O 
o fe} 
3 & 
© O e 
© 5000- 
e A 
= a 
O- T T 1 
ie} 20 40 60 


Time to release precipitation (mins) 


Fig. 3. The effect of cloud depth on the time to release 
precipitation. 
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Even in this case, where the initial concentra- 
tion of smoke particles must have been very 
high, there was however no evidence of 
overseeding, or of the clearing of a lane in 
the cloud. It would seem from this and other 
similar experiments that such clearing is 
difficult to achieve with silver iodide seeding, 
at least with the moderate quantities used in 
these experiments. 


5, Conclusions 


From the experiments reported it seems 
clear that silver iodide is an effective agent for 
inducing precipitation in supercooled clouds 
whose temperatures are less than -5°C. 
The time taken to release precipitation is 
usually some 20 to 25 minutes after seeding 
in the case of a cumulus cloud, but may be 
longer in the case of layer clouds. The question 
as to whether the precipitation will reach 
the ground, clearly depends on the quantity 
released, the height of cloud base above the 
terrain and the dryness of the air between. 
In most cases however any reasonably deep 
cloud that was seeded produced rain that 
reached the ground. It is of interest to note 
that these conclusions are very similar to 
those of Squires and SMITH (1949) in relation 
to the effect of seeding supercooled clouds 
with dry ice. 
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By ROBERT R. LONG, The Johns Hopkins University, Baltimore 
(Manuscript received March 2, 1956). 


Abstract 


A theoretical discussion is given of the solitary wave. Part I is concerned with the wave on a 
water surface; Part II considers the solitary wave at the interface of two superimposed liquids of 
different density, bounded above and below by rigid surfaces. 

In Part I higher approximations are obtained for the speed of propagation, and limits to the 
wave speed are derived from the momentum theorem. The results include a correction to the 


second approximation of Weinstein. 


In Part II the first two approximations to the wave speed are derived. The equations of the 
second approximation yield limits on the wave amplitude. 


PART I 


THE SOLITARY WATERS WAVE 


I. Introduction 


The solitary wave of elevation on a water 
surface, shown schematically in Fig. 1, was 
first studied systematically by RusseLL (1844). 
His empirical formula for the speed of propaga- 
tion, c, is 


=g(h + a), (1) 


were g is the acceleration of gravity, h is the 
depth of the water at a great distance from the 
crest, and a is the height of the crest above the 
level h. This may be written 


F2=1 +4, (2) 


where F? is the Froude number, c?/gh, and « 
is a/h. 

Equation (2) was obtained theoretically by 
BOUSSINESQ (1871) and RAYLEIGH (1876) as 
a first approximation. Later McCowan (1891) 


1 This research was sponsored by the Office of Naval 
Research under contract N-onr-248 (31). 


in an attempt to improve the approximation, 
derived the following expressions defining the 
speed: 

tan 6 


B * 
x = sn 6 tan = B(x + a). (3) 


F2 - 


It is easily verified that this agrees with the 
Boussinesq-Rayleigh approximation to the 
order «. In the same paper, however, McCowan 
chose as a “much closer approximation” 


’ 


tan P 
Fie, 
B 
aes 3 I 
a= Sy sine 8 (1 +32) tan f(x +a) (4) 


To the order &?, this yields F? =1 +a - 23/6002. 

It is clear from McCowan’s method that his 
coefficient of «? is too large in absolute value. 
However, WEINSTEIN (1926) obtained F?= 
=1+x& —21/20x°, using methods rigorous to 
the order «?. In addition to the conflict with 
McCowan’s result, the Weinstein formula does 


Tellus VIII (1956), 4 


SOLITARY WAVES IN THE ONE- AND TWO-FLUID SYSTEMS 


not agree with the results of the present paper. 
These two considerations led the author to the 
discovery of a numerical error in the last step 
of Weinstein’s development. His corrected 
result is 


P=1+0-—a2+0(e). (5) 


= 


Numerous other studies have been made of 
the solitary wave. McCowan (1894) investi- 


gated the highest wave possible. Assuming a 
solution similar to that in his first paper, he 
concluded that the maximum wave would 
occur for 6 approximately equal to 1 in F?= 
=tan ß/ß. This yields F?=1.56 and, from the 
energy equation, «=0.78. 

More recently PACKHAM (1952), using a 
technique developed by Davies (1951), obtain- 
ed an approximate formula for the wave speed 
similar in form to that of McCowan. Finally 
KELLER (1954) and FRIEDRICHS and Hyers 
(1954) developed systematic approaches for 
first and higher approximations to the solu- 
tion. In the latter paper, the existence of the 
wave for sufficiently small « was established 
using methods differing from those of Lavren- 
tieff! who is credited with the earliest proof 
of existence. 

The purpose of the present paper is to derive 
higher approximations to the solitary wave. 
By direct integration F? is obtained as a power 
series in « as far as the «* term. An application 
of the momentum theorem yields one addi- 
tional term. The total result is 


F=1 u Bu ee 
20 70 27 
AT i +0 («°). (6) 


77 000 


The momentum theorem is further employed 
to obtain upper and lower limits on F? for 
values of « up to that of the maximum wave. 


1 See FRIEDRICHS and HYERS (1954). 
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2. Higher approximations to the solitary 
wave 


We assume the existence of a wave of 
permanent form with a single crest (Bean): 
The height of the crest above the horizontal 
bottom y=o is h+a. The free surface lowers 
symmetrically on both sides of the crest to 
the height y=h at |x|= co. The fluid motion 
is two-dimensional, irrotational and friction- 
less. The density is constant in time and uniform 
in space. 

If we translate the axes of our coordinate 
system with the wave, the flow becomes steady, 
the wave motion superimposed on a basic 
current, c. The dynamic equation at the free 
surface is 


q:2° =<? —2¢(y,—'h), (7) 


where q is the fluid speed and “s” denotes 
quantities on the free surface. In view of the 
irrotationality and incompressibility, a velocity 
potential p(x, y) and a streamfunction y (x, y) 
exist such that y+ iy is an analytic function of 
the complex variable x+iy and vice versa. 


Therefore 


[2 - 29 (ys 


the suffix “s” denoting derivatives at the free 
surface. Defining nondimensional quantities, 


DEI after Yas 
Si nr NN a Open er 


(9) 


we may write (8) as 


[n=2p(ye—1)| Ie. ~)s 4 CLO} 


We now assume solutions of the form used 
originally by RAYLEIGH (1876): 


FOR) y'2 Fay’ 
ol = “ Aid) 
4! 
2 ! FO) y'8 F6) y : 
y = FU y =: zt 5! — (12) 


F is a function of g’ alone, and F is its n' 
derivative. Choosing y’=o at the bottom of 
the channel (y’=o), we have y'= -1 at the 
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free surface and we satisfy the kinematic condition that De=o at the crest where e=«& 


conditions at both surfaces. Defining F = -ı 
—n(q), at the free surface, 


Dx; = -—1-cosD:n= -1-n+ 


LM (14) 


where D denotes the operator of differentia- 
tion with respect to @’, and Qn is the indefinite 
integral of 7. Setting e=y;-ı we find from 
equation (10) 


(1 — 2pe) [(1 + cos D- 9)? + (De)?] = 1. (15) 


Using (14), 
cos D-y=D cot D-ée= 
(2) (4) 
ir Ne ee 
5° 0245227945 
Baer 22ng(2n) 
VIN Doit (16) 


where B,,-, are the Bernoulli numbers. Equa- 
tion (15) becomes 

{1 — 2pe) [(1 + Dcot De)? + (De)?] =1. (17) 
This equation may be solved by successive 
approximations by observing that ¢ is a quanti- 
ty cf the order & (it is equal to « at the crest 
where x =o), and assuming that the 2m" deriva- 
tive, 62”, is o(æ"+1), It is evident a posteriori 
that such a solution exists provided the succes- 
sive approximations converge. Assuming this 


we begin by retaining in (17) terms of the 
order a2, i.e., 


2e(1 -p)+e?(1 -4p) ht. (a8). (18) 
Integrating | 

De)? 63 

ZT ~p) += (1 ~ 4p) +0(a4), (19) 


where the constant of integration vanishes by 
the condition that De=e=o at |x|= 00. The 


requires 


LL F2 1 + a +o (xt). 
P 


(20) 


Equation (20) is the Rayleigh-Boussinesq 
approximation. | 
To continue the approximation we substi- 

tute (20) into (19) to obtain 
(De)? = 3 ea — 3 83 + 0 (x). 


(21) 


It is now useful to express equation (17) by 


solving for D cot D-e and expanding all 
terms in series. This yields 


e) 64) a ¢(6) g(8) 

= = 
3 45 945 4725 
22 303 2 2 

= de TRE mine poe + 

2 2 2 2 

4,4 24202 565 

+ er „Der „ar +0(a®), (22) 


where terms of o(xf) are omitted. For the 
next approximation we retain terms of o(œÿ) 
in (22), substituting p=1 -& + bga?. 


g(2) e(4) 
= € e(1 — &« + b,æ?) - 
3 45 ( 2%) 
3 5 (De)? 
ee - 2a) TMS +o(a4), (23) 


We note that we can substitute for (De)? in 
(23) by using (21), and commit only the same 
error, o(a!). Furthermore, differentiating (21) 
we obtain 

135 


= ae? + =o +o(xt), (24) 


e(4) Be One 3 


which may be used to eliminate e in (23). 
We now have a polynomial on the right of 
(23). Integrating, the condition at the crest 
yields b, =21/20 or 


A À 
Pols + 0 (a), 


in agreement with Weinstein’s corrected results. 
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The above development illustrates com- 
pletely the approximation procedure. At each 
stage all derivatives in (22) except e may be 
eliminated by using the results of the previous 
approximation. The resulting equation may 
then be integrated by quadratures. Each step 
is increasingly complicated but the difficulties 
in proceeding to higher approximations are 
merely numerical in nature. The method has 
been extended to obtair the « and a! coef- 
ficients in the series for p or F2. The results 
are contained in the following equations: 


(Dore (sa as + Soe ern at) + 


2 I 400 
26 
iG 3 De ar +2 997,8) + 
70 


8 D 
vet (84 Ta HE at) + 
43.2 40 
DE 
5 25 
1509 
+ ¢6 EE |) de 7 
€ ( Re) ox); (26) 
ey 3 083 
= —- De , 34 mee 5 
ir 35” 2 800 me) 
(27) 
3 


Re ni Sr (oe) 
20 70 175 
(28) 


All quantities in (27) and (28) are obtained as 
by-products of the investigation in the next 
section. This provides a check on the accuracy 
of these two equations. 

We note that the coefficient of e? in equation 
(26) must be positive. This may be used to 
derive a maximum value of «. To the order of 
the Weinstein approximation, the first two 
terms yield &max =.800. Equating the full coef- 
ficient to Zero, &max=-813. These values are 
slightly larger than that of McCowan (max = 
=.78) and slightly smaller than that of APTÉ 


(1953) (max = 844). 

3. Use of momentum theorem for additional 
approximation 
The momentum equation is 


Ys 


J (p’ + ou?) dy = constant, (29) 
0 
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where p’ is the fluid pressure, u is the horizontal 
velocity and @ the density. The constant in 
(29) may be evaluated at a great distance from 
the crest where the pressure is hydrostatic, 
ys=h and u=c. Then, substituting for p’ from 
the Bernoulli equation, 


Ge) 


Be Frege. (30) 


and integrating the known terms in (29), 


Ys 


Bs —v?) dy=g(y,- h)?-22 (y,- h) +e ys. 
(31) 


In particular at the crest y„=h+a, v=o. In 
addition we may express # along the vertical 
at the crest as 


u=u+w penal, (32) 
; h+a 
Hence 
h+a 
J'wdy =0, (33) 


from continuity considerations. Equation (31) 
becomes at the crest 


1+% 


w? a 
a m 
[= pi A I +® 


0 


(34) 


We may express the perturbation velocity 
along the vertical under the crest as 


Ww y 2 y 4 
—=d,t4 a) 
pr (+) (5) 


in which ag, a,,... are proportional to the 
2nd, 4th,... derivatives of u with respect to 
y (or x) at the point (x=o, y=o). Consistent 
with our original assumption regarding de- 
rivatives of €, dyn =0 (@"*1)(n>o). If we neglect 


tn. MCE) 


Ay, Ag... we May compute a, and a, from the 
conditions 
1+& } 
J'wdy =0 (36) 
0 
} 
We wa ! 
=—- =(I -2pa) - 
CR CET + ( pe) a 


where w, and u, in (37) are evaluated at the 
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point (x=0, y=h+a). Obviously (37) equals 
- 02/2 + 0 (0). This and (36) lead to 
a? 
dp = ae +0 (a), 
(38) 


de = x22 +0 (a3), 
Substituting in (34) we find 


(39) 


p= rat a? Fo (dr). 


This constitutes a simple independent deriva- 
tion of the Weinstein approximation. 

We will now show that (34) and (35) may 
be used to derive one higher approximation 
for p or F? than was obtained by the successive 
approximation method of the previous section. 
Thus we have from (11) 


ax ,, 
ne mt 0), (40) 
or 
u, I 
c Tee 7 (41) 


where u, is the velocity at (x=0); =o) and 7, 
is evaluated at x=o'=0. Ran derivatives 


of (41) at (v=0, y=0) ar 
20). = = 


92 
dx" Ê = dy" 


24 
ve (1 + a (42) 
94 3E di 2) = 10 none ne 
ax’ \¢ dy"\ ¢ “ +1)? (I +%0) 
4! a 
7 ee (43) 


The constants dy, da, a, in (35) may be com- 
puted from (41) - (43) to order «> by using 
the final approximation of the previous section 
and substituting for 7 from (14): 


Qn = cosec D+ € 
(2) 7 


+ (4) + 
6 360° 


DICO CEE oe 
: (44) 
rn ST Non 1217 

I§ 120 604 000 


ed) + 09 (aS). 
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Equations (36) and (37) then determine a, and 
a, to this order. The result is 


Das) - 
er Ma : +o(). (45) 


Substituting in (34) we determine p to the 
order &. The method verifies the previous 
results and adds the « term. The final result is 


6 
Pap ga ees gee 
2 70 EIS 
2427 à A 
77 000 ” 62 (46) 


4. Further applications of the momentum 
theorem 


The momentum theorem, equation (34) is 
useful to obtain certain limits on the value of 
F? for any value of «. For example, the left- 
hand side of (34) is positive. The resulting 
inequality can be written 

F<rre (47) 
an upper limit to F? derived originally by 
STARR (1947). This result is in harmony with 
the negative coefficients of terms of greater 
order than in equation (46). (We may guess 
that the remaining coefficients are all negative.) 
From equation (7), evaluated at the crest, 
F?>2x. Combining this with (47), 

PSRSEUNERE T, (48) 
also derived by Starr. 

Obviously (48) are weak conditions for the 
speed and height of the maximum wave, since 
the condition that the integral over w?/c? be 
positive is a weak condition for a wave of 
large amplitude. This observation is in accord 
with independent indications that the maxi- 
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mum « is about .80. Setting F? = 2%nax in (46) 
we obtain &max =.801. 

It is of interest to note that (48) holds for 
the maximum elevation and Froude number of 
any steady wave disturbance that vanishes as 
x—>- oo. In particular it may be applied to 
the so-called “undular jump” observed in open 
channel flow (ROUSE, 1938) in which a super- 
critical flow (F?> 1) with horizontal free sur- 
face changes downstream to a series of waves 
about an increased mean level. Such a system 
is observed without visible turbulence only 
when F? < 2. Above this, the waves break to 
form a turbulent hydraulic jump. This connec- 
tion between experiment and theory may not 
be real, however, for BENJAMIN and LIGHTHILL 
(1954) believe that the undular type of jump 
must always involve some energy loss. 

A lower limit on the wave speed is obtained 
in the following manner: we assume the 
existence of a unique solitary wave for a given 
value of a. The symmetry of the problem then 
requires that u be symmetrical about x=o 
(the vertical at the crest), v antisymmetrical 
about x =o, and v=o at x=o. Again denoting 
quantities at the free surface by the suffix “‘s” 
we have v,=q, sin®, where © is the angle of 
the surface to the horizontal. Assuming y, 
decreases monotonically from the crest to |x| 
= co, O=0 left of the crest and O <0 right of 
the crest. Therefore v,=o in the interval 
(—co<x <0), v,=o in (o<x< oo). Since 
v=o at the bottom of the channel and along 
x=0, v(x, y)=o everywhere left of x=o, 
v(x, y)<o everywhere right of x=o; for 
v(x, y) is a harmonic function and cannot be 
a maximum or minimum within the fluid. 
We have therefore shown that 


’ 


(49) 


At the free surface 
du à I (4, 7) [2 dy ou 7] 
— | =q2— =} — iso 
a 1 (9, y). "Loy dpa], (50) 
Since u-iv and x +iy are functions of P+iy, 
Ou F dx dv, all 
le = ; si 
st) ds dy. dy, dp; dy, (sx) 
Using u,=4;cos 9, v= 4; sin © and the pressure 


equation at the free surface, 
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) 2 
=) = 2 sin 0 cos 0 ( +828), (52) 


ox R 2 cos © 
where r is the radius of curvature of the curve 
of the free surface, given by 


1 dO 
Bey (53) 
In the interval (— co < x < 0) the derivative 
of u with respect to x (52) will be negative if 
r is positive. If r is negative (in the vicinity 
of the crest) (du/dx), will remain negative left 
of the crest unless q;/|r| > g cos 2 O/2 cos ©. 
Obviously this cannot occur unless the wave 
has a large amplitude. However, for the maxi- 
mum wave, in the vicinity of the crest and 
along the vertical below the crest, q?=g(h+ 
+a-y), according to Stokes (LAMB 1932). 
The radius of curvature of a streamline is 
given by the vorticity equation 


og 4 

ee (54) 
This leads to 

qe _£ 

ia Gy 


Thus the quantity in parenthesis in (52) falls 
to zero at the crest in the maximum wave. 
Elsewhere along the surface of the maximum 
wave r=o and (du/dx), is negative left of the 
crest and positive right of the crest. It seems 
probable that for smaller waves q;/r+g cos 
20/2 cos © is everywhere positive. Assuming 
this 

+ <oin(-æ<x<o) 

— < — O0 © A = > 

abel 
9 (56) 
*\zoin (0 <x <0) 
OX ee a | 


In addition Zu/0x=0 at x=0 from symmetry 
considerations, and du/0x= -dv/dy=<o at 
y=0 left of the crest, 24/90x=0 at y=o right 
of the crest. The harmonic nature of du/dx 
then insures that 0u/0x <o everywhere left of 
x=0, duJox=0 everywhere right of x=o. 
Therefore 

Ou Au 


—>0 or — <0 at X=0. 
Pe dy? 


(57) 
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Finally we can also show that the derivative 
of u with respect to y is actually negative for 
(o<ysh+a) at x=o: Suppose du/dy=o at 
a point y, in this interval. Then, either 4/dy =0 
in the interval o<y=y, or du/dy is negative at 
a point Yo (0 < Yo < yx). If u/ay =0ino=y=y, 
the derivatives of u with respect to y all vanish 
at the point (0,0). It is easily shown that this 
would require #= const, v=o in the whole 
field of flow. This is impossible for a finite 
wave. On the other hand if du/Ay <o at yo, 
the mean value theorem shows that d*u/Ay? > 0 
somewhere in 0 < yy) = y, contrary to (57). 

These results require a velocity profile of 
the form shown in Fig. 2. To the order of the 
Weinstein approximation, equations (35) and 
(38) define a parabola. 


Fig. 2. Form of velocity profile u(y) along vertical 
at crest. 


Referring to Fig. 2 it is evident that 


+a h+a 
h+a-y\? 
ford swe f(x) = 
y* ye 
aoe Wa) (58) 


Since y* = (h+a)/2, 


h+a 


2 
fu < = (h +a). 


y* 


(59) 


To compute the remaining contribution to the 
integral over w?, we define first 


aa 
y= fwdy, À = fwdy = fdy. (60) 
0 0 0 


If w, is the perturbation velocity at the point 
(0,0), 


* A 
Fwedy = fwdy <w,A, (Gr) 
0 0 
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since w, is the maximum value of w in the 
integrand. Furthermore 


Ete Ting pense 
2 2 


(62) 
so that 


* A? h + — y* 
Judy <wA< a Swe an < 


(63) 


Combining (59) and (63) we obtain 


h+a 


/w?dy < 5/12 w2(h +a). (64) 
0 


Using equation (34), inequality (64) leads to 


(65) 


B+C 


ne 
= D” 


where 
B= 220 +610 ax +764 a? + 374 a 


"a 
C= 120(1 + (x ~ La) (66) 


D = 340 +340& +289? 


For « greater than 6/11 we have F?>2x. We 
conclude that for the maximum wave 


Fy? = 1.0909, Om=.5454. (67) 


This is slightly stronger than a result by STARR 
(1947), &,=.5000. Computations based on (65) 
and (66) are shown in Table I, and the corres- 
ponding curve in Fig. 3. 

The numerical results of this paper as shown 
in Table I specify the speed of a solitary wave 
with extremely little uncertainty for values of 
& less than .30 and with reasonable accuracy 
up to .45. Above this the limits leave much 
to be desired, although further considerations 
of the general properties of the solitary wave 
would probably permit considerable improve- 
ment. For example, it seems quite likely that 
an upper limit to the value of the integral in 
(64) is given by a parabolic velocity profile 
for w(y). If so, the coefficient in (64) becomes 
1/5 instead of s/r2. If this could be proved, a 
lower limit is obtained up to «=5/7 instead 


Tellus VIII (1956), 4 


SOLITARY WAVES IN THE ONE- AND TWO-FLUID SYSTEMS 467 
Table 1. 

Es aie Canin D On 0000 .0500 1000 .1500 .2000 .2500 .3000 .3500 -4000 
= (EQENO) matin. sc 1.0000 1.0499 1.0985 1.1488 1.1976 1.2461 1.2940 1.3415 1.3880 
u > (Eq. 66).... 1.0000 1.0497 1.0988 1.1471 1.1943 1.2403 1.2845 1.3291 1.3649 
= (gerd) lean 11.0008 1.0500 1.1000 1.1500 1.2000 1.2500 1.3000 1.3500 1.4000 
N Oe oe oe -4500 .5000 5500 .6000 -6500 -7000 -7500 -8000 -8500 
FA(Eq: 6)....... 1.4340 1.4790 1.5288 1.5657 1.6073 1.6473 1.6854 1.7217 1.7557 
A= (ERGs 66.)22. 1.3984 1.4193 1.1000 1.2000 1.3000 1.4000 1.5000 1.6000 1.7000 
P= (Kgs Amann 1.4500 1.5000 1.5500 1.6000 1.6500 1.7000 1.7500 1.8000 1.8500 
a phenomena of this type may occur in the 


| 


re) a 
180 200 210 


100 110 120 130 140 150 160 
F 


170 180 


Fig. 3. Curve (a) is the lower limit of the wave 
speed, (b) is the curve of equation (5). 


of 6/11 and the uncertainty in the speed is 
greatly reduced. Part of the proof may involve 
showing that all derivatives of w with respect 
to y are negative at the point (0,0). 


Part II 


THE SOLITARY WAVEIN À TWO- 
FLUID SYSTEM 


1. Equations of the two-layer model 


The mathematical techniques of Part I may 
be applied to a study of a solitary wave at the 
interface of two superimposed layers of immis- 
cible liquids of different density. Such waves 
have been observed in experiments by KEULE- 
GAN (1953) and Lone (1955). Recently Ax- 
DULLAH (1955) has offered some evidence that 
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atmosphere. Keulegan! derived the first ap- 
proximation to the speed and form of the wave, 
using the approach of Boussinesq (1871). In 
the present analysis a second approximation is 
obtained. This permits a discussion of the 
maximum amplitude of the wave. 

The model adopted for this investigation is 
shown in Fig. 4. The two fluids are in a channel 
of infinite length, bounded by two rigid planes 
at y=o and y=H. We assume the existence 
of a wave of permanent form with a single 
crest (or trough) at the interface of two 
immiscible superimposed liquids of density 
0, and 05. The height of the crest above the 
horizontal bottom is h+a. The fluids are in- 
compressible and frictionless. 

The dynamic equation at the interface is 


01 (Gi — 1?) — 02 (@ — 12) + 
+2(01 -— 02)g (ys —h) =0, (1) 


where q is the fluid speed, g is gravity, and y, 
is the height of the interface above y=o. We 
choose solutions for the velocity potential and 
streamfunctions similar in form to those of 
Parti 


; ye 4 
qa MO MO ET, (2) 


3 5 
MMOYMOT HMO... (3) 


à 4 


(4) 


y _n_n@ Wo), yw YB 
4 2! 4! 


1 Keulegan’s derivation was for the special case of 
Ao/[o small and u, = uy. 


Fig. 4. Theoretical model of a two-fluid system. 


— H)3 
y, = NO @-H)- No NE aL 


fe ar oars, Nr 


+ NSW 
5! 


(5) 


where M, N are functions of x. These solu- 
tions satisfy the kinematic conditions at y=o 
and y=H. The kinematic condition at the 
interface yields 
Vs 
par el nr ds 


(6) 


re (7) 
Defining the quantities 
Ys BER? 4 M® a 
h wa ; uy We 
NO ii 
i wry ‘; EL (8) 


equations (6) and (7) may be written 


n® 4 (4) : 
fe Ne healt se) ge 00 CE 
(9) 
=(2) 
C(1-r)= -e-lert ul r+e)s - 
EU) 
PPT OA LS (10) 


where superscripts now denote differentiation 
with respect to x/h. We observe that & is a 
quantity of order «= a/h (it is equal to « at the 
crest). We assume that the n'" derivative of e, 
n 

n and & are of order x 2. With this assump- 
tion we may solve for n and ¢ from (9) and 
(10) by successive approximations: 
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.e(2) gi)? 7 
An Gi ae (4) 
n= -e+te-e Be m + 
+0 (at), (11) 
DU &2 2 &3 8 ea 
Gin "Gn Gr eee 
ees (1 —1)8+0(a4), (12) 
3 > 460 


where terms of order «* are neglected. We 
now compute q? and q in equation (1) b 
differentiating equations (2) and (4) and sub- 
stituting from (8). We obtain for example 


2 
i I=2H+7?2 —7® (1 +26) + 
1 
(4) 
Zn +42) 0 (a) (13) 


Using (11) we may eliminate 7 in (13). This 
procedure yields 


2 Qe 
Z-I= SP ae ME des 


2 2 
+— €) - eel) + 0 (at), 
45 


(14) 


2 av 280 x > 483 5 

u3 I-r (I-r? (I-r® 
2. ea) 2 

+ =e (1 —r) - —+—e (1 - r)8 — 
3 ( ) 3 45 


(15) | 


= = ce) + 0 (a4). 


Substituting in equation (1) 


à ts 
ou? | — 26+ 36% — 4e + —e(02) - = QG) + 
3 3 


2 2 2€ 
+ — e() - = eet) | - es - + 
45 3 1-r 


3€? 


(os 


+260 (x = 7) — Let + 
3 3 


46 
(ae 


2 
+— e() (1 -r}8 - + et | + 
45 3 


+2 (01 — 0,)ghe = 0 (af). (16) 
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2. First approximation 


We solve (16) by retaining first, terms of 
order &?. The simplified differential equation 
may be integrated by quadratures by using 
e instead of x/h as the independent variable. 
If we define 
(17) 


An = QU = Oats (x a rye 
this integral may be written 


OS. @ he h 
dB = et — 268 — (0, — 04) g — et + o(at). 
1 


3 ay ay 
(18) 


In (18) the constant of integration vanishes by 
the condition that e® =¢=0 at |x| =. Defin- 


me 


F2 u F uz 
4 ue 02) h A > 02) h 
1 1 
bor 0, 
R= a (19) 


the condition at the crest, e®)=o, e=«, yields 

sR aR 
al a) + 0 (a2). 
(20) 


Since the right hand side of (18) must be 
positive 


F?(1 — a) + 


sR 


: ER 
ile. 


Subtracting (20) from (21) we find 


) =i. (or) 


F?(1 —-e)+ 


R? F2 
E>O if Fi Rp > 
o if F sR? F2 fe 
Es 1 = (x = R)? < 


These are the criteria for the wave to be one of 
elevation or of depression. If the density 
difference is small s~ 1. If, in addition, the wave 
advances at a speed c on the interface of two 
fluids at rest at infinity, #,=u,=cand FR, =F,= 
=F, Then, from (20) 


F?=(1 - R) E | +0(a?). (23) 
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It this case (22) becomes: 


: I 
SO he = 
2 


(24) 


: I 
eect) he 
2 


Relations (23) and (24) were derived by 
KEULEGAN (1953). 

In the general case (22) shows that the 
existence of a wave of elevation is favored by a 
large density difference, high velocity in the 
upper layer in the direction of wave propaga- 
tion, and by a shallow lower layer. If the density 
of the upper fluid is zero, only a wave of eleva- 
tion is possible and the wave speed is given by 


(25) 


in agreement with the Rayleigh-Boussinesq 
approximation. Equation (25) also holds if 
the upper fluid is at rest. 

If h is finite and the total depth H of the 
fluids is infinite, (19) becomes 

F}=1+a+0(«?) (26) 

Thus a wave on a very shallow lower layer 
behaves like a water wave in a gravity field 
(01 - 95) /0ı- | ; 

Equation (18) may be integrated again to 
yield the profile of the wave: 


e = à sech? ae (27) 


R) 


4 ets 
as R 
u sF2R? | 


sa Aa: 


(28) 


3. Second approximation 


To continue the approximation we write 
(18) as 


(29) 


EU) = 3 (e2x — &?) + 0 (x). 
1 


Differentiating 


ay 


Using (29) and (30) we may eliminate all 
derivatives in (16) (except where &? occurs 
alone) and commit only the same error o(«f). 
The resulting equation may be integrated by 
quadratures. This yields 


eu) 2a 
4 — = a; a, + a_3e4 - — a, 
3 45 4 
9 135 pa 
2 geet — 43 gq, 4.235 ga) 4 M2 (gay, — et) — 
2 2 8 a, 


— (01 — @2) ghe? + o (a). (31) 


The condition at the crest yields 


DEAR pas oe | 
PC aa Fe (TR) 
Bi Ds Re | is a] 
+a oy RFF 050 Fa Rs 


= 3 
Bas, Ar en 
lt s=o, and’u, =u,=c, we find 
Cc 2 
Sa es 3 
A I + 3 + Oo (a), (33) 


which agrees with the results of Part I. 
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Fig. 5. Limits on wave amplitude. 


Limits to the wave amplitude may be obtain- 
ed by noting that to the present order of 
approximation the right-hand side of (31) must 
be positive. This may be written 


In the case s= 1, F?= F3 this leads to 


| Ob laa = «| = 


ARGeR 
3 


The curves of this inequality are shown in 
Figs: 
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Abstract 


By means of thermistors the temperature field is measured in cold air from the sea entering a 
strongly heated arid surface. The temperature at 1 m level at first increases and then decreases 
when the air moves further inland. From these measurements the field of the coefficients of 
eddy conductiv ty is calculated. The result shows that the eddy conductivity not only increases 
upwards but also in the horizontal direction being 3—4 times larger half a km from the shoreline 
than at the sea. In the Appendix a formula is deduced for determining the coefficient of eddy 
conductivity. In doing this the definition of Az is somewhat modified considering the fact that 
the scales of temperature and potential temperature may differ, the difference depending upon 
the reference pressure of the potential temperature. The error of neglecting the effect of the 
reference pressure of the potential temperature in the definition of A; is discussed. 


I. Introduction 


In local meteorological phenomena the eddy 
conductivity plays an important rôle. It depends 
on the distribution of local obstacles, on the 
existing temperature field, which may cause 
convective phenomena, and finally on the wind 
field. On the other hand the eddy conductivity 
influences the temperature and wind fields. 

In our investigation the interaction between 
the fields of eddy conductivity and temperature 
resulted in a very peculiar effect not earlier 
observed as far as the authors know. Studying 
a rather strong and cold see breeze entering a 
shore line with a hot surface several kilometers 
inside, we have found that due to the heating 
from below the temperature at the level 1 m 
above the ground at first increases as the air 
moves inland. But several hundred meters in- 
land the temperature again drops at this level. 
This fact may be due to the increase of convec- 
tion as the air moves further inland. 


If the local meteorological phenomenon to 
be investigated is small enough (of the dimen- 
sion of centimeters and decimeters or some 
few meters) simple measurements can give 
information as to the temperature, wind and 
eddy conductivity. In such cases the pheno- 
menon of transition of energy from the sur- 
face to the layers of air a few meters above 
can be studied. 


When arriving at the problem of determin- 
ing energy and eddy conductivity conditions in 
spaces of a dimension of several hundred meters 
or some few kilometers in the horizontal and 
some ten and hundred meters in the vertical, 
the problems of measuring all the quantities 
necessary for a correct determination of the 
eddy conductivity become overwhelming. The 
only possible way of facilitating the measure- 
ments in a field study of eddy conductivity is 
to find out ideal places of measurements where 
characteristic phenomena occur, and where 
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the influence of local obstacles is a minimum. 
If the study of eddy conductivity is restricted 
to local strongly marked phenomena along 
some border line the mathematical problem 
may be simplified to a two-dimensional one. 


2. Location of the measurements 


In our investigation we have tried to find 
a place along the sea shore where the problem 
of determining the eddy conductivity is to a 
great extent reduced to a two-dimensional 
problem. Fig. 1 shows the localities in the 
neighbourhood of Gudesjö at the westcoast 
of the island of Oland in the Baltic. The coast- 
line is mainly north-southerly but at our 
observational site there is a wide bay with the 
shore line running northeast-southwest. Just 
at the shore there is a 20 m wide band of flat 
stones rizing only a few decimeters above the 
water, then there is a steep stony slope rising 
5 m above the water. Then follows a wide 
sparsely vegetated arid plain, where only some 
juniper bushes grow. Four masts were placed 
in a straight line along a small road, not quite 
perpendicular to the shore line. The effective 
distance of the masts from the shore line is then 
reduced but this is partly compensated for by 
the fact that the wind was mainly southwesterly 
during the time of observation. The path of 
the wind overland from mast to mast has its 
maximum value when the wind direction is 


fo} 


230% 


big. at: 
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The purpose was to measure the temperature 
lapse rate and the variation of it as the air 
moves inland and from these values and from 
approximate wind measurements at the third 
mast in the 10 m level, some conclusions will 
be drawn as to the variation of the eddy 
conductivity with the distance from the shore. 


3. Apparatus of measurements 


The four masts were of a type, very suitable 
for meteorological observations. They were 
made by thin aluminium tubes and of a very 
small weight, thus easy to transport. They were 
constructed by the Royal Swedish Board 
of Civil Aviation for air port illumination 
purposes and were manufactured in sections 
of 5 meters, easy to put together at the place 
of observation. 

Fig. 2 shows a vertical plane through the 
four masts. In this figure the thermistors are 
marked with small dots. In all there were 16 
thermistors of type Western Electric 17 A, 
which were coupled to a Honeywell-Brown 
6-point recorder. By means of auxiliary syn- 
chronous relays it was possible to record up to 
20 points in one cycle on the recorder. On a 
special recorder also wind speed (10-minutes- 
mean) and wind direction were obtained. 

From the records, the vertical mean-tem- 
perature field has been evaluated for periods 
of one hour. In order to compute the values of 
eddy conductivity in different heights, some 


47m above 
* sea level 
O 50 190 m 


The Oland-expedition 1954. Observational site. |, = effective length. 
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No heating 


assumed obove 


this sine 


Fig. 2. Vertical section with path of integration. (See the text). 


mathematical deductions have to be made. 
These may be found in the appendix where 
the case of mean motion along the vertical 
section is treated. Some other simplifications 
are also introduced. 

In the decuction a somewhat unusual defini- 
tion of A; is introduced, this because the scale 
of potential temperature depends on the stand- 
ard pressure in question. As c, refers to ab- 
solute temperature the following definition of 
A; has to be used 


AC] 


Un Aza 


Olly 


where T/Q is some sort of a scale factor. 

U, is the vertical mean transport of energy. T 
is the absolute mean-temperature © is the 
potential mean-temperature referring to any 
standard pressure. 


4. The field of eddy conductivity 


By means of the equation (2) in the appendix 


it is possible to deduce the mean value A! at 
the heights of measurements in the space 
between the different masts. At first mean 
values of @ have to be determined from the 
series of records obtained. In general to meas- 
urements per point were obtained in an hour 
of 1 measurement. From the @-values, values of 
Fr — were obtained in the regions 
Ox dz 

between the masts. From the wind records u 


was obtained and so the integral could be 
evaluated. The integral corresponds to the 
upward energy flux divided by c,. Now some 
of the energy gained from below is transported 
above the top of the masts, especially for the 
masts at greater distance from the shore. This 
vertical transport is smaller near the shore-line 
and vanishes above a certain level (cf. fig. 2). 
As we do not exactly know the integral limit 
z=H we put the energy flux equal to zero at 
the top of mast nr ı and integrate down to 
surface level. We may then introduce an error 
and the value of A may be too small. If we 
assume the energy transport at the top of the 
first mast not to be equal to zero, a rough 
computation shows that the computed A,- 
values will be modified in such a way that if 
the energy transport at the top is assumed to 
be half of the energy transport at ground 
surface the computed values of A, become 
twice as large. However, no significant differ- 
ence could be found between the water 
temperature and the temperature at the first 
mast at the level 18 m. The value at the surface 
level of the integral is then assumed to be the 
same also between mast nr 2 and 3 and 4 as 
radiational and ground conditions were equal 
everywhere. From this integral value we 
determine the other integral values at the 
heights of measurement and obtain finally the 
A’-values. The paths of integration are shown 
in fig. 2. 

The 1st and 2nd of July 1954 were days with 
great insolation. Fig. 3 shows an exampel of 
a vertical cross-section and also the effective 
length I. (Fig. 1) from the shore line of the 
different masts. The effective length is taken 
along the direction of air flow and replaces 
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Fig. 3. Mean temperatures at the 4 masts at various heights and mean values of vertical eddy conductivity at 
different levels between the masts. The measurements are made in the period 11.30—12.30 GCT, July 1954. 


x in the equation of Az. The figures in the 
space between the masts are the evaluated 
values of A‘. Five vertical cross-sections based 
upon five hours of measurement were obtained 
in this way and Fig. 4 shows the mean-values of 
A! obtained in these five cases. The increase of 
A’, with distance from the shore is striking, 
the values being 2-4 times soo m from the 
shore in comparision with the value at a 
distance of too m from the shore. The value of 
A =0,2 cgs-units at the surface level is very 
uncertain because no accurate measurements of 
vertical temperature lapse rate were made there. 

Sources of errors. There are many sources of 
errors in the determination of the values of 
A.. The most important ones are listed below. 

1) Time-mean-values. The temperatures have 
to be determined with an accuracy of at least 
0,05° C. If the time-mean-value of temperature 
does not include a sufficiently number of 
measurements one single measured extreme 
value could damage the whole time-mean- 
value. In our case we had to use a time interval 
of more than one hour to get the effect of 


extremes eliminated. In such a long time the 
mean-conditions, especially if the wind changes, 
may cause errors in determining all the quan- 
tities in the equation. To overcome this diffi- 
culty in all cases a faster multirecorder has to be 
used. However, in our measurements wind 
conditions presumably were sufficiently stable. 

2) No measurements of vertical transport at 
the top of masts. 

The transport of energy above the top level 
of the first mast causes an error in the determi- 
nation of A}. This error has already been 
discussed and may be overcome by using some 
kind of instrument measuring the vertical eddy 
exchange at the top of the mast. Measurements 
above the top heights of the masts will also 
give information as to the correction of the 
integral value. Such measurements must be 
made by means of captive balloons or other 
air-borne means, e.g. radiosondes. However 
such measurements are very expensive and do 
not give an accurate value of the temperature 
field without numerous and systematic measure- 
ments. 
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Fig. 4. Mean values of vertical eddy conductivity between masts at different levels. 
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5 
3) The measurement of < . The measure- 


ment of the vertical temperature lapse rates 
in the lowest layers are made from graphs 
drawn by means of the measured absolute 
values at the different heights. Instead of 
measuring the temperature values separately, 
some suitable temperature-difference instru- 
ment would increase the accuracy of measuring 
the gradient. 

4) The three-dimensional problem. When the 
wind is entering in an oblique angle to the 
shore line the problem is a three-dimensional 
one. We have reduced it to a two-dimensional 
one by using the concept of effective length 
and assuring that surface conditions are uniform 
over a large area. In order to obtain accurate 
values of Az this procedure is not sufficient, 
we have to measure also the effect of transport 
of energy perpendicular to the line of the 
masts. For this purpose it would be necessary 
to use a system of double rows of masts with an 
intermediate distance of about 100—200 m. 


Appendix 
By LARS RAAB 


The equations of energy flux. Regard a non- 
movable volume in a streaming media with 
the volume V and a surrounding surface S. 
An infinitesimal volume element is denoted 
with dV and the corresponding surface element 
with dS. 

The energy gained through radiation in the 
total frequency range is per unit mass 


Tf Cee By av 
0 


where «, is the coefficient of absorption and E, 
the black-body radiation in the frequency 
interval v to v+dv. For the volume element 
dV we obtain 


odV fa, E, dv 
0 


We denote the absorbed radiant energy per 
unit mass and unit time with U, and we get 
instead of the last expression: 


oU,dV 
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The energy gained through molecular conduc- 
tion of heat is the surface integral of the heat 
current U, (positive outwards) into the volume 


— /[ U,dS 
$ 
which can be written (Green’s formula): 
7 (PUsx oUsy DUsz = 
JT ( pe Dy Pz Le 
OU; 
= [fav 


(Einstein’s summation symbol is used here.) 
The energy gained in the volume element 


dV is: 


AU. 
ne 


U, is the heat flow in the direction of the 
i-axis. The energy gained per unit time in the 
volume element dV through radiation and 
conduction is then 


The internal energy will then increase, a 
process which is influenced by the fact that the 
air expands, 


We introduce the potential temperature and 
obtain 


Td0 dT dp 
2 Od... (de one 
2Q iT do 
en gen 


The equation of energy flux is then 


T do 
PQ dt 


ee 


u -odV 


Divide with @ dV and we get the energy gain 
of unit mass 
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In case of turbulent motion the gain of energy 
is divided into a mean-gain of a certain time 
and space dimension and a gain resulting of 
the eddy motion. All quantities have to be 
divided into a time-mean-value and a devia- 
tion from it, e.g.: 


e=0+0 
9-09 +0! 
and so on 
We get 

(Son pe U 

OR IX; 
an T+T ER 

PEt On Gh? Ob Ne 
20 | , 207 
* OX; ue OX; 


The fluctuations in U, and U,; are neglected 
= US=0 
The pressure fluctuations are also neglected. 
Then we have 
TAC" 
mm tan 
To 


Furthermore we assume homogenous condi- 
tions in the direction of the y-(2)-axis. In the 
practical case described in this paper the shore- 
line is assumed to be parallel with the y-axis. 
The eddy exchange in this direction is then 
neglected and we may thus simplify the equa- 
tion to 


es Crew at 8 90 , 00° 

Br 0 . 0 a On Ot 
re) | 6 ,90 _20 

or areas oped CONGO 


ox Ox ox ox Oz 


— 90° ’ ©) yo 5 
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Repeating the mean-value-process we get: 
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The fourth term on the left hand side corre- 
sponds to the eddy vertical transport of energy. 
Quite formally it can be looked upon as a 
mechanism of the same kind as heat conduc- 
tion. The eddy transport in the direction of the 
x-axis per unit time and unit surface is formally 
written: 


Cr = GA 


— (no summation) (1) 


* Ox 
where A, is the coefficient of horizontal eddy 


conductivity. In the definition of A, the right 
side of the formula is multiplied with the 


factor 3 Often this factor is neglected which of 


course if permissible if 1 =@, that is studying 
eddy conditions at surface pressures. If the 
formula should be correct at any pressure we 
should remember that c, refers to conventional 
temperature in computing energy gain in an 
air parcel at constant pressure. The scale of © 
may not be the same as the scale of conven- 


00 
tional temperature and hence + has to be 
x 


multiplied with 7 to give a correct value of the 


eddy energy transport. Compare the equation 
expressing the energy gain in the volume 


element dV. Here the factor — is obtained 
through deduction of the formula. 
The fact that a is equal to a scale factor 


can be seen from the formula of potential 
temperature 


172008 
p Fi p* Zi 
0 
were py is the standard pressure. 


478 


If p and p, are constant we get 


eee (p = constant) 
[oO ah 

1 © de 

cae Tye 


which is the ratio between the two scales of 
temperature. 

Now we use the new definition of A, to 
compute the energy gained in the volume 
element dV: 


Tike? 20 pl) 90 
ee Jr: 


2 (+2) dV 


0z dz 


. 00 
In cour case 1s — = ©. 
dy 


If we want to have the energy gain per unit 
mass, we have to divide with odV and we get 


2|- 
Fe 
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We now make the following assumptions: 
U, =0. 


We count quite formally the molecular trans- 
port of heat into the eddy transport of energy 
and denote: 


Beirat, tind 
Cp 
AAA ORR: 
Cp 
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where À is the coefficient of molecular heat 
conduction. Furthermore we assume stationary 
conditions: 

00 £3 

GE en o and w=o 

ot 


We then get the following equation 
11 97/00) , 07 ,/20\ eo 
© E (4 =| T7 (4: 2)| ER 


Some estimates can be made of the order of 
magnitude of the two terms on the left hand 
side in this equations assuming A, being equal 
to A,. The second derivatives of © with 
respect to x and z are of the magnitudes of 
1074 and 1071 degrees C per m?. The deriva- 
tives of A, with respect to x and z are of the 
magnitudes of 1071 and Io cgs-units, obtained 
by neglecting the term with A,. Hence the 
second term is 100 times larger than the first 
one justifying the neglect of it. We thus may 
use the simplified formula: 


Integrating from the height H downward to 
the height z we obtain 


see F336 
pis JO (2) 
… 00 Ë : 
if —=o for a sufficiently high value z= H. 
dz = 


The formula (2) is then used to treat the 
problem of variation of eddy conductivity 
normal to a coastline with different heating of 
the surface air-layers. 

The measurements described are made at 
normal surface pressure. In this case T=® and 
the use of the more complicated formula (1) 
is from a practical point unnecessary. However, 
as soon as © and T are different there will be 
an error in A, if old formula of the type 


90 


* Ox 


Of a one 


are used. This error (denoted by e) will be 
equal to e=(@ -T)/@. This means that when 
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studying eddy conductivity at a pressure of 


900 mb instead of 1000 mb & = -10 4) we 
0 
have an error in A,=3 % when using the old 


formula and referring © to 1000 mb. At p= soo 
mb the corresponding error will be=15 %. 
Hence the necessity of using formula (rt) in 
studying upper air eddy conductivity becomes 
quite clear. 
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Studying only surface phenomena at normal 
pressures near the sea level (p between 900 
and 1100 mb) the error in A. =3 % may be 
neglected and the simple formula may be 
used without introducing any serious error. 
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Variations in Insolation, General Circulation and Climate’ 


By H. WEXLER, U.S. Weather Bureau, Washington, D.C. 


(Manuscript received December 19, 1955; revised May 14, 1956) 


Abstract 


Two primary causes of variation in insolation—changes in solar output and in atmospheric 
transmission — are discussed, particularly with respect to basic changes in the mid- and high- 
latitude circulation pattern, Sunspot maxima in winter, on the average, are associated with 
greater tendency toward circulation blocks over the Northern Hemisphere and favorable 
conditions for heavy snow in regions krown to have been glaciated in the past (Wisconsin) 
ice age. An explanation, based on an assumed greater insolation during sunspot maxima, is 
proposed. But summer weather conditions during sunspot maxima do not appear to favor 
preservation of the previous winter’s accumulation of snow and ice. Hence an alternative 
explanation, based on decreased insolation accompanying periods of high atmospheric turbidity 
is proposed. The unequal response of continents and oceans to the decreased insolation causes 
heavy snow blankets to be deposited in areas known to have been glaciated and the lesser 
insolation in summer helps preserve a remnant of the snow until the next winter’s contribution. 


Introduction 


The human race is poised precariously on a 
thin climatic knife-edge. Only 50° F represents 
the difference in temperature between the 
tropics and polar regions. The average tempera- 
ture of the world today is estimated to be only 
15° F warmer than that of the last ice age when 
thick ice covered large portions of North 
America and Europe. A long-continued warm- 
ing or cooling of even 5° F can seriously affect 
the agricultural economy and threaten the 
survival of a nation. In light of these figures 
the reported warming of the earth as a whole 
by 1.5° F in the past century, due mostly to 
warmer winters in the past 50 years (WILLETT, 
1950), must be considered with great serious- 
ness as an indication of a climatic change which 
if long continued could have a crucial influence 
on the future of the human race on this planet. 


1 Much of the content of this paper was presented in 
talks before staff members of the Institute of Meteorology, 
Stockholm, 8 October 1954, and the Meteorological 
Office, London, 12 October 1954. 


What causes climatic unrest? Many theories 
have been advanced, ranging from variable 
energy emitted by the sun to changes in com- 
position of the earth’s atmosphere, to drifts 
and upheavals of continents, to changes in the 
orbital eccentricity and axis of the earth. The 
last two types of explanations are ignored here 
mainly because they cannot account for recent 
climatic fluctuations. If explanations can be 
found for these recent changes, it would appear 
likely that the same causes, differing perhaps 
in degree and duration, might explain some of 
the past glacial epochs. 

The theories treated here are based on 
changes in radiation received at the earth’s 
surface caused by variations in the sun’s radia- 
tion or in passage of this radiation through the 
terrestrial atmosphere. These theories are par- 
ticularly attractive because the unequal heating 
of the earth’s surface and atmosphere is the 
basic cause of weather, and it therefore appears 
natural to associate climatic changes with 
changes in radiation reaching the ground (in- 
solation). 
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Fig. 1. Normal monthly change in poleward gradient of 
effective solar radiation; units, 10°? ly min (10° lat.)=1 
(after Harris, 1953). 


In this paper, two hypotheses and supporting 
observations will be presented describing how 
variations in insolation can cause changes in 
general circulation and climate for various 
time-scales. 


Hypothesis I — Increased insolation 


Confining our attention at first to the winter 
season in the Northern Hemisphere, it will be 
argued that stronger insolation—whether caus- 
ed by increases in solar output, atmospheric 
transmission or solar declination—will steepen 
the poleward gradient of the absorbed radia- 
tion; this in turn will increase the poleward 
temperature gradient of the continental atmos- 
phere, increase the speed of the continental 
westerlies, and thereby increase the likelihood 
of large-scale instability and breakdown of the 
westerlies in the form of meandering planetary 
waves and cut-off anticyclonic and cyclonic 
eddies in the eastern portions of the oceans 
and western portions of continents (blocking). 
The new circulation patterns will produce new 
temperature and precipitation regimes for 
these regions. 

The inereased insolations and their time- 
scales are: 


months—increased insolation after the winter 
solstice 

years —increased insolation accompanying 
the maximum of the sunspot cycle 

decades—increased insolation caused by de- 
creasing atmospheric turbidity in 
the past 40 or so years. 
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Fig. 2. Normal monthly change in sea level pressure, 
mb (after Harris, 1953). 


Months Time-Scale 


Taking up first the months time-scale, Figs. 
ı and 2 show the normal month-to-month 
poleward gradient changes of “effective” solar 
radiation absorbed by the atmosphere and 
the earth’s surface (taking into account average 
albedoes) and corresponding month-to-month 
changes in normal northern hemispheric sea- 
level pressures, each being averaged around 
its respective latitude circle (Harris, 1953). 
From December to March an increase in radia- 
tion gradient occurs at middle and higher 
latitudes which is followed a half-month to a 
month later by increases in pressure at higher 
latitudes and decreases at lower latitudes—a 
usual indication of blocking. Following the 
reasoning given earlier, it would appear that 
the springtime maximum of blocking observed 
in the North Atlantic and Pacific Oceans (Rex, 
1950) Owes its existence to the steepening of 
the poleward gradient of absorbed radiation 
over their respective up-wind continents—North 
America and Asia. In autumn the increase in 
pressure at higher latitudes with no immediate- 
ly previous increase in radiation gradient can 
be explained by the thermal formation of 
shallow polar anticyclones at high latitudes as 
winter darkness sets in—as contrasted to the 
dynamically caused deep anticyclones of late 
winter and spring. 


Years Time-Scale 

The next time-scale, that of years, will be 
considered within the context of the 11-year 
sunspot cycle. During periods of maximum 
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Fig. 4. Average temperature 

difference (°F), sunspot maxima 

minus sunspot minima, January 

1900—1939. (Patterns over the 

oceans and large portions of 

Asia are based on too few data 
to be reliable.) 
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Fig. 5. Average precipitation difference, sunspot maxima minus sunspot minima, 
expressed as percentage of the average (normal) precipitation amounts, January 
1900—1939. (Patterns over the oceans are based on too few data to be reliable.) 


sunspottedness it is known from ionospheric 
observations that some ultraviolet components 
of solar radiation become much larger than 
during sunspot minima (MENZEL, 1948). There 
is less convincing evidence that with high 
sunspot numbers, the integrated spectrum of 
solar radiation received at the bottom of the 
atmosphere is also larger than during sunspot 
minima. In absence of adequate data, it will be 
assumed here that an atmosphere of normal 
cloudiness and transmission, when irradiated 
by a spotted sun, will have a higher value of 
insolation than is the case for a less spotted sun. 
Thus, as greater or lesser amounts of insolation 
accompanying the alternations of the sunspot 
cycle are received in winter, this would result 
in increases and decreases in poleward radiation 
gradients, with resulting changes in poleward 
temperature gradients, speed of the westerlies 
and downstream blocking frequency after the 
westerlies leave the continents. Let us now see 
if, in the winter, blocking patterns are found 
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for the sunspot time-scale and also if they 
are more characteristic of sunspot maxima 
than the minima. ; 

In Fig. 3 the average pressure difference 
between Januarys of sunspot maxima and those 
of the preceding sunspot minima are shown for 
the four cycles, 1900—1939,—the maxima and 
minima periods being defined as the respective 
greatest and smallest three consecutive January 
totals of sunspot numbers for each of the cycles. 
Thus, there are 4x 3=12 maximum sunspot 
Januarys minus 4x3 =12 minimum sunspot 
Januarys which were used in constructing 
Fer 

The pattern is one of marked symmetry of 
pressure difference roughly about the Pole, 
with three positive cells of increasing magni- 
tude being met as one proceeds westward from 
eastern Siberia, and three intervening negative 
cells of increasing absolute magnitude west- 
ward from Hudson’s Bay. Each of the positive 
cells has a negative area to the south—charac- 
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teristic of blocking—with the strongest block 
found in the eastern North Pacific Ocean, the 
next strongest in the North Atlantic, and the 
weakest in western Siberia. 

The pressure differences were tested for 
statistical significance. In order to cut down on 
the amount of work only values at grid points 
of every twenty degrees of longitude around 
60° N and 30°N were used. The standard 
deviations of January sea-level pressures at 
these grid points were secured from SCHUMANN 
and Van Rooy (1951) and a normal deviate 
test was made. The test indicates that the 
pressure differences are larger than would be 
expected. For example, the probability of 
obtaining a value of 9.5 millibars at 60° N, 
140° W is only .0024, the probability for the 
value at 60° N, 120° W, is .0028, at 30° N, 
40° E, it is .0308. Of the thirty-six grid points 
tested, seventeen had normal deviates greater 
than 1.00, whereas only 11.4 are expected by 
chance. These probabilities are for values 
selected after looking at the data so they 
cannot be interpreted in the usual way. But 
even if this factor is considered, the pressure 
differences appear to be larger than expected by 
chance. 

The strong Pacific block is probably a result 
of the increased thermal westerly winds over 
the vast continental complex of Asia, Europe, 
and North Africa, where the transformation of 
increased poleward gradient of absorbed radia- 
tion into increased atmospheric temperature 
gradient should also be the largest.! The next 
largest block in the Atlantic is downwind from 
the next largest continent, North America. 
The third and smallest block over western 
Siberia is probably a resonance phenomenon 
caused by the two oceanic blocks and the 
tendency of the large-scale planetary waves to 
be three in number in the Northern Hemi- 
sphere. 

Note from Fig. 3 that the 3 blocks are almost 
120° of longitude apart. But it is interesting to 
see that the Pacific block is located 50° to 70° 
of longitude downwind of the Asiatic coast 
while the Atlantic block is only 20° to 40° 
downwind from the North American coast. 


! Because of the larger thickness of the water layer 
over which heat from the absorbed radiation is distributed, 
the ocean surface temperatures respond only very slightly 
to changes in insolation as compared to land surface 
temperatures. 
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Although the phase of the blocking wave train 
is fixed mainly by land-ocean thermal con- 
trasts, planetary wave dynamics must enter in 
an important way to cause an eastward shift of 
the entire wave train. 

The pressure-difference patterns found by 
sunspot differences can be used to predict the 
temperature-difference patterns by realizing 
that the lines of constant pressure-difference 
represent wind-difference stream-lines. In this 
way the increased advection of cold air from 
Canada southward into the United States 
should bring with it lower temperatures during 
sunspot maxima; Fig. 4 shows that during 
sunspot maxima an average January in Mon- 
tana is 8.4° F colder than one during minima. 
Likewise, the strong advection of maritime air 


into the interior of Alaska, northeastern North : 
America and most of Greenland should account” 
for the higher temperatures observed there \ 


during sunspot maxima. Following the same / 
line of reasoning, during maximum sun- | 
spottedness central Europe is warmer, northern | 
Europe colder. 

The significance test used for the pressure- 
difference chart could not be applied to the 
temperature-difference chart because standard 
deviations of monthly temperatures are not 
readily available. The temperature differences 
for seven sunspot cycles were secured for two 
cities, Havre and Duluth. A t-test was made to 
determine whether the mean of the seven 
values differed significantly from zero. The 
probabilities were .12 for Havre and .10 for 
Duluth. Since these two cities were chosen 
from among many the probabilities do not 
appear to be unusual. Perhaps seven values are 
too few to test whether the temperature 
differences are real. 

Even the precipitation-difference patterns 
shown in Fig. 5 can be explained by quite 
simple and conventional meteorological reason- 
ing from the pressure and temperature differ- 
ence patterns shown in the two preceding 
figures. For example, the drier region in the west 
coast of Canada and northwest coast of the 
United States can be explained by the increased 
“chinook” or foehn winds shown in Fig. 3, 


! In six out of seven sunspot cycles, from 1877 to 
1949, the January temperatures for North Dakota and 
Montana are lower in sunspot maxima than in sunspot 
minima. 
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and also the southward displacement of the 
usual Pacific storm tracks by the blocking in 
the Gulf of Alaska. This southward shift of 
storm tracks combined with the greater lift 
of air caused by the higher dome of polar air 
in the western and central United States also 
accounts for the much greater precipitation in 
that region, mostly as snow. From Fig. 3 it 
can be seen from the much stronger polar 
continental anticyclone over central Canada 
and the slightly stronger Bermuda anticyclone 
that an intense Polar Front would exist along 
the Appalachian Mts. The temperature differ- 
ence chart in Fig. 4 bears this out. This would 
be a region of strong cyclogenesis with heavy 
snow deposits in the Great Lakes and Lauren- 
tian Regions. 

Referring again to the circulation change 
encouraging foehn winds along the west coast 
of Canada and southeastern Alaska, such winds 
would allow more insolation, warmer and 
drier air, and decreased precipitation to prevail 
during maximum sunspottedness, encouraging 
rapid glacial recession in that region. But dur- 
ing minimum sunspottedness, the circulation- 
difference pattern would be reversed, more 
storms would move inland from the Gulf of 
Alaska, encouraging heavy snow deposits and 
thus favoring glacial growth and expansion. 

A correlation between sunspot number and 
fluctuations of termini of outlet glaciers of the 
Juneau Ice Field in southeastern Alaska has 
been set forth (LAWRENCE, 1953) which adds 
support to the deduction described in the 
preceding paragraph. Through a span of 200 
years Lawrence inferred from geological evi- 
dence that at each sunspot minimum the gla- 
ciers expanded and built end moraines. This 
cyclic fluctuation of glaciers was superimposed 
on a gradual recession which in the past 30 
years has proceeded at an accelerated rate. 

The year-by-year study of the economy 
of the Kärsa and Stor Glaciers of Northern 
Scandinavia (AHLMANN, 1953) shows a net 
loss during the 1940’s and early 1950's, but 
with temporary gains during the 1947—48, 
1948—49 seasons, which years, incidentally, 
coincide with the recent very pronounced sun- 
spot maximum. Glaciers in Glacier National 
Park, Montana, which had experienced rapid 
decreases in area and volume between 1902 and 
1940 suffered only moderate losses between 
1945 and 1950, with at least one glacier, 
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Grinnell, undergoing an increase (DIGHTMAN 
and BEATTY, 1952). 

Both of these recent Scandinavian and Mon- 
tana glacial increases agree with the sunspot- 
glacier relationships deduced by the summary 
of the four sunspot cycles shown in Figs. 3, 4 
and 5. However, according to this relationship 
the 1948—49 budget year should also have 
been a period of glacial decline in southeast 
Alaska but instead it was a period of increase 
(AHLMANN, 1953). 


Decades, Time-Scale 


Let us consider now the past half century or 
so of warming which has been particularly 
marked in the North Atlantic subarctic in 
winter. For example, in absence of sunlight the 
average January mean temperature at Spits- 
bergen from 1912 to 1937 has increased by 
24° F. From long-wave radiation reasoning 
alone, such a temperature increase would re- 
quire a secular increase of the downcoming 
atmospheric long-wave radiation by 30%. An 
increase of atmospheric carbon dioxide content 
by 10 % could not account for such an in- 
crease of atmospheric radiation. (Nor would 
it account for cooling trends observed over 
large areas; ‘such*as*the 388085" Er drop in 
January temperatures in the Great Basin of the 
United States from 1910 to 1952.) An explana- 
tion must therefore be sought in terms of 
wind changes—bringing air of more southerly 
origin into this region of extraordinarily large 
horizontal temperature gradient. 

The suggested explanation runs along famil- 
iar lines upon realization that the series of 
catastrophic volcanoes that plagued the North- 
ern Hemisphere during the 19th and early 
20th centuries (e.g., Tomboro in 1815, Kraka- 
toa in 1883, Peleé, Santa Maria, Colima in 
1902—04) terminated in 1912 with the eruption 
of Katmai. Each of these violent volcanoes 
ejected several cubic miles of finely divided 
dust into the atmosphere, which slowly settled 
out over the years. Daily solar radiation meas- 
urements, first available at the time of Krakatoa, 
showed a three-year decrease of direct solar 
radiation of 10 % to 20 % as far away as 
France. Assuming that half of the scattered 
radiation is lost to Earth, this means a drop of 
5 % to 10 % of radiation from the sun and sky. 

From 1900 to 1938, clear weather measure- 
ments of radiation from the sun and sky made 
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in the United States and Europe have shown 
an increasing trend, averaging 0.1 % per year. 
From 1915 to 1938, 19 of the 24 years have 
been above the 39-year normal (WEXLER, 
1953). Likewise, from 1926 to 1952, the annual 
average of the direct solar radiation intensity 
measured at the Smithsonian Solar Station at 
Montezuma, Chile, has increased by 2 %. 

These increases in solar radiation in recent 
decades can be interpreted either as evidence 
that the atmosphere is slowly ridding itself of 
volcanic dust and has become more transparent 
to solar radiation, or that the sun itself is 
emitting more energy. In view of the fact 
that the major rise of solar radiation intensity 
occurred after the 1912 Katmai eruption, it is 
likely that the first premise is correct. The out- 
going long-wave radiation from earth to space 
would not be affected significantly by the 
presence or absence of the small, micron-sized 
volcanic dust particles in the atmosphere. 

It is of interest to note that in Willett’s 
compilation of world temperatures (WILLETT, 
1950) the 1885—89 pentad, which follows the 
Krakatoa explosion, has the lowest winter 
temperature in the past century. Following a 
steady rise in the next three pentads there isa 
sharp temperature drop in the 1905—09 pentad, 
after the eruptions of Peleé, Santa Maria and 
Colima. This is followed by a continuous rise 
to 1940 although there appears to be lesser 
rises in the 1915—19 pentad (Katmai?) and 
the 1935—39 pentad (Chilean Andes eruptions 
of 1932?). 

The increased insolation observed in clear 
weather would tend to increase the global 
temperature over the years unless there were 
compensating increases in world cloud amount, 
which is not possible to determine since secular 
changes in cloud amount have been compiled 
for only a few stations. However, because of 
evidence that volcanic dust particles serve as 
efficient ice and condensation nuclei in the 
atmosphere (SCHAEFER, 1950) and are even 
suspected of increasing precipitation downwind 
from volcanic eruptions (IsoNo and KoMABA- 
YASI, 1954) it is likely that the continued vol- 
canic dust “fall-out” of the past decades has 
served to decrease average world cloudiness, 
thus increasing even more the global insolation. 

Applying the same reasoning used earlier, 
because of a cleaner atmosphere an increase in 
global insolation in winter would not only 
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Fig. 6. Poleward gradient of total (sun and sky) radia- 
tion against latitude for a clear atmosphere, p = 0.8, and 
for a turbid atmosphere, p = 0.6, winter and summer. 


cause an average world-wide warming, but the 
resulting increase in poleward gradients of 
insolation and continental temperatures would 
cause greater blocking in the eastern portions 
of the Pacific and Atlantic Oceans, and a 
greater meridional austausch there which, be- 
cause of the spherical shape of the earth, would 
give the largest temperature rises at higher 
latitudes, as has been actually observed in 
recent decades. 

An important link in this chain of reasoning 
is given in Fig. 6. There it is seen that if the 
transmission, p, of the atmosphere is increased, 
the poleward radiation gradient is also in- 
creased significantly in winter, but not in 
summer. 


Sunspots and Glaciation 


In this section the advantages and short- 
comings of the sunspot hypothesis of glacia- 
tion will be described under the assumptions 
that the sun has long-period variations of the 
order of glacial periods and that the sunspot 
difference patterns of pressure, temperature and 
precipitation (Figs. 3, 4 and 5) are also charac- 
teristic in geographical distribution of the 
glacial-interglacial difference patterns, but 
differ in duration and perhaps in degree. In 
this way many features of the latest of the 
Pleistocene glaciations—the Wisconsin Ice— 
can be accounted for. 

The Pacific, Atlantic, and Siberian blocks 
would provide the winter conditions favorable 
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for (1) the simultaneous growth of three major 
ice sheets, one over most of Canada and the 
adjacent United States, another over Scandinavia 
and the third over northwestern Siberia—these 
are each associated with contemporary winter 
conditions favorable for glaciation in those 
same regions (see Figs. 4 and 5); (2) the greatest 
southward extent of the ice in the United States— 
the largest block being located downwind of 
the largest continent, Asia, intensifies the flow 
of polar air from the Canadian cold-air source 
region into the United States, which combined 
with the greatly increased precipitation in the 
same region creates and nourishes the snow and 
ice fields in central and eastern North America 
and extends the southern boundary of this ice 
to latitudes far south of those found elsewhere 
(38° in the United States compared to 50° in 
Europe and 60° in Asia); (3) distribution of ice 
in Central and Southern Europe—in Central 
Europe the increased warmth brought about 
by enhanced flow of air from the ocean would 
raise the winter mean temperature to above 
freezing and inhibit formation of snow and ice 
fields except in the highlands; in particular, 
the Alpine ice growth would be favored by 
the low temperature and high precipitation 
regimes of the cut-off Mediterranean cyclones 
which usually accompany the Atlantic block; 
(4) the milder climate in central Alaska (PEW, 
1953)—the Pacific blocking eddy brings into 
the interior of Alaska warmer maritime air 
from the Pacific Ocean; (5) origin of the “Baffin 
Type” glacier caps found in Peary Land and in the 
Canadian Archipelago—the increased flow of 
maritime air from the Atlantic would favor 
increased precipitation and glaciation in those 
areas; this favorable winter condition for 
glaciation is lacking today and nourishment is 
supplied by “... immediate refreezing of 
much of the melt water of summer.” (Bairp, 
1952). These glaciers are thought to be sur- 
vivors of the glacial period before the climatic 
optimum of several thousand years ago. 

It might appear puzzling at first sight that an 
ice age could be caused by an increase in the 
sun’s energy of a few per cent; since this 
should be expected to increase world tempera- 
ture if all this increased energy could penetrate 
to the earth’s surface and become available for 
heating it and the atmosphere above. However 
the increased insolation would create a new 
circulation pattern which would promote the 
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formation of extensive and persistent cloud 
covers in middle latitudes. 

In Fig. 7 there are plotted against sine of 
the latitude the present annual cloud amount 
and an assumed cloud amount at the time of 
the Wisconsin Ice maximum, both for the 
Northern Hemisphere and both weighted 
according to the present solar radiation re- 
ceived outside the earth’s atmosphere in each 
latitude zone. In drawing this latter curve aid 
was derived from charts of present world 
cloudiness and the January sunspot maximum 
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Fig. 7. Present annual cloudiness (Cp) 

and albedo (Ap) and assumed Wisconsin 

Ice Age annual cloudiness (C,,) and albedo 
(Aw), for the Northern Hemisphere. 


minus minimum charts given in Figs. 3, 4 and 
5. Over the zone 20° to 70° latitude, comprising 
most of the areas affected by glaciers and more 
frequent cyclones, average cloudiness at the 
Wisconsin maximum was assumed to be 70 % 
compared to the present value of 50%. Be- 
tween 20° latitude and the equator the average 
cloudiness was assumed to be 46% compared 
to the present value of so %. Both of these 
changes are in agreement with the equatorward 
shift of belts of maximum and minimum 
cloudiness accompanying the equatorward 
movement of the westerlies as the glaciers 
develop. The high cloud cover at middle 
latitudes during the Wisconsin maximum 
would be increased by the warmer and moister 
equatorial air (favored by decreased cloudiness 
equatorward of 20° latitude) meeting the cold 
continental air farther north and creating more 
frequent and vigorous cyclones, with their 
extensive frontal and convective cloud forma- 
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tion. Using values of cloud and surface albedo 
given by HoucHton (1954), the Wisconsin 
cloud amount curve was transformed into an 
albedo curve, which, together with Houghton’s 
present albedo curve, is also shown in Fig. 7. 

The critical dependency of global tempera- 
ture on albedo is such that an 8° C drop would 
follow a rise in earth albedo from its present 
value of 35% to 43 % (WEXLER, 1953). It has 
been estimated that a 11° C spread represents 
the difference in world temperature from 
interglacial to glacial periods, or 8° C from the 
maximum of glaciation to the present. It 
appears however that this is an estimate of 
average temperature changes in certain glaciat- 
ed regions and not necessarily that of the earth 
as a whole since it was based mainly on observa- 
tions of recessions and advances of snow lines 
in mountainous regions once covered by 
glaciers. Other regions such as central Alaska 
were undoubtedly enjoying temperatures high- 
er than their contemporary glaciated regions 
during the Wisconsin Ice Age and perhaps 
even higher than their present day tempera- 
tures (PEWE, 1953). Taking these warmer areas 
into account would bring the world tempera- 
ture increase since the Wisconsin maximum 
down from its estimated value of 8°C to 
perhaps half that value, say 4° C. This increase 
of 4° C in the Northern Hemisphere could be 
accounted for by a decrease in albedo from 
38% to 34% corresponding to a decrease in 
cloud amount from 64 % to 54 %, asshown 
to Figs 7.7 


In order to build thick ice over the years, it _ 


is not only necessary to have an abundance of 
winter snow but the summer ablation must not 
eliminate entirely the winter accumulation of 
snow and ice. Although the sunspot difference 
charts for July show cool, cloudy conditions 
which would favor diminished ablation of the 
critical southern portions of the glaciated areas 
(WEXLER, 1953) these conditions cannot be 
explained by a reasonable atmospheric circula- 
tion model caused by increased insolation. 
There is a slim possibility, however, that a deep 
snow blanket, once laid down over a large 


1 After the above paragraph was written the author 
learned that Wundt had earlier estimated that a world 
albedo increase of 3% (from 44% to 47%) would 
account for a 4° C drop of the present world temperature 
to a value estimated to be characteristic of that at the 
last ice maximum (WUNDT, 1933). 
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area, can influence the summer climate in such 
a way as to prevent its own total disappearance. 
Since the principal causes for summer ablation 
of a snow field are turbulent transfer of heat 
downward from air to the surface and heating 
by solar radiation, the snow will tend to 
protect itself from these two wasting processes 
by various devices: high albedo of snow, crea- 
tion of a turbulence-inhibiting surface inversion 
as warmer air from the snow-free areas moves 
over the snow field, formation of stratus clouds 
and fog over the melting snow which will 
increase the albedo even more. 

Apart from this possibility, the principal 
weakness of the increased insolation hypothesis 
for glaciation described here is its inability to 
account for the cool, cloudy summers necessary 
to preserve a remnant of the preceding winter’s 


blanket of snow. 


Hypothesis II — Decreased insolation 


An alternate way to explain changes in 
general circulation and climate would be to 
assume long periods of decreased insolation 
caused, for example, by intense volcanic ac- 
tivity which loads the atmosphere with vast 
quantities of floating dust. 


Volcanoes and Glaciation 


A pall of volcanic dust coming from a single 
volcano, Krakatoa, lowered the incoming solar 
and sky radiation by an average of 5 % to 
10 % for three years. It is therefore reasonable 
to expect that during a period when volcanoes 
of the Krakatoa type might be erupting at 
yearly intervals, insolation would be reduced 
by even more than ro %. Since the long-wave 
radiation sent upward by the earth and its 
atmosphere would not be appreciably affected 
by the micron-sized dust particles, a reduction 
of 10% of the incoming solar radiation should 
eventually lower by some 7°C the mean 
temperatures of the earth’s surface and its over- 
lying atmosphere. This would encourage a 
greater proportion of solid precipitation in 
winter which, in favorable regions, would lay 
down extensive snowfields; during the shorter, 
cooler summers, these might not be entirely 
dissipated and so would provide a foundation 
for the next winter’s snow and, if long continu- 
ed, would build up thick ice-sheets. This 
reasoning is based on the assumption that the 
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present winter precipitation amounts in the 
regions known to have been glaciated would 
not be appreciably decreased by the presence 
of volcanic dust. Let’s examine this question. 

The availability of the volcanic dust as 
possible sublimation, freezing or condensation 
nuclei might appear at first glance to favor 
increased precipitation if other natural nuclei 
were below the optimum concentration. How- 
ever, the following factors would tend to 
decrease winter precipitation: (1) lesser intensity 
of the atmospheric circulation — the lowered meri- 
dional gradient of insolation caused by greater 
turbidity (smaller transmission coefficient, p, 
(see Fig. 6)) would decrease the intensity of 
the atmospheric circulation; the reduced wind 
speed would lower the evaporation of water 
vapor from the surface and the smaller meridio- 
nal motions would reduce frontal activity and 
precipitation; (2) greater stability of the atmos- 
phere—caused by greater absorption of solar 
radiation by the dust filled upper layers of the 
atmosphere and lesser absorption by the earth’s 
surface and lower layers of the atmosphere; 
(3) lowered temperature of the sea surface caused 
by decreased insolation—this would tend to de- 
crease evaporation but would take many years 
to accomplish because of the thick layers of 
the ocean required to be cooled. For example, 
if a dense volcanic dust pall should decrease 
insolation 20 % below normal, then the average 
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temperature of a 100 meter thick surface mixed 
layer of ocean would drop 2.6° C in one year; 
but if the cooling were extended through the 
average depth of the oceans, 4000 meters, the 
average temperature drop in one year would 
be only 0.06° C. 

All these factors considered together, would 
seem to indicate lesser winter precipitation 
during periods of prolonged volcanic dust palls 
over a homogeneous earth. But for the Northern 
Hemisphere of the actual earth, with its massive 
land and ocean contrasts, one of the major 
reasons — that of decreased meridional motions 
— cited above in support of decreased precipita- 
tion would no longer be valid. For, in a period 
of reduced insolation, the continental tem- 
perature would drop much more rapidly than 
the ocean temperatures and initially (that is for 
the several decades it would take for the oceans 
to cool significantly), the winter contrast in 
temperature between continent and ocean 
would increase markedly. For example, as 
shown in Fig. 8, a 20 % decrease in insola- 
tion following the summer solstice leads, after 
the usual time lag, to a 5° C decrease in the 
average temperature of the lower three kilo- 
meters of air over land compared to only 2° C 
over the oceans. The increased temperature 
contrast between continents and oceans should 
introduce fundamental changes in the winter 
circulation and weather patterns over the 
Northern Hemisphere. 

An attempt was made to construct a new 
700-mb contour chart for January during 20 % 
reduced insolation by the following method: 


a) The present day average insolation chart for 
North America was constructed, using the 
data for the U.S. from Fritz (1949) and 
data from Canada from MATEER (1955). 


b) A new chart was then constructed showing 
the insolation received at the ground during 
January when only 80% of the present 
cloudless day insolation was observed at 
latitude 40° N. 

c) The two sets of insolation isopleths were 
plotted on the same chart and arrows drawn 
to indicate the southward displacement that 
present isopleths would have to undergo 
to become coincident with those for 20% 
reduced insolation. 

d) Since it appears likely that over the snow- 
covered interior of the continent the thick- 
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ness (proportional to average temperature) 
between 1000 mb and 700 mb depends 
primarily on the insolation received (con- 
duction of heat from below being negli- 
gible), the present day normal thickness lines 
1000 to 700 mb, were displaced southward 
by exactly the same distances that the in- 
solation isolines were under 20 % reduced 
insolation. 


e) Under the assumption that the 1000 mb sur- 
face remains unchanged (a poor assumption 
but an inescapable one at present) the new 
700 mb chart was constructed by subtract- 
ing the new thickness pattern from the pre- 
sent 1000 mb normal heights, and the result 
is shown as Fig. 9. No changes in the 700 
mb pattern were made on the west coast 
because there the maritime air was assumed 
to be in equilibrium with the ocean sur- 
face whose temperature would be but 
slightly lowered during the first few years 
of reduced insolation. On the other hand 
changes in the 700 mb pattern were made 
on the east coast and beyond that to the 
mid-Atlantic Ocean because the air in this 
area, coming from the colder North Ameri- 
can continent, would definitely be colder, 
despite its travel over the warmer ocean. 


In Fig. 9 an intense cyclonic trough is found 
over the eastern coastal area, with marked 
confluence of polar and tropical air over the 
eastern half of the United States. A strong 
anticyclonic ridge is found over the eastern 
North Atlantic Ocean with a tendency for 
cyclonic trough over North Africa—a signature 
of a blocking pattern. A similar blocking pat- 
tern is observed on the west coast of North 
America. 

This essentially thermally derived pattern is 
bound to undergo dynamic modification as 
the changed flow patterns over the Eurasian 
and North American continents react on each 
other. It is probable that calculations and 
experiments now getting under way will shed 
light on the changes to be expected, but on 
the basis of qualitative reasoning it appears that 
the changes may be as follows: 


a) extension southward of the lobes of the 
polar cyclone over the North American and 
Asiatic continents, that is, the two main 
troughs in the westerlies now located near 
the east coasts of North America and Asia 
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Fig. 9. Estimated average January 700 mb contours, 
20 % reduced insolation. 


would be intensified and perhaps displaced 
slightly westward; 


b) the warm, moist air from the oceans, not 
much changed in temperature from before 
the dust pall, would encounter at the coastal 
areas deeper and colder continental air 
masses than formerly and deposit heavier 
precipitation, mostly as snow; in North 
America this would be particularly true 
along the Gulf and Atlantic Coasts although 
there would be significant over-running of 
moist Pacific air from the northwest over 
the western slope of the dome of polar air; 


c) in the United States the southern ends of 
these elongated polar troughs would be 
vulnerable to “cut-off” by the warmer 
Pacific air thus causing closed cyclones aloft 
and heavy precipitation, mostly snow, 
through the Great Lakes region and Lau- 
rentians ; 

d) the air moving around these low-latitude 
“cut-off”? lows would tend to create anti- 
cyclonic circulations to the northeast, in a 
manner typical of the “figure-eight”” con- 
stant vorticity trajectories so often seen on 
the daily synoptic charts; 


e) in North America, these dynamically creat- 
ed anticyclones would be located near the 
Baffin Strait and southern Greenland, and 
its induced influx of moist maritime air 
would produce the heavy precipitation 
regimes needed to form the glaciers in 
Baffin Land mentioned earlier; 


f) the cold northerly flow found on the east 
side of these high latitude anticyclones would 
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Fig. 10. Estimated January temperature and precipita- 
tion patterns for 20 % reduced insolation. Tempera- 
ture departure shown in lower chart, precipitation classes 


in upper. 


become moist and unstable in its passage 
over the North Sea and bring cold, snowy 
weather to the west coast of Scandinavia; 

g) from vorticity considerations the same flow 
of air would induce a westerly trough 
centered near the Urals which, by the 
transport of warmer moister air over the 
deep dome of Siberian polar air, would lead 
to heavy winter snows in northwest Siberia; 

h) again from vorticity considerations a ridge 
would be found over central Siberia and a 
deeper than normal trough on the Asiatic 
Pacific Coast; this trough would encourage 
formation of a broad ridge of high pressure 
in the eastern Pacific. 


From thermal considerations it is likely that 
both the Asiatic and North American east 
coasts would have more intense polar troughs, 
the stronger trough being found over east 
Asia. One might then wonder why this strong- 
er trough would not favor glaciation in that 
area even more so than in eastern North 
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America, although it is known that North 
American glaciation far exceeded that of east 
Asia. The answer appears to be in the un- 
obstructed flow of moist air from the Gulf of 
Mexico northeastward over the eastern United 
States where as storms develop, the moisture 
can be precipitated over areas known to have 
been glaciated. In contrast, because of moun- 
tain barriers, there is no such easy access of 
moist air from the Indian Ocean and Bay of 
Bengal into comparable latitudes in China. 

The estimated temperature and precipita- 
tion patterns for January, under 20% reduced 
insolation, as computed by the objective tech- 
niques developed by the Extended Forecast 
Section of the U.S. Weather Bureau, are 
shown in Fig. 10. 

In summer, the dust pall would not only 
reduce insolation, the westerlies would be 
displaced farther south than normal, thus elimi- 
nating the high-level summer anticyclone 
found normally over the mid-West and re- 
sponsible for accentuation of the intense heat 
found there. The winter snow blanket would 
thus have a better chance of surviving a summer 
characterized by lower temperatures, reduced 
insolation during clear weather, and increased 
cloudiness caused by increased storminess as- 
sociated with stronger westerlies displaced 
southward. Calculations indicate that a mid- 
July weather chart during 20 % reduced insola- 
tion would be similar to that of a present early 
June chart insofar as the net radiation (insola- 
tion minus outgoing radiation) accumulated 
since January is concerned. Actually since 
January and subsequent months would be 
cooler because of the reduced insolation, the 
mid-July weather chart would probably be 
similar to that of the present mid-May. 

After 25 years, say, of 20% reduced insola- 
tion caused by a dust pall maintained by wide- 
spread volcanic activity, there would have been 
enough energy loss by this means alone to have 
cooled the oceans by an average of 1.6° C 
throughout their entire depth. Throughout 
this period, the Arctic ice pack would spread 
and thicken and in winter the continents would 
remain colder than the oceans, thereby en- 
couraging the laying down of thick blankets 
of snow as described earlier. In summer, as the 
oceans become steadily colder, the contrast 
between the warm continents and cold oceans 
would increase, which would encourage mon- 
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soon-type invasions of the cooler maritime air 
deeper into the continent, with its attendant 
cloudiness, a condition favorable to summer 
preservation of the snow. 

If the volcanic activity should cease and the 


dust settle out, the winter contrast between 


continent and ocean would be considerably 
reduced. For the oceans would warm up only 
very slowly, while those portions of the 
continents not covered by snow and ice would 
respond quickly to the increased insolation. 
The continental polar troughs in the westerlies 
would be less elongated and located nearer to 
the east coasts; the westerlies would thus be- 
come more zonal, maritime air would pene- 
trate deeper into the continents, and replace- 
ment of the continental polar air by this milder 
air, plus the increased insolation, would dis- 
courage heavy snow accumulation; in North 
America, the Pacific air moving down the 
Rocky Mts. would produce chinook condi- 
tions which would add to the warmth and 
dryness of the air. Any snow that did survive 
into the summer would meet a speedy demise 
because of increased clear weather insolation, 
the emergence of the upper level continental 
anticyclone and decreased cloudiness. 


Comparison of the sunspot and volcanic 
hypothesis of glaciation 


In summary, it appears that the decreased 
insolation hypothesis of glaciation has the 
advantage over the increased insolation hypo- 
thesis, mainly because the cool, cloudy, stormy 
summers would be more favorable for pre- 
servation of the winter snow, which, under 
either hypothesis, should be equally abundant 
in the regions known to have been glaciated. 
In both hypotheses, winters would be of the 
“low index” type, characterized by atmos- 
pheric blocks and their strong meridional 
motions. 

However, in putting forth two alternate 
hypotheses — one based on increased insolation 
and the other on decreased insolation — which 
might explain heavy winter precipitation over 
areas known to have been glaciated in the 
last ice age, an apparent paradox has risen. 

The first hypothesis seeks to explain a block- 
ing pattern by an increased westerly flow 
emerging from the continents to the oceans — 
the increased zonal winds being caused by a 
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stronger meridional gradient in insolation over 
the continent. The evidence, though to a large 
extent empirical, being based on analogy with 
the increased insolation after the winter sol- 
stice, has some theoretical support in Rossby’s 
criterion for the initiation of the horizontal 
hydraulic jump (Rossby, 1950). An excellent 
synoptic case history of fast-moving westerlies 
breaking down into blocks in the eastern North 
Pacific and Atlantic Oceans from February 22 
to March 9, 1949, is described by CRESSMAN 
(1950). The number (three) and locations of 
blocking anticyclones are determined empiri- 
cally by the sunspot difference charts as illus- 
trated in Figs. 3, 4 and 5. 

The second hypothesis seeks to explain the 
initiation of a winter blocking pattern by the 
formation of intense upper troughs over the 
eastern portions of North America and Asia — 
these troughs being caused by the increased 
temperature contrast between continent and 
ocean under reduced insolation. In contrast to 
the first hypothesis, there exist no empirical 
or theoretical guide-posts which will enable a 
convincing estimate to be made of the resulting 
general circulation patterns. But it appears 
reasonable to expect that the appearance on the 
eastern coasts of deeper than “normal” troughs 
should, by vorticity considerations, cause a 
larger amplitude planetary wave pattern down- 
wind, giving rise to blocking anticyclogenesis 
in the eastern portions of the North Atlantic 
and Pacific Oceans. But it need not necessarily 
follow that the resulting blocking pattern is 
similar in all respects to that caused by increased 
insolation. 

An interesting conclusion is that if, in winter, 
both increased and decreased insolations lead 


to greater incidence of oceanic blocking pat- 


terns, then under “normal” insolation the 
blocking patterns should be less frequent. 

In summer the decreased insolation hypo- 
thesis, in contrast to that based on increased 
insolation, would call for a “high index” type 
of circulation, with its protective features for 
snow. In summary, it appears that the de- 
creased insolation hypothesis is to be preferred 
over the increased insolation hypothesis be- 
cause its induced circulation patterns can cause 
substantial winter snow cover in the right pla- 
ces but, in addition, its summer circulation 
patterns permit more snow to survive. 

The two hypotheses as outlined here are each 


Tellus VIII (1956), 4 


VARIATIONS IN INSOLATION, GENERAL CIRCULATION AND CLIMATE 493 


dependent on the unequal thermal response of 
continent and ocean to changing insolation. 
Since the Southern Hemisphere, particularly 
in middle latitudes, is essentially maritime (less 
than 6 % of the area between 30° and 65° South 
Latitude is covered by land) one might expect 
that this hemisphere would react quite dif- 
ferently to changing insolation than would 
the Northern Hemisphere. The winter-time 
low index blocking patterns characteristic of 
the Northern Hemisphere during either maxi- 
mum sunspottedness or long continued ab- 
normally high atmospheric turbidity would be 
absent in the Southern Hemisphere. Thus, the 
large continental ice sheets found in the North- 
ern Hemisphere at elevations only a few hun- 
dred feet above sea level, would be absent in 
the Southern Hemisphere. Only the highlands 
type of glacier would be found on the Southern 
Hemisphere mountains and plateaus. This 
latter type of glacier suffers its principal abla- 
tion by insolation in summer (AHLMANN, 1953). 
Assuming then that the winter-time nourish- 
ment of the highlands glaciers in the Southern 
Hemisphere remains essentially the same, the 
variations of glaciers would depend primarily 
on the insolation received in the summer. 
A lesser insolation would preserve more of the 
glacier until the next winter’s snows, thus 
causing glacial growth, while a greater insola- 
tion would waste the glacier more than could 
be replenished the following winter, causing 
glacial decline. Thus, the remarkable recession 
of highlands glaciers observed over large areas 
in the Northern and Southern Hemispheres in 
the past decades might well be accounted for 
by the increasing trend toward higher insola- 
tion noted since before the turn of the century, 
particularly after 1915 (WEXLER, 1953). 

The question might well be asked, however, 
whether the polar ice accumulations in both 
hemispheres have also undergone simultaneous 
declines. In the Arctic there is evidence (AHL- 
MANN, 1948) that from 1924 to 1942 in the area 
north of Siberia alone, the summer ice pack has 
diminished by 1,000,000 km?, or 12% of the 
area of the total Arctic ice pack. It is not known 
whether there has been a similar decline in the 
Antarctic — either in the pack ice or the inland 
ice. But botanical evidence gathered on nuna- 
taks by the N orwegian-British-Swedish Expe- 
dition of 1949—52, 65 miles inland and 170 
miles from their base at Maudheim, Queen 
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Maud Land, indicated no recent! subsidence 
of the inland ice (ScHyTT, 1955). If this is 
characteristic in general of the inland ice of 
Antarctica, then one would be forced to the 
conclusion that the recent glacial recessions, 
so characteristic of the Arctic, particularly in 
the Greenland, Iceland and Scandinavian areas, 
do not have a counterpart in the Antarctic. The 
extensive glaciological and meteorological 
observations planned in Antarctica during the 
International Geophysical Year of 1957—58 
should throw considerable light on this ques- 
tion; if it turns out that the Antarctic inland ice 
has not generally subsided in recent decades or 
centuries, this would lend support to the thesis 
propounded here, namely, that the climatic 
variation in the polar regions of the two hemi- 
spheres need not necessarily move in parallel, 
and that because of the greater land-ocean 
contrast in the Northern Hemisphere, vari- 
ations in insolation would be expected to have 
a much greater influence on the circulation of 
the Northern Hemisphere, and therefore on 
the climate over large areas in that hemisphere. 


Future changes 


What can be said of the future? Will the 
present warming trend continue? It would 
appear from the reasoning given here that after 
many years of volcanic dust “fall-out” the 
atmosphere cannot become much more trans- 
parent to solar radiation, so that from this 
viewpoint alone, there should be no further 
significant increase of insolation and tempera- 
ture. However, to add to atmospheric storage 
of volcanic dust, there have been small vol- 
canoes since 1912, such as those of 


(1) April 1932 in the Chilean Andes which 
decreased direct solar radiation at the 
Smithsonian Solar Station at Montezuma, 
Chile, 800 miles to the north, by 3% 
during the next 3 months, although there 
was no effect at Table Mt., California; 


(2) June 1951, on Fogo Island of the Cape 
Verde group which decreased insolation 
by 8 % and visibility as far away as Florida 


for 4 days; , 


1 According to a verbal communication from lichen- 
ologist George Llano, it might take as much as 500 years 
for lichens to become established on newly exposed 
Queen Maud Land nunataks. 
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(3) July 1953, Mt. Spurr (near Anchorage‘ 
Alaska) volcano, following which direct 
solar radiation values at Table Mt., Cali- 
fornia, decreased by 3 % to 5 %, and those 
at the Blue Hill Observatory, Milton, 
Massachusetts, by 5% to 10% for 6 
months, one of the two longest periods of 
radiation deficit noted at the latter place 
since the inauguration of the solar observa- 
tions 20 years previously. 

Since even relatively minor volcanoes can 
decrease the solar radiation intensities received 
at the earth’s surface by significant amounts, it 
appears likely that prolonged periods of ex- 


plosive “Krakatoa” type volcanoes could have 
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larger effects and encourage persistent glacia- 
tion in the manner and places indicated above. 

There remains the future output of the sun. 
There is lack of agreement among solar 
physicists whether the energy from the sun 
in its present stage of development can vary by 
more than one or two per cent estimated 
during current sunspot cycles. 

Thus the problem of predicting the future 
trend of climate of Planet Earth would seem 
to depend on predicting the future energy 
output of the sun and the future transmission 
of the atmosphere to solar radiation. The latter 
prediction would depend on the ability to 
foretell future volcanic eruptions of major 
intensity. 
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Abstract 


Oscillations of surface pressure which, within the observational error, occur at the same time 
over very large oceanic regions, have been described by Frolow in the tropical Indian Ocean and 
the Caribbean. Similar oscillations in the equatoral Pacific are reported and the widespread, 
presumably planetary, extension of the pressure variations emphasized. Statistical investigation 
suggests that simultaneous falls of pressure along the equator are associated with solar disturbances. 


Introduction 


The success of weather forecasts based upon 
the analysis of moving patterns in the surface 
pressure field and the explanation of the diurnal 
and semi-diurnal pressure waves by the re- 
sonance theory incline most meteorologists 
to the view that any short-term variation of 
the surface pressure at a given place can be 
entirely accounted for by the movement and 
development of “pressure systems” of various 
scales, including the diurnal and semi-diurnal 
waves and orographic perturbations. However, 
most experienced synoptic workers in the 
tropics (and, on rare occasions, in high lati- 
tudes) have noticed curious pressure oscillations 
at the surface which are not explicable in 
advective terms; over very wide regions the 
pattern may remain the same for several days 
but the absolute value of the pressure may rise 
and fall several millibars during this time. The 

1 Present address: c/o The Meteorology Group, Army 
Electronic Proving Group, Fort Huachuca, Arizona. 

The research described here was carried out under 
Contract No. AF 19(604)-546 between the University 
of California and the Geophysics Research Directorate, 
Air Force Research and Development Command, United 
States Air Force. 
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most extensive investigations of these pressure 
changes seem to be attributable to FroLow 
(1942, 1951); he described some remarkable 
and wide-spread stationary oscillations in the 
Caribbean and the western Indian Ocean. 
Recently, Pyre (1954) described a similar 
change in surface pressure over the Pacific; 
during his study of the diurnal pressure varia- 
tion over Oahu, he had occasion to refer to the 
synoptic trend of pressure at four stations on 
the island. He says, “All four synoptic trend 
graphs show the pronounced fall of pressure on 
the ısth and 16th” (of June 1953) “which 
seems not to have been associated with any 
specific feature on the synoptic maps up to 
soo mb. The sea level isobars retreated north- 
ward on these days apparently uniformly along 
the entire southern edge of the Pacific high, 
and then moved back southward again on the 
17th”. He mentions that the surface pressures 
fell not only at the Hawaiian stations but also 
at Midway and Johnston Islands. These stations 
are sufficiently far removed from one another 
and from Hawaii to make it improbable that 
the oscillation was due to a moving pressure 
system. 
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In Frolow’s investigations, more particularly 
those devoted to stationary oscillations in the 
Caribbean (FROLOwW, 1942), a rather compli- 
cated smoothing of the data was undertaken 
in an attempt to eliminate the effect of moving 
perturbations in the pressure field. In the 
Caribbean these can produce pressure varia- 
tions at a station as large as, or larger, than 
the stationary oscillations which he wished to 
study. It is possible that the smoothing process 
introduced a fictitious wave, having a period of 
about 6 days, into the pressure curves for all 
his stations. At all events his early researches 
emphasize the periodic character of the varia- 
tions, a feature not obvious in the Pacific. Since 
Pyle described an aperiodic stationary pressure 
variation and since it was revealed by a method 
of analysis quite different from Frolow’s, we 
thought it desirable to investigate the varia- 
tions in the tropical Pacific on a scale much 
greater than any previously attempted. This 
appeared all the more necessary, because Fro- 
low’s work had not attracted the attention its 
importance warranted. 

The enquiry, then, was planned to fall into 
two parts: first, the investigation of the reality 
of stationary pressure waves in the tropical 
Pacific; secondly, if the reality of the oscilla- 
tions was established in the first phase, the 
exploration of possible explanations of the 
phenomenon. 


The Data and the Method of Analysis 


Selection of the material for study was based 
on the following reasoning: if widespread 
stationary perturbations like those suggested 
by Frolow exist in the Pacific, they can best 
be discovered by eliminating, or reducing to a 
minimum, all moving and local disturbances of 
the surface pressure field. At the same time, 
treatment of the raw data, e.g. by reduction to 
sea level or smoothing by running means, 
ought to be avoided as far as possible, since 
fictitious standing waves or the magnification 
of effects due to orography might be introduced 
by such handling. The following principles of 
selection were then derived: 


(a) Only stations on atolls, where the islands 
themselves are quite tiny and not above a few 
metres in height above mean sea level, were 
selected. This eliminated orographical effects 
and reduction errors. 
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(b) Only atolls far from any continent or 
large mountainous island were selected. Two 
stations (Tarawa and Funafuti) were about 
2000 nautical miles away from the continent of 
Australia. All others were more than 2000 
miles distant from any large land mass. Thus, 
possible “monsoon” effects, such as might be 
advanced to explain Frolow’s results in the 
western Indian Ocean and the Caribbean, were 
automatically eliminated. 

(c) The diurnal and semidiurnal pressure 
waves were eliminated by using as the primary 
datum the departure of the surface pressure at 
each synoptic epoch from its mean for the 
month. For example, the observed pressure at 
Tarawa at 00.00 G.C.T. on each day of 
March 1953 had subtracted from it the mean 
00.00 G.C.T. pressure for March 1953; simi- 
larly the 06.00, 12.00, and 18.00 G. C. T. mean 
pressures for the month were subtracted from 
the individual observations at the relevant 
times. The result is what we call here the 
“pressure departure” which could be plotted 
against time throughout the month, giving, in 
the ideal graph, 124 sequential observations in a 
31 day month. Some stations reported less 
frequently than four times a day and there 
were, of course, gaps in the record at individual 
stations. Nevertheless, the method of treatment 
eliminated the diurnal and semidiurnal waves 
from the final plot even in an incomplete 
record and at the same time left the stationary 
(if any) and moving pressure perturbations 
completely unaffected. 

(d) The moving perturbations were elimi- 
nated, or reduced, by natural means. This 
resulted from selecting only stations within a 
belt lying between the parallels(approxima tely) 
10° N and 10° S and east of 165° E. It had been 
shown some time ago that the moving vortices 
and atmospheric waves of the Marshall Islands 
produce surface pressure waves of approxi- 
mately 1/10 of the amplitude of the semidiurnal 
pressure wave. One of us (OHMSTEDE 1955) 
had also shown that simultaneous pressure 
oscillations affecting the entire Marshall archi- 
pelago occurred during the fall of 1952 and 
that these were of an amplitude 10 to 20 times 
greater that those associated with passing vor- 
tices. When it is remembered that the period 
studied by Ohmstede was one of exceptional 
typhoon activity and that many of the vortices 
passing the Marshall Islands were, without 
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Fig. 1. The atoll stations of the Central Pacific. Blackened circles indicate “parallelogram’” stations. The three 
surrounded by an additional circle were used in the study of 1955 data. Open circles indicate “fringe”’ stations. 


any doubt whatever, the early stages of ty- 
phoons, the advantage of restricting the present 
study to atolls in and east of the Marshall 
Islands can readily be seen. No smoothing of 
the pressure departures was therefore under- 
taken and the effect of moving disturbances, 
such as they were, have not been eliminated 
from the graphs for 1953 and 1954. 

The use of these principles of selection gave 
the network of observing stations within the 
parallelogram on Figure 1. A latitude belt of 
17° was covered and the greatest distance 
between any two stations was that between 
Eniwetok and Takaroa (approximately 4000 
nautical miles). In order to study possible 
interference effects of moving pressure systems 
in higher latitudes, as a contrast to the changes 
within the parallelogram, additional fringe 
stations were selected north and south of the 
parallelogram; these also are shown on Fig. 1. 
The treatment of the data for these stations was 
exactly the same as for stations within the 
parallelogram. 
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The majority of the stations shown on the 
map possessed mercurial barometers and the 
standard corrections were applied. Whenever 
possible, our data were taken from the pub- 
lished Daily Weather Bulletin—Pacific Island 
Section (DIRECTOR, NEW ZEALAND METEORO- 
LOGICAL SERVICE: 1953, 1954, 1955). These 
Bulletins lacked the observations for certain 
synoptic epochs, the deficiencies differing in the 
different years—some of these deficiencies, 
together with the complete observations for the 
Marshall Islands and the fringe stations north 
of the parallelogram, were supplied from our 
teletype record, as received in Hawaii. We 
estimated that the pressure departure graphs 
were accurate to 1/10 millibar and that the 
transcription and transmission errors were less 
than 2 %. 

In this account of the investigation we are 
concerned mainly with establishing the chief 
properties of the oscillations and with exploring 
in a very preliminary way, means to explain 
them. A total of eighteen months’ data was 
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chosen for this pilot project; six months were 
chosen by lot from 1953 and the early part of 
1954—these months sufficed to establish the 
reality of the oscillations and the data for 1955 
were then examined for clues leading to 
explanatory hypotheses. 


The stationary oscillations 


The pressure departure graphs for March, 
May, July and December 1953 and for January 
and February 1954, the six months chosen by 
lot, are displayed in Fig. 2, 3, 4, 5, 6, and 7 
respectively. Christmas Island is missing from 


MARCUS ASE 


FRENCH en P Pan ER 
FRIGATE! - A . - 
‘SHOAL > e 


| + 
ee pes Oe + 
WaxE 1 — ES = = = 
| 
| | A RT: 
JOHNSTON — Te — 3 


TAKAROA 


PALMERSTON 


3. Daily pressure departures. May 1953. 


these graphs, as the station was not furnished 
with a barometer until later. The record for 
Takaroa is very broken; it is therefore omitted 
from the mean for the parallelogram; shown 
at the bottom left of each figure. 

Examination of these graphs shows that 
stationary oscillations of the type described by 
Frolow undoubtedly occur in the equatorial 
Pacific. Note, for example, the variation of the 
barometers between 21 March 1953 and the 
end of that month (Fig. 2). Taking the mean of 
the departures within the parallelogram smooths 
out the very slight moving disturbances affect- 
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ing the stations during this period and facilitates 
the study of the details of the oscillation. It is 
surprising how many of these details are re- 
peated at cach station. Even Takaroa, with its 
fragmentary record, confirms the existence of 
the disturbance. The fringe stations also con- 
firm its existence. However, more intense mov- 
ing disturbances also affect those stations, so 
that the details of the stationary variation only 
appear in the their mean (bottom right). 

The question of the horizontal extent of 
stationary oscillations now arises: are such 
changes restricted to the parallelogram and its 
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fringes, or are they more widespread? The 
answer to this question is illustrated by con- 
sidering the oscillation 25—27 May, 1953 (Fig. 
3). The surface synoptic map for 00.00 G.C.T. 
25 May 1953 is shown on Fig. 8, that for the 
corresponding hour on the 27th on Fig. 9. Both 
these maps were derived with the minimum 
analysis on our part. The North Pacific analysis 
was copied, with kind permission, from the 
routine forecast maps of the United States 
Weather Bureau, Honolulu, the South Pacific 
being derived from the published surface anal- 
ysis of the Commonwealth Weather Bureau, 
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Fig. 7. Daily pressure departures. February 1954. 


Australia. The two analyses fitted remarkably 
well in the equatorial regions, but some ana- 
lytic splicing was necessary here and there, for 
which we take responsibility. 

48-hour isallobars for the region between the 
parallels 35° N and 35° S were derived by 
graphical subtraction on the light table. The 
results are shown on Fig. 10, which reveals 
that the stationary oscillation was much more 
extensive than the observations from the 
parallelogram could indicate. We have not 
investigated the full extent of any of the oscilla- 
tions on our graphs. It is possible that they 


are more widespread than the Pacific stations 
alone can show—unfortunately the equatorial 
and sub-equatorial data in the Indian and 
Atlantic Oceans, which would illuminate this 
problem, are not available to us. 
From the study of the months chosen by lot, 
alone, one can conclude that: 
(a) stationary pressure oscillations, like those 
described by Frolow, occur in the Pacific. 
(b) the oscillations are not large, by high 
latitude standards, the isallobaric change 
being usually 1 to 3 mb per day, and 
the total change amounting to 10 mb. 
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Fig. 8. Sea-level pressure analysis. 00 G.C.T. 25 May 1953. 


(c) although the changes are not large, they 
are, along the equator, larger than any 
other pressure variations except those due 
to the semi-diurnal pressure wave. 


(d) Pressure data from the equatorial atolls, 
Tarawa, Canton Island and Christmas 
Island probably constitute the best ma- 
terial yet available for the study of the 
stationary pressure variations. 
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Explanatory hypotheses 


In considering stationary oscillations in the 
Caribbean, FROLOW (1942) says, “We have two 
basic possibilities: either these disturbances are 
produced by a terrestrial cause or they result 
from cosmic activity. In considering extra- 
terrestrial influences it is natural to contemplate 
the sun’s activity. We know that variations in 
its activity show in numerous phenomena.” 
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Fig. 9. Sea-level pressure analysis. 00 G.C.T. 27 May 1953. 


He then shows a graph in which the pressure 
variations are correlated with the relative sun- 
spot number, commenting. The first part 
of this curve seems to be correlated with the 
other diagrams; but in order to obtain this 
correlation . .. we were obliged to assume an 
advance in phase of the solar varition of about 
36 hours. But the similarity of the sunspot 
variation with the pressure variation is not 


very great, especially in July.” 


On the other hand, Frolow is aware of the 
fact that, “Pressure waves resulting from the 
uniform resonance of the atmosphere could 
not be the same as ones caused by the sun and 
connected with the varying periodicity in 
sunspots” and also, “Because they are synchro- 
nous within a large region they cannot be 
explained by simple advection of air masses.” 
In reading the entire section of Frolow’s work 
entitled “Hypothesis and Consequences” one 
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Fig. 10. Sea-level isallobars (25 May 1953—27 May 1953). 


gets the impression that he inclines toward an 
hypothesis linking solar activity with the oscilla- 
tions chiefly because other explanations, attri- 
buting the changes to resonance or the advec- 
tion of air masses, fail to account for the geo- 
graphical features of the oscillations, but that 
he is unable to offer good evidence of a causal 
connection between the solar and terrestrial 
variations. 

It is now generally recognized that the rela- 
tive sunspot number, while valuable as an 
index of solar activity in studies extending over 
many years, and involving seasonal and annual 
data (usually smoothed) is of little use in work 
correlating day to day changes in the sun’s 
activity with corresponding changes in a 
terrestrial measure. For this reason, before 
rejecting the cosmic hypothesis entirely, as we 
were inclined to do on the same grounds as 
made Frolow cautious, we thought it advisable 
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to check the most recent pressure data (1955) 
against a more modern index of solar activity. 
For the latter we chose a measure of solar radio 
noise, since some little success has attended 
efforts to correlate this parameter with changes 
in zonal shear in the upper troposphere over 
the northern Marshall Islands (PALMER 1956). 

Three equatorial atoll stations were chosen: 
Tarawa, Canton and Christmas. As before, the 
pressure departures for each station were com- 
puted, but now the mean of all pressure 
departures for each station for each Greenwich 
day was taken, giving a “mean daily pressure 
departure” at the three stations. Finally the 
average of the three mean daily pressure 
departures was taken; this we have called the 
“daily equatorial pressure departure’. It seems 
to be the best index of the stationary pressure 
oscillations in the Pacific, short of taking means 
for the entire equatorial zone from well ana- 
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Fig. 11. Chree graphs: Average equatorial pressure departure for the thirteen days before and nineteen days 
after a solar “outburst””. On the left the entire sample of so cases in 1955; on the right two sub-samples (Group 
1 shown as dashed line, Group 2 as solid line) each of 25 cases chosen by lot. 


lyzed pressure maps, a proceeding not only 
laborious but also open to the objection that 
judgement on the part of the analyst would be 
involved. 

Days upon which “solar outbursts” occur 
were reported as a routine measure by the 
High Altitude Observatory and the National 
Bureau of Standards, Boulder, Colorado, 
U.S.A. between 14 January 1955 and 13 Octo- 
ber 1955 in the “Preliminary Report of Solar 
Activity” (TROTTER, 1955). 

We tabulated the days on which the report 
specifically mentioned, in plain language, that 
solar outbursts or a series of bursts or a noise 
storm had occurred at Boulder on one or other 
of the three frequencies 167, 200 or 460 mc/s. 
The reports of solar outbursts ceased to appear 
in the Report after 13 October 1955; for the 
period 14 October 1955 to 31 December 1955, 
therefore, we depended upon teletype reports 
of solar radio noise disturbances transmitted to 
us in Hawaii by courtesy of the National 
Bureau of Standards. While these messages 


enabled us to complete our list of “outburst 
days” during 1955, the list is not quite homo- 
geneous, since the teletype messages include 
reports from radio observatories other then 
Boulder and sometimes the information from 
Boulder is not available. While the list, there- 
fore, specifies days upon which, somewhere 
on the earth, radio telescopes observed major 
solar disturbances, it is not exhaustive; un- 
doubtedly there were other days during which 
disturbances occurred and of which we have 
no record. The important point, of course, 
is that these days were selected, not by us, but 
by independent authorities working on prob- 
lems unrelated to equatorial barometric pres- 
sure. 

The 50 outburst days of 1955 were taken as 
zero days for the application of a superposed 
epoch, or Chree, test. In this test the equatorial 
pressure departures for each of the so outburst 
days were summed and averaged; similar 
averages for the equatorial pressure departures 
for all so days one day after zero day were taken; 
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the process was repeated for a range up to 19 
days after and 13 days before the outburst 
days. The result is shown as the left-hand graph 
of Fig. 11. The graph distinctly shows a connec- 
tion between solar activity and the equatorial 
pressure departures such that the pressure is 
lower than the mean for the month after solar 
“outbursts” and higher before. 

The exhibition of graphs of this type usually 
elicits from critics the claim that the results are 
due to “chance”. It is always possible to justify 
this claim since there is always a finite proba- 
bility that sampling fluctuations can lead to a 
fictitious order in a sample, the order being 
absent from the parent universe. The question, 
then, is not whether the order is or is not due 
to chance, but rather what the probability is 
that it is so attributable. Student’s test may be 
applied to the data of Fig. 11, and it shows that 
the chances that the pressure departures of the 
sample for the days -ı to —13 inclusive and 
those of the sample o to + 19 inclusive were 
drawn from the same (normal) universe by 
chance is less than 10°”. It might be objected 
that statistical tests based directly or indirectly 
upon the assumption that the means of the 
pressure departures form a normal universe, 
and are independent, are unjustified, since the 
pressure departure on any given day may be 
expected to show a high positive correlation 
with that of the day before and the day after. 
But, far from vitiating the conclusions to be 
drawn from Fig. 11, the argument from the 
persistence of the pressure departure enormous- 
ly enhances the value of the statistical analysis. 
On the basis of persistence alone, no-one would 
expect the pronounced pressure fall that is dis- 
played as occurring between days —3 and +s. 

This matter was investigated more closely 
by dividing the so cases into two groups of 
25 each. The groups were selected by lot, the 
last two digits of the numbers in the right 
hand column on page 149 of the 1956 Edition 
of the Honolulu Telephone Directory being 
used as a set of random numbers. Chree ana- 
lysis of the two groups (Group 1 and Group 2) 
are shown on the right hand side of Fig. 11. 
Although the samples are rather small, the 
general character of the correlation between the 
pressure changes and the solar disturbances is 
preserved in each—pressure is lower than the 
monthly mean after a solar disturbance and 
higher before. 
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The significance of the pressure changes does 
not depend, however, upon this minor statisti- 
cal enquiry; if it did, it could be dismissed as 
belonging to meteorological curiosa, of which 
the literature of the past hundred years is full. 
Since the war, however, statistical investiga- 
tions of surface pressure changes in high lati- 
tudes and their correlation with a variety of 
measures of solar activity (Duet, B., and 
Due Lt, G., 1948; Craic, R. A., 1952; SHAPIRO, 
R., 1954, 1956) and the theoretical work of 
HAURWITZ (1946) and WEXIER (1950) show 
beyond doubt that there is some connection 
between disturbances on the sun and subsequent 
pressure changes at the surface. Unfortunately, 
the changes in high and middle latitudes are so 
overlaid and obscured by large interdiurnal 
changes due to travelling pressure systems that 
data gathered over many years and subject to 
laborious analysis are required to demonstrate 
the connection. In the equatorial region travell- 
ing disturbances of the pressure field have the 
least effect on the barometer. We should be 
able, then, to study the solar—induced pressure 
changes synoptically, something which has so 
far proved to be impossible in high latitudes, 
and a short run of observations is sufficient to 
delineate the chief features of the connection. 


Relation of the Surface Changes to the 
Circulation Aloft 


Clearly, much more statistical work, as well 
as synoptic analyses of individual situations, 
will be necessary before the properties of the 
solar—trophospheric relationship that lead to 
the surface pressure falls could be suggested 
even in the barest outline. At this stage, how- 
ever, it may be profitable to see what can be 
deduced from an assumption that the rela- 
tionship that holds in 1955 will be found valid 
in other years and that future statistical work 
will set the empirical connections described 
here upon a firm foundation. 

First, as has already been remarked, the rate 
of pressure fall from day zero to day five, the 
maximum, is still very small. If the upper level 
mass divergence that causes it were confined 
entirely to the layer above the tropopause 
(approximately too mb) observed measures of 
the divergence in that layer over the Marshall 
Islands during atomic tests are more than large 
enough to account for the pressure changes. 
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Fig. 12. Chree graph: Average meridional shear of the 

zonal component of the wind at 30,000 feet between 

Wake and Eniwetok for the thirteen days before and 

the nineteen days after a solar outburst. Entire sample 
of so cases. 


Thus the theoretical speculations already re- 
ferred to, which postulate rapid ozonospheric 
divergence at the sub-solar latitude during 
solar flares or other heightened emission by the 
sun, would be quite compatible with the pres- 
sure changes described here. 

On the other hand, one of us (PALMER 1953, 
1956) has suggested that there is a major 
effect of solar disturbances on the vertical 
component of the vorticity in and near the 
quasi-stationary pressure trough of the upper 
trophosphere in the Central Pacific. This trough 
varies in position and intensity with the season, 
but in general the effect of changes in its 
intensity, and in particular the formation in it 
of upper level cyclones in the layer 25,000— 
50,000 feet, can be followed by studying the 
day to day changes in the meridional variation 
of the zonal component of the wind at 30,000 
feet between Wake and Eniwetok; this para- 
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meter gives a very rough measure of the 
“cyclonicity” of the upper atmosphere in the 
northern Marshalls—Johnston Island sector. 
We have continuous measurements of the 
parameter from 1 January 1953 on and sporadic 
measurements, chiefly confined to periods dur- 
ing which tests were being carried out in the 
Marshalls, back to 1946. It is therefore a simple 
matter to apply the Chree analysis to the data 
for 1955, using the so outburst days already 
chosen in the pressure study. The results of 
the analysis are shown in Figure 12. There is a 
well-marked tendency for the shear at 30,000 
feet to become cyclonic three to eight days 
after a solar outburst. Results from other work 
indicate an average lag period of five or six 
days. 

In comparing Fig. 11 and 12, we are tempted 
to suggest that the surface pressure falls along 
the equator are a reflection of upper level 
divergence affecting first the ozonosphere and 
later the upper troposphere, culminating after 
five or six days in the first minimum of pressure 
on Fig. 11 and the attainment of cyclonic 
shear in the upper atmosphere over the northern 
Marshall Islands. The second pressure minimum 
at about plus fifteen days would then have to be 
attributed to upper level divergence associated 
with movement of the upper level cyclone 
(or cyclones) and its subsequent development 
in higher latitudes. But the pressure changes 
are so widespread that, if there is anything at 
all in this speculation, one ought to predict that 
extensive upper air data will show that the 
formation of upper level cyclones is triggered 
by solar events, not only in the Marshalls— 
Johnston Island sector but also other parts of 
the tropics between the five and fifteen degree 
parallels, in particular, in the South Pacific and 
in the region between Central America and 
the longitude of Hawaii. Upper air data from 
these regions is so far inadequate. or entirely 
lacking; it is to be hoped that the upper air 
observing program planned for the equatorial 
regions during the International Geophysical 
Year will facilitate the checking of this predic- 
tion. 
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Some Remarks on the Sudden Cosmic Ray Increase on 


Februari 23, 1956 


By E.-A. BRUNBERG and D. ECKHARTT. The Royal Institute of Technology, Stockholm 


(Manuscript received June 8, 1956) 


Abstract 


The north and south pointing GM-telescopes at Stockholm measured different cosmic ray 
increases during the solar flare on Feb. 23. It is shown that this difference is due to particles 
of momenta greater than 10 GeV/c which could not have travelled along straight line paths 


from the sun. 


A first report of the Swedish registrations of 
the large cosmic ray increase on February 23, 
1956, has been given by SANDsTROM (1956) 
and ECKHARTT (1956). This increase was 
associated with a solar flare which was observed 
to have started at 3.30 h GMT (FORBUSH, 1956). 
A few important conclusions could be drawn 
from the Stockholm registrations measured 
with GM-telescopes in different zenith angles 
in the north—south plane. These conclusions 
are supported by the results of the directional 
telescopes at Rome (BACHELET a. CONFORTO, 
1956). 

In the following discussion we have used 
the data from the unshielded four-channel 
directional telescopes at Stockholm (magn. 
latitude 2,,= 58° N). The telescope measuring 
in zenith direction will be denoted by Sthlm z. 
The north and south pointing telescopes 
measure in directions inclined 30° to the zenith 
and are denoted in an analogous way by the 
subscripts N and S. The station is situated at 


an elevation of 46 m above sea level. The 
total amount of absorber above the instru- 
ment is 1047 g/cm? which is composed of 
1039 g/cm? of air mass (equivalent to a mean 
air pressure of 1023 mb during the flare) and 
8 g/cm? of material of the roof. 


It must be pointed out that a check of the 
pulse rates of the individual channels revealed 
that one group of GM-counters changed and 
increased its sensitivity between 19 h and 20 h 
GMT on the 22nd. This group is common to 
one of the four channels for both the north 
and the zenith directions. After 20 h the 
particular group remained on the new and 
constant level of sensitivity. Consequently only 
the data after 20 h have been used. Before this 
time, the sum of the three properly working 
channels of Sthlmy has been taken and multi- 
plied by 4/3. 

The pulse rates of the different telescopes 
are given in Fig. 1. 
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Fig. 1. The uncorrected hourly pulse rates measured by the directional telescopes at Stockholm together with the 
mean air pressure during the corresponding hours. Small figures indicate opening angles of the different telescopes. 


The North—South Asymmetry 
The pulse rates from Sthlmy and Sthlms, 


uncorrected for pressure, are compared. Since 
particles arriving from these inclined directions 
have traversed atmospheric layers which have 
approximately the same structure, we assume 
that the absorption effects are the same for 
the two directions. Comparing the pulse rates 
of the two telescopes therefore means com- 
paring the flux of primaries above the at- 
mosphere. 

In Fig. 2a the ratio between the pulse rates 
of Sthlmy and Sthlms together with their 
probable errors are plotted. Dots represent 
values derived from all the four northern 
channels, crosses represent values computed 
from the three properly working channels 
times 4/3, as mentioned above. It is seen that 
before the 20—21 h interval on the 22nd there 
are differences between dots and crosses due 
to the difference in sensitivity of one group 
of counters, as mentioned above. Using the 
values corresponding to crosses before the 
particular time, we get the same value of the 
ratio before and after the flare. 

Fig. 2a shows that during the hours follow- 
ing the solar outburst more particles were 
recorded by Sthlmy than by Sthlms. The 
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Fig. 2 a. The ratio between the hourly pulse rates of the 
north and south pointing telescopes at Stockholm. 


Fig. 2 b. Thé ratio between the bi-hourly pulse rates of 
the north and south pointing telescopes at Stockholm. 


Fig. 2 c. The ratio between the bi-hourly pulse rates of 

the directional telescopes at Bergen, Norway. These 

telescopes are preferably orientated in the east—west 
plane. (See text.) 
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slight decrease immediately before the onset 
of the flare is not outside the statistical errors. 

In Fig. 2b the ratios between Sthlmy and 
Sthlms, computed from bi-hourly values of 
the pulse rates, are given which allow a 
comparison with the values measured at 
Bergen, Norway, (Am = 61.5° N), Fig. 2c. 
The Bergen telescopes are orientated in the 
ENE—WSW plane below about 28 cm of 
lead. This is essentially an orientation in the 
east—west plane (TREFALL a. TRUMPY, 1956). 
While the Stockholm bi-hourly values clearly 
show the asymmetry in the north—south 
direction, no such definite conclusions can 
be drawn from the Bergen values concerning 
the east—west direction, as has already been 


pointed out by Trumpy and Trefall. 
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Fig. 3. The ratio between the relative increases 
in pulse rates of the north and south pointing 
telescopes at Stockholm and Rome. 


In Fig. 3 we show the ratio between the 
relative increases of pulse rates for the two 
directions with respect to the mean value of 
the two hours preceding the event. The ratios 
for the Rome telescopes (Am = 42.5° N) are 
computed from the published 15-minute 
records. These telescopes are inclined 30° to 
the zenith and have slightly larger apertures 
than the Stockholm telescopes. It is clearly seen 
that during the two hours after the beginning 
of the event, relatively more particles are 
measured by the north pointing telescope af 
both the stations. 

Discussing this north—south asymmetry ond) 
must take into account the deflection of thé 
primary particles in the earth’s magnetid 
field. The initial direction of the primary 
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particles far away from the earth will be called 
the asymptotic direction in the following. In 
Fig. 4 the asymptotic directions for those 
primary particles which are measured by the 
north, south, and zenith pointing telescopes 
in Stockholm are given. Let us represent the 
asymptotic direction by a vector. The end 
point of this vector coincides with the center 
of a globe. 

Each curve on the globe is then the locus 
of the foot points of the afore-mentioned vec- 
tor, which gives the asymptotic direction for 
different particle momenta (BRUNBERG, 1956). 
The zenith angles of arrival are denoted on 
each curve (32° N, 32° S...) and the numbers 
indicate the momenta in GeV/c. As an example 
take the point 18 GeV/c on the curve 32° S. 
This point has in the geographic system the 
coordinates Y=30°+18°=48°E, ®=38°N. 
(The longitude of Stockholm is 18° E.) This 
means that a particle with a momentum of 
18 GeV/c which is approaching the earth from 
an asymptotic direction, defined by the two 
angles Ÿ and @, will arrive at Stockholm from 


„the 32° S direction. 
© Itis seen from Fig. 4 that Sthlmy and Sthlms, 
. the asymptotic directions of which are repre- 
_ sented by the curves 32° N and 32° S for pri- 


mary momenta higher than 12 GeV/c, scan 
considerably different directions in space. 
Particles with momenta around and below 
10 GeV/c, however, come from nearly the 
same directions independent of the measuring 
directions of the telescopes. This depends on 
the focusing effect of the geomagnetic field 
which has been discussed by Äsrröm (1956). 
Below 10 GeV/c the asymptotic directions of 
32° N and 32°S are still nearer each other 
than are 48° N and 48° S, which are plotted 
in the figure. Note that on the 32° N curve 
even momenta up to 20 GeV/c lie within the 
main bulk of asymptotic directions for lower 
‘momenta, whereas the 32° S curve runs away 
to the west. 

Hence, the fact that Sthlmy records more 
particles than Sthlms indicates firstly, that 


‘there are particles with momenta higher 


han 10 GeV/c in the spectrum of the extra 
lare particles. Secondly, either the intensity 
i.e. the number of particles) or the energy of 
he extra flare primaries, or both, are different 
“ the different asymptotic directions. The 
same conclusion is arrived at when con- 
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Fig. 4. Photograph of a part of the globe on which the 

angular coordinates of the asymptotic directions for the 

Stockholm telescopes are drawn. The dotted line defines 

asymptotic directions that have a 70° distance in angle 
from the direction towards the sun. 


sidering the geomagnetic deflection for the 
telescopes at Rome. 

The conclusion that there are particles with 
momenta higher than 10 GeV/c in the spec- 
trum of the solar flare primaries is supported 
by the observed increases at the equatorial 
stations of Huancayo (MEYER a. SIMPSON, 1956) 
and Ahmadabad (SARABHAI et al, 1956) 
where the cut-off momentum is around Io 
GeV/c. 

Another important conclusion that could 
be drawn from the observed north—south 
asymmetry is, that the flare primary particles 
did not travel in straight line paths from the 
sun to the earth as has usually been assumed 
(FIROR, 1954, SARABHAI et al., 1956). Consider 
Fig. 5 in which the position of the earth and 
the asymptotic directions for the Stockholm 
station are given in relation to the sun at the 
beginning of the flare. From this picture it is 
evident that particles with momenta 10 (20) 
GeV/c and measured by the Stockholm north 
pointing telescope have arrived at an angle of 
about 70° (55°) to the earth—sun direction. 
Particles with momenta higher than 10 GeV/c 
and measured by Sthlms have been deflected 
still more (compare Fig. 4). Considering the 
telescopes at Rome, one finds the angular 
distances between the asymptotic directions 
and the direction towards the sun generally 
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Sun 


holm 


Fig. 5. The angular coordinates of the asymptotic directions for arrival at Stockholm together with the direc- 
tions towards the sun at the beginning of the flare. 


smaller except for those particles of higher 
momenta that are measured by the south 
pointing telescope. 

This deflection is of particular interest as it 
shows the existence of magnetic fields in the 
interplanetary space. From the above men- 
tioned values for deflection angles and corre- 


sponding momenta, a value of 0.6 : 1075 gauss 
is computed for the mean magnetic field 
through which the particles must have passed 
in order to be measured by Sthlmy if 
emitted from the vicinity of the sun. 

The authors are indebted to Prof. H. Alfvén 
for many valuable discussions. 
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Current Data on the Chemical Composition of Air and Precipitation VII 


(For further information see EGNÉR, H., Eriksson, E., Tellus 7, pp. 134—139, and Tellus 8, p. 285) 


Precipitation April 1956 (D 604) Air April 1956 (L 604) 
mg/m* Si. g/m (= kg/km?) 
: =! [018 
Lu ed jos lia gal ee ‘ 
S Cl fe rm Na | K | Mg | Ca a x|a| S | Cl ee Na | K|Mg| Ca 
2 | 2 = a 
23 2 fe) 2 2 ce) I 6} 6.1 12 4 * * * k * * * 
25 4 o [6) 5 2 2 6| 5.7 8 677.3| 70:91 071 2711 Aol Ors" 277 
22 15 2 4| 19 3 4, 2 6.5 2 TS) 35| ©0), T7 7.00.74 ©7320 
17 124 3 5| 25 4 5| 108} 6.7| 67 53 22 1.31 2.7 2.510 3-211.5:67 5:9 
I 6 > 5 9 > ~ = * * * 525 31S 23) 208 2.200335| 69 
12 17 I 2| 18 3 Al 282 6.8. 74 IQ} 0.0| 0-0] 2.4| 2.4| 3-7] 12. | To. 
4 15 I 2) TT 7 2: 6} 6.5 “à 37| 1.7| 0-0; 0.5| 0.5] 0.7] 0.2] 0.9 
29 EF 2 30 17, ANIME MG IE T2 Lt} 2-6] 1.6| 1-0} 1.0] 1-3] ©.5| 3.9 
18 6 I Ey In fe) 37 140051 35 12 
9 13 I I 4 4 5r1Leı 75315380] 78) 0:6 02 nA] 0.83e2.517 0.21 21.0 
22 8 2 3 5 5 2 Ole 5.7 o 9 O16) 273) EA 7.52.21, 0.21 4:9 
16 9 3 6 9 4 AIRIS TI 5:8 5 2 0:0, WS 723 0.5], 22 0:8) 287 
NE 13 3 8 5 4 Sie SLs 520 Ol 3m 8:81 0:4) alll 10:6| tse. sean 
10 6 I I 6 2 3 9) 5.8 9, 171 20:0] 0.01 0:71.0:4|12.00 0.531202 
18 9 I o 9] Io Sie 500 6:81 7A 27) 2227| 9:08 2.4 0.9 0:310,0.202.7 
13 9 3 4| 13 4 aim TSH 5:9 8 20) 0/0 5.0) Ed 2131 0.00 0:0| 2357 
Io 15 Al rai 23) Tr 9| 255| 5-8} I4| 190} 22 0.0| 4.0] 0.7] 0.7) 0.7] 13 
I2 ÊT 2 8 7 3 30 22) 4:8 Ol SS) 92-4) S272) 5228) eno. seo LO 
18 8 2 31 20 5 Ale 1210 6:3] 24: 14 —| — 
10 7 3 5| 14 6 4] 37! 6.8) 67| 41! 1.0) 8.9) 3.2) 1.4] 3-3] 5.2| 14 
II SOM ny Sa APN ey Sa GST GE 55 dl QT A | Ales Sie ee ee 
13 123 5 4| 80 6| 36] 15| 5.0 o} 58| 5-9] II. ZO), A: 0:00 1.5148 
33 son 6 2; 5 Ale SO} 52 (6) wo ar Ho GO| 5 Kr ar Be 
20 7 2 3 6 Sn | Sr ©)  =25| 2.8 30:0 1.51 (0:7) 0.8.0.2 25 
34 Ey, 4 6} 20 6 SIMSOMG OT RAIN EN), — — 
54 75l 19] 40| 40 8 8| 30; 4.4 o| 38 
20 521, BEIN 300 42 6 za sai 6-0], ı7| 42 82! 0:0 6.0) 0.7.0.6 0-31 5:8 
25 16 SI, bz 9 8 5} 55} 5:9] ı0| 24 
33 20 8 250 12 6 4| 32) 4-8 ol 25 8 % “ 
2 2 Hol Tale 72 6 5| 122] 6.7| 69] 34 = 2 % es 
24 25 6} 11 9| 62 7| 26] 4-5 Ole AZ) 2:7) 0:5, 2-8) :0.3|07.017.0:0| 7e 
19 34 8} 20| 25 Si 131, 56) 5-8 8| 40 
31 45| 120 31 29) 25] Lol) Cok 5-0 o 37| 5:0] ©0| 25| 1-0], 7.8. 0.21 5.6 
8 7 4 9] 16 5 Ty) 15 5-2 2055 
QE Io 6| 10 8 6 3| 30] 4.8 Ol AO} 21:61 72 Ir SNS Sr 2 
3 42 I ı| 21 * * * * * * * * * * * * 
33 6608 170 26103651 221 447 460 44 o| Irol 0.6) 10. | 0.9} 6.5| 0.0) 12. | II 
II 9 2 4 Sr £2 Z| 266.3; 34| 23| 127 3.8) 3.0) 0.7| (0-6) 2.2| 19 
2 o I 6 5 * * * * * * * ok * * * x * 
8 3 I Te LO 5 2 18) 6.41 25 20| 6.8 0-0] 0-5} 0-4) © 5 19 2.9 
52 20 5 ZI Lak 8 Sim Gos 22 13 “ ag x * = 
15 46 OW) 2210.28 7 INTIMES fo) 2| 9.7| 6.8] 4.0] 1-4] 0.0] 0.8] 8.1 
16 60| 14] 40) AL 15 9| 134] 4.7 (6) 86] 20. 3210 5B 2.2) a 7.210. 0-3, cm 
II 7o| 38| 409] 52| 16] 15| 130| 4.4 oe 157| 7.4 23.-| 5.9) 10. 16). 0.210.907 
30 104 2 le OSs 5.7 fe) 29| 29. 4.6| 2.8] 2.0] 0.3] 2-3] 6.5 
18 50 8 Se 36... IR role mo. 297 4619 2-8 2:70 2.512.0.0|0 0:5 3.0 
52 1261, AA A2lE844 231025101360 14.5 (6) 51| 28. | 16. 428 1:30 0.01 771.71,26 
26 920.22 9 35| 25] II 156| 43 fe) 75) 2 Tie | 3281, 3.70 1.41 06-0079 
* * * * * * * * * * * LR Bat 7 il Ze 1.9 5.5 21 
89 61| 44| 119| 39] 36] 36] 290] 5.6| ro] 39| 26. | 4.6| 10 on ol Geel) 221 
119 57| 42| 68] 44| 24} 11| 136| 4.8 (6) 25| 12. 6.5| 4-3} 0.6| 0.0| 0.2| 1.3 
100 28| 481 84] 161 241 141 1421 4-5 Om 4 65 IMT-0l0.7IN0'0INonIES "x 


— No sampling 
* Sample discarded 
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$14 NOTES 


Precipitation May 1956 (D 605) 


~ 


Air May 1956 (L 605) 


2 mg/m” 3 falls ee g/m? (= kg/km?) 
s| cl | &|@|Na| K|Mg| Ca S [ela] S| a|g|Na| K | Mg| Ca 
Z| a en A 

Ri val 32) vs Es Ay) Ziele. a6)" 5g” Oh Er Ciera PR RSS 
Ki 20217 BZ I fo) 5 2 2) 50) 5-9 6| T2| 1.9| 0-0} 2.7| 0-8] 2-4) 0.01 57 
Ar 15] 12 18 2 122 2 al 321 6:7| 63| 23| 2.3] 0.0] 163) £-6| 2.1. 0-.0/04-0 
Oj 8 23 35 al Sigil), ZOO 5| 75 * *| 89] 6.0] 0.0] 2-7] 1-0] 1-4] 1-3} 8.9 
Ro 41| 36 29 Se > 6 6| 40] 6.3] 12] 14| 4-9| 9-4] 1-9] I-T] 1-1] 1.4| 8.7 
Of 728 18 2 323 8 4| 42 & #| 60] 0:0| IE | 48) 2.417 2:0] No.5) so 
Br | oll eo All ol ei. ol a8) eal © 9| 12] 3-5} 0.0] 1.0| 0.7) 0.3) 6.7| 3-0 
AF By! ads) 31 3 EC 3 2| 37 6:3| 197| 12) 2090|, 32 22] 220 5-7 OO moe 
ÄH 
Sv Ta 217 Ol) at Pole 25" ssl ei 7a) GS ton 34 
Ra PE|| a5 20 3 TES 7 3| 2724|. 6.410 306] -30| 1.2] ©o| 3-11 0.9| E.5| 0.7] 5.6 
Äm Tele 7 2 del) Mir 3 2\ 38 6.21 22/22) 71.81, 0.018521 rer 0) Odie 
Sa il 2X) 12 3 2| 15 6 8| 55| 6.6| 18] 34| 1-4] o.0| 1.6] 0.6] 0.6] 0.1] 0.0 
Ul 16| 24 2 4 Fi LO 5 2| 47| 6-2] 5] 24] 4.1] 0.5] 1-4] 0.5] 0.5) 0-5} 0.9 
Er Fa) 14 3 3] 16 7 3) 42| 5-3 ©] I6| 2.8| 2-5| 1-6| I-0| 1-5) 1-4) 0-4 
St 9728 ay 3 I] 15 7 61° 511 5-8] rz) 39] 7-5| 2-3] 2.2] 0-9) 0.61 0-8) 5-5 
Fo 19| 48 Pail 6 O22 6 4| 32] 5-8 6] 40) 15) 3:8) RSI rer 
Kvı IO} 246 55 or 251351 27 817320|75:2 o| 240| 72. 0.O|II. 2.0| 2.9] 2-9] 43- 
Kv7 &| 55 19 4 6| 14 5 3} 30] 4.6 o| 66] 12. 0.0| 4-1| I-I] I-I] 0.9| I4- 
VN E2\. 25 24 4 5| 18 3 5| 28) 6.1] 28) 36 
La a xR 16 6 ©| 16 6 4| 22| 5.8 I] 14 * ÿ rs 5 * = 
Bo nn RS aka 6 3| 81 6| 16| 45| 6.7| 96 2| 6.3) 51. | 1-8) 31. 1:.7| 13-7] 2% 
Vi I5| 42| 333 7 4) 195 9 27| 4X) 47 o| 115| ro. |280. | 2.7| 27. | 5-3} 15. | 19- 
Fa 7123 34 4 2 2721|, 270 5| 64 * #1 SE 1.8| 0:0] 251 17.210.001 mae aera 
Fl 321 77. 44| 12] 119) 38) 33 9} 56] 7-1] 179) 54] 3-8} 3-7] 2-5] 1-6] 0.5] 3-4] 22. 
Am 34| 48 74 I 6| 45 FIN 122 420 5:0 © 2 
Si 40| 77) Sof 13] 19| 57| ı2| 17] 40] 45] o| 39 
al 33| 65 99| 13| 24| 68 étre 5764.09 Oo} 39] 1.0) 8-7] 4.41 4.0] 0.0] 1.0| 4.2 
Sö 29| 84 44| II] 90} 35| 22 9| 47| 6.6] 107] 48 
Sm Et 26 22 5 3 27| 20 Lol) 47, 641 50| 4al Sigh) 27| 3:2 1 Go IS CS 
Sy 6| 18 18 9 I 7 6 4| 64| 6.6 *| 62) 8.51. 78.51 3300 tae ro cire: 
BH 8 35 26 7 4| 19 8 6| 46| 5-3 ol 62| 2:7] 2:5) 23) 0:8] 0.2] 00.3182.8 
Ha 19} 41] 41; 8) 14) 27) 7| 9 46] 5-9] II] 39 
Al 22| 70) 55) 14] 19) 43) 14] TO) Iı2| 6.2) 35) 62) 5.6] 1.2) 5.12 1.9] 0.5) 1.9] 94 
Hi 8 12 init 3 5 9 I 4 2 533 © 27 
Äs 23| 46 22 OT ITS 8 51 321 4:9 O| 34| 14. OL) 3.2) 8] FO. Ti 272] 070 
Va 23 8 o fe) I o o 3 2| 5-8 6 2 =“ % * * * * 
Li 56| 96} 900| 14] Iz] 470] 23] 66} 51] 4.5 oj 85] 3-6} 48. | 1.5| 74. I-5| 3.2] 9-3 
Ka 29| 29 13 5 Sl, 23 7 AlgES2 536 6) 14| 4.1] 3.4| 2.3] 0.8| 0.0] 0.9| 18. 
Ku 37| 34 18 6} 10| 16 2 4| 29] 5-4 oO Ta x x * * * * * 
Jy 280227 54 6 Yu LO 5 4| 60| 5.8 O|[ 22| o.2| 0.6| 1.1] 0.6| 0.4] o.1| 4.7 
Tv 14| 16 10 3 Alte 2 sl Sal eyes QT 
Od Ill 37 49| 12 733 I 6] 56] 5.6 61 60| 5:3] 14. | 72] 3.0] 0:8, 051255 
As 1018270 971 14) 29| 63) 26) 36| 67| 6.0 az] 7 2.7| 4.8] 4.5] 2.2| 0.7| 4.6] 22. 
Ty 15| 42 60| Io] 13| 44| Lo] ı3| 74l 6.6| 68 57 4:7) TOs | 72 7-.20022.51 6:0) 25. 
Ab Fe LS 60 O Oo} 50 2 9| 49 Me Sl ae 0.0] 4.I| 3-1] 0.0] 2.7| 20. 
Ed 25| 40| 124 4 2| 79) 48| 19] 77| 6.7] 106} 50] 6.7] 9.4113. | 4.1] 1.6] 0.8] 0.0 
Le 14| 67 32 el | 6| 28) 94] 5:5 60.73 
x “= 45 2 II de 25 9 81.0412 5:96.76 62 17. 5:5| 3:5| 2.8| 0.7] 4.7| 12. 
pa er a : N +? er By: ay 74 6.0 19} O61) 6.8) 4.5] 6-7] 3-9] 1.4] 10. 7.1 
Ci aes 2 4 3 A, = HOG 6.0] 14| 31/14. 8-1! 7.5] 1.1) 0-6] 0.7 5.0 

2 3 531) 22) 229321 43l7272]|7.021018.5 O| 13] 9-4] 4-1| 4.2] 3.5] o.7| 0.7| 9.8 
Fe 95| 134 59180) E47 1248 2.18 2073 10.321 75% O| 21| 4.8) 3.7| 2.5] 1.9| 0.9] 0.6| 3.6 


— No sampling 
* Sample discarded 
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Precipitation June 1956 (D 606) Air June 1956 (L 606) 
a mg/m? 3 |. g/m? (= kg/km’) 
© 
° [mm Gal pH lus 5 ö A 
S Cio 2 Na| K | Mg} Ca a RlQ| Ss | CI cl Na} K | Mg} Ca 
epee es re) tie ee) ie i lee) ba 
Ri 81 727 7 | 
Ah RSS) Ge Ol ast) ° 5-0), 14 9] 3-3] 4-3] 1.4] 2.2] 0.6) 0.3] 6.4 
Ki 731 223 27 5 2 | | EL Gee fo) ©} 5.0) 0)ETrS) 0] 202 "or 
Ar 57| 38 22 5 5| 31 8 O50} OS) 829] DL! 0 0 00 Er 80.7202 500.07 
Oj 83} 93] 46] 13] 28! 34] 2ıl 3] 1231 6.2) 29| 18 | 
Rö 48| 38 16 6 SE 27, 12 5l 56] 6.0 13 « * * * * * “ 
Of 68| 73 32 7 8 210 3012. 1312 74 6.0 17 10/2420] .0:0|.3.9| 2:3187.510.0:0|1,.:0:0 
Br 83] 41 10 co) 34 8 8 4| 29] 5.0 ol T2| 0.11.23] 2.o| 0.011 0.01.54 00 
a 102| 70 47 6| 10| 26| 22 2| 61| 5.5 fe} 8] 1.3| 0.0] 0.9] 0.7] 2.0] 0.5] 2.8 
H en. 
Sv 96] 32 151 ro) zT AS 9 Sie 371160 A321 LO 
Rä TA 020 33 20 Ar 29 4: 610707 76:6 273 271..4:0| 43.0|) 25010 1:08 0:0]2.0:.01 0533 
Am 85] 56 23 7A ee 7) Viet | ma 9| 45| 5-7 Z| - 201 32| 20! 1261”1.3100.6, 2.1) #7 
Sa 130] 40 no) 12 5 161 29 310 861 ser Ol 12) 7224| 70.0), 1.5| .0:6|0.0.010.0.8| 24.9 
Ul 89] 60 2 Pal) ISIN T2| 2 5| 49| 4.8 or 51 00 rl .0:.81°04104:5|,.0:6 
Er 62| 89 331 2ER 22710251451 151, Onl 5.6) To TO} TS 2:6 Tes) 1.8004 ..24|055 
St 107| 70 34| 14 Ale 22517. 231 12131267 5.4 fe) 141827 T0 2.8 0.91. 0.41) 0.5|| 44 
Fo 178] 69 73| 65! 46) 51| 32 8] 140] 4.7 o| 22 = cS * * Fe * 
Kvı 46| 637 781% “TAF #2 45| 20| 787| 6.0] 16] 161| 42 0.069. 9600.17 7.8 Tr. 
Kv7 49] 91 31 9 zu 23) 14| 1512 83| 4.6 ol 3188.41 72,8], 3:0] tae 0.518 7,9105: 
VN — — — 
La 35| 47 61 9147221, 37| 42 To) 03 CM 56-7733 5 5 * * = = * 
Bo 60] 71| 264| 14| 15| 160| 15) 23] 81| 5.6 4| 34| 0.8| 24 2.5| 18 a) ZEN ni 
Vi 44| 108} 894) 15] Io] 481] 28] 59] 80] 4.7 o| 105] 4.8] 19 1.8] 13 0.610. 2.41. 0.3 
Fa 103| 49 45| 15 AZO 44107517821 6.7 4012 TO 3.2|05.3,28, rs ord ON me! 
Fl FE 78 Sol el Ar aa del TSP ol Sr 201 9253) 3:71 771, 7.500903 .79 
Am 106| 77 67], 141 29] 53| 17) 1717 a5. 5.0 o| 16 
Si 7021705], 17921 729) 171707| 351 pz) 64| 4.6 o| 28 
Pi 761.220 1067|, 191 37 162) TG TI 49|| 4.8 Ol ZO) 5 1-3)\ 3.51 320 1-6) 035) oa] 23 
S6 55| 68 48] 12] I9] 30} 20 8} 49] 4.6 o} 311 —-} —| — — —| —| — 
Sm 2| 66 23) 72 2| 29| 35 8 76| 5.5 or ZO} Ol C201) 2.51, 0.816. OT Ouran 
Sy BT Asia 10) ss, as Ay) (ON 3 él ARLON CAN EE 
BH 49] 59 34 12 4) 30| 35) 13] 57| 5.3 or 2212.0) 2155| .2221 7.05. 0:010.0:3 2.188 
Ha —| — —| — 
Al 52| III AZ eral Sr DOW 22072710 5.8 71 30 2720| 1726| Ir) 0.401.929. 
Hi 44| 54 29] I0| 14| 66| 12 97 Azul 6.4 1719| zo — 
As 106| 82 ASIN 13) 1353| 28) TE 6| 45} 4.8 OS 5.202,8 2.107.510 0.8 77.2 075 
Va 47| 54 19 7 8} 13] Io 7 or 5.1 O| 19 4 = À + H 5 
Bt 53| 142] 1440| 12 5| 716} 50| 61| 83| 4.6 0723| 87114 TO 59.01) 2-01" O10} 20. 
Ka 44] 51 22 8 SITE | ET za a 7 Bi Lo 4 = = # = * 5 
Ku 47| 58 14 6 5 LO| To Te 531, 4.91 Ol) 16 = “à # : x + 
Jy 731072 43 7, 2a 64) Te BI SA 5-0 © cS .70410.2,4181951.0.716. 0.26 7:3|.556 
Tv Bauer gs) fOr TON 2) 20, 0177 57] Oza] 0) OP A Ble Ale Fl 
Od 711 95 Solel) 3281440 2710215 0670149 OMS 155227) 43|02.800:.681.4| 220 
As 63] 95| 118| 64] 62] 80] 50] 18] 68] 6.3] 29] 32] 0.8) 0.2) 4.1] 2.2] 0.2] 2.0] 6.1 
Ty 30] 64 337 16) 251 35) 13 9| 100} 6.4| 39] 36] 0.6| 21 SOLL 10.042 Metco 
Ab 34| 42 75 I o| 69| ızl 14| 95) 6:51" 58| 32] 8:3] 4.3] 3.8| 3.3| 0.0| 1.3) 2m 
Ed 58; 84 84 4 1| 56 ZN 421011 5:0 0251734 | 324 25S 0.5/8 0.0 E74] 79 
Le 73] 136] 127| 18] 24| 51] x18] 36) 173} ‚4.8 OW 361202 720) P45 7 MO YO) 03-2] 10: 
Ro 0612) EAS || 31] 23 20 29 E70 5 Oo AO TA 231. 23|00.782:0| 122 
NA Re 13h 221013817 270 2319134) 0.2) 0651 39) 3.3) 327) 4-7) 2.4 0.97 1.9) 13. 
Au A723 38 210227108 TE) ol 571072840541, 5-90. 2.72 0.50 7.2722 
Gü 78| 98 3302150. 28 0.251 217) 228120054 ©. 179.107] 24010320 7.210.0.3|03%7 in 
Fe 84| 129 36| 16| 29| 24) 14 9| 121| 4.9 00223122351 0:0] led! 7.91. 0.080.71247 
Sc ARI 772 ol LO 98 Lol) 16 1461 4.6 fo} IEC 


— No sampling 
* Sample discarded 
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a 


Precipitation July 1956 (D 607) Air July 1956 (L 607) 
© mg/m? Sı- “g/m? (= kg/km?) | 
rg © 
8 WC I cl: (Nes) Re iM | “Ca S|xlal S| Cl] | Na] K |Mg| Ca | 
Zu = A er | 
| 
Ri Bide ce) TEE Du aus EUR A aie ep BP Fe Cr a Se | 
Ki 69037 I4 I I 9 7 21052, Ossie 4 8 2361 4.3 29) 2502.20 200 
Ar 5222 II 4 3 9 7 al By 93 ate 8| ©o| 0:0], 5 1.22.0.7180320.0 | 
Oj 26| 19 31 3 8 14| 15 31° 531, 10.4) 40 2a 2.0 3 220]! (0.3/5 00/0) ro | 
Rö 60] 47 21 @| Bid] 75) 22 A 38 032 | En re, * “à 4 Hi = i = 
Of og Go) 22) All al) Re sea) u | © zZ we OR i = : | OS 
Br TO 32 9 4 Ally tea. 6 3108220 53 fo) 7| 2.01 2:91 2.0 0.6 0.512031 75 
AP les v lies), a 2 ale eh CDS EIRE. 
AH 4020) oa ae ea el al 25 Gurl] © ze z * + z + 
Sv 87) 33h 20) 63h ‚2 rl Sh enOoh on 162) 33 
Ra 5826 37 4 Oo] I8| 14 3181607|M6:8 92 27] 0:6) 2:6) 2,6) Nr 41 1.8002] 70% 
Aan 4) 58) 39) 22) Sl Bh 73) 7 ale 3256 Gol) ol He. kei Ta Fee : 
Sa 69] 49 I5 5 Ol Lo 8 6} 49] 5.3 Oo EE} 30| 2:0] 147) 0.81 0.6), 7.8) 89 
Ul 86] 66 | AG) ara eG ee 5| 36] 4.9 6 15| 24] 0.4| 1.5] 0-4) O:6/ oO FIG 
Er 51| 47 14 3 All Hi) 71 25 520 fe) 13] 2 293s 251 0:98 AsO Moses 
SIE 66] 50 15 4 |e Oe 6] 53] 5.4 fe) TE) 02 20" 670.808) 0/70 
Fo 72| 45 14 5 Bi 25 7 Si Asi) 15,5 fe) 9) 2A SRO) (222) EIS 
Kvı 50] 540 50 8} 52) 20} 36] 23] 498] 4.3 O| I14| 29. 3:9 706.4] 1-2]0751 DER 
Kv7 55.02 28 ee 3A 131012 7| 36| 4.6 Ol) 25| Ist! 00 3.8) 2.207.831272:608 
Vk 
La 50] 43 84 A A | ES 8 561 6.2 14, 2763.00 0.0 39 9:0, 2378022 
Bo Sul 83 2239 1410. 2711156 1213| 24170 02 6 «251 TS) CSI 125]j8 2.31 0.6 Wo sw rH 
Vi 33] 98| 2140| 12! 11/1140) 46| 148] 105] 5.4 Ol) 2081. 7.9778 1.9| 43. 2.410 5:81 RR 
Fa 50| 34 23 2 gl sae 6 Sl 355 Q| FTo| 4.0] 0.0) 3%), 2.8|,2.71.7.61 06 
Fl | 42 63 A SAE aes 8| 390] 5.6 Oh £2) Le TOM ng Era Oo) 
Am 102| 86| 195 9 Me TO AG! 5G, [6) I4 
Si 91| 72] 310} Io 6.270381 281 20 500 5% 00728 
BI AS IMC SCC TE TO 205) ro SSI NC CI 44l 551 56! 4.0], 3:21 72.210065] 06 | 
SÔ 24| 37 63 4) 1719| 37| 22) TT #58 .6.5| GO 34 | 
Sm 20154 22 5 DIM TO IEEE O2 AG 0331 To 24| O.5| 0.0 246 2.0107.21.04(73 | 
Sy 521 38 38 7) xt] 15] 28) Fly Fel Gal aol xe} IR sal * | 
BH 44| 46 99 8 4, 60 18) AIS 0! 28] 2.4| 2.8 1% 0.7) 0.3]. 1.1] 44 
Sk 42| 58 927279231 53] | DOW 72 2 735g Er 28 
Al 30] 60 88, 10, Toll on 231, 7) Sales 201 37 * * = =. * > 
Hi 22 027 51 5 24) 36 6 6} 34] 6.4 2 2 
As 80] 76 AST eS Er 15 8| 54| 4.8 fe) 20) Ind} O10}. 3.3) Lo TS pees) oe. | 
Va 58} 65 21 6 Bl) Le 16 SE ei] 2 Oo} 15 * * * * * * * | 
ET 103| 112] 1260| 23| 31| 790| 40| 87| x10] 4.6 O| 72/18. |267.| 2.41156.| 6.9] 22. | 20. 
Ka 38] 25 18 fe) Slush Le 7,172 6,5| So 2, * * * * * * * 
Ku 7350 8 I 3 7 7 32 510 3 5 * * * x * x * 
Jy 42532 2 I 2 7 6 Al) 4310 6:01 To TO} 1.3) 3.0| L.ı) 0.7| 0.6| Tr] 7.1 
Tv 56| 47 20 2 | | ar 5| 36] 5.6 O 10| 1.7| 9.5] 2.0] 1.6| ro) 4.5) 27. 
Od 5S 77 78 188221020961 2) 14| 98| 5.5 0: 35] 0-4] 2:4] ©:5| 3.6] 0.97] 1.2] 5.6 
As 38] 96) 603] 13] 25} 390] 36} 222] oo! 6.0 T4 821 2.8] 2.31 Sly G6. Sie 3.4| 8.1 
Ay 3210 91217510 1917733103491, SAIT Os 6.8) 26017 752] 2a Br on TANT OMS 
Ab 127101722 2130 I Se | 5 a). 291 7.0) 5 Al ane 
Ed 142| 140 Sl el | || | 26 Oo 22] 5.6) 0.01 221 2.810058 Dal | 
Le 72| 146 83| 2 245, Tel 3310238) 9:3 fo) 29} 17. 9.4 5:6] 2.5], 2.31. 9,8170. | 
Ro GI| 101 87], 2003002677. TO TOI 5.0 fo) 28] IL. 52 .35| 3:5] Plot 1730| 228 | 
NA 97| 159| 398] 20] 16] 212| 34| 36 175|. 6.0| 17 33) 9:31 754 565) 2 TITI ONE 
Au 125| 173 42] 29| 54| 16| 19] 29] 150] 5.8 8 77] 0,9, 2.61 585 Tr Elo ll 5.5 
Ba 145| 100 30 LOW tris 30 |e Salle ess o 9| 08:0) 48 29] LEERE aoe 
Fe 162| 136 45| 23] 20) 23| 16) 13] 123] 5.9 8 TO} 2.5| 2.6) 1.4) 1.2] 0.711 0.6) 80 
SC 68) 113| 346| 22| 24] 202| 32| 26] 159 5.71 IO} 44) 1.6| 8.5] 8.0] 7.5] 1.8] 3.6] ro. 
BV 180| 33 66| 35) 59] 94] 38] 66| 308] 5.5 ol 23 * * * * * * * 
Ho ee, el) Sa 
Bn 541 378| 11] 40] 65| 40| 28| 63 7041 6.4| Sol 62|27 5.8| IL. 8.2| 3.9| 2.8] 31. 
— No sampling Uppsala, October 20, 1956 
* Sample discarded Gunnar Brodin 
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NOTES 517 


The chemical network has been extended with new stations i Belgium and Holland in collaboration with 


DT: Grandjean (Institut Royal Météorologique d’Uccle, Belgium) and Dr. Hauer (Koninklijk Nederlands 
Meteorologisch Instituut de De Bilt, Holland), 


The analytical determinations are carried out by Dr. J. Bouquiaux and Mr. M. Mertens at the Institute of 
Publique Health, Brussels. 


Coordinates and altitudes of the stations are given in the table below. 


Precipitation March 1956, D 603 Air March 1956, L 603 
2 mg/m? =: A 4: ae ug/m? (= kg/km?) 
S |mm A | 4 pial PERS a 
S Cl oO ch Na} K | Mg| Ca Q x|a| S Cl = Na | K|Mg| Ca 
A zZ an! zZ 

31| 70 56 4 61721 sl I er Air fo) 51| 63. 9-1 0.1:0183:0.72510,.0.807.4 
SA 26] 56 50 3 3 ZT 3} 59| 4.6 ol 44/17. 80142-01221 08207) 20.2 
B 133| 150] 144 7 9 Bien Le 7| 83] 4.6 010 222155.01,.23> 82.216 3:00 WO pee 
D 55] 50 52 3 Al 27 7 an Ze) fo) 18} 16. 2 01S). O} a O.8 E70 en) 
DB | — = le ails) el | Tee 

April 1956, D 604 April 1956, L 604 
U 28] 116 59| 3838| 67)" 271 12 3] 114| 4.6 02.641572 ite 1702|, 2.07170|70.724.2 
SA 26| 83 el = 351758 8 4| 68} 4.3 o] 57/15 6.5], 5.01, 7.6) 20.420:.4| 241 
B 68| 139 78| 66| 106} 48| 16 7\ 2152| 5.2 ol 39) 14 9:2| 3:9| 3.7] -0.6| 0.7| 4.7 
D 39| 36 43\ 431 7219|. To 8 96| 4.9 (o] 47| 11 54105018 2271 OO) Onl 74 
DB 32| 102 7D eae 47| 38| 11 4| 78| 4.4 Ol 541262 || 722 7.8] 0.410 :01M 02170 

May 1956, D 605 May 1956, L 605 
U 28; X 49] 40} 40 21 261 2% 1651 74-9 ol 73| 38. ee re 8 
SA So 75 2 2a a ME 7| 86) 5.4 OSG) CIN 222117220 Er 0|N Oss Onl RS 
PSS) 23H 1111 O212-03! «641: 26 9] 187| 5.0 fo} 2015-3 19. [5.9 3:91,7.61.00.% 
D Ay eZEOle 102)" 83| Loo} °45|.26| "ı5| T9T| 15.3 O| 40|14. 4.8] 4.0] 1.5| 0.6] 0.8] 2.0 
DB 28| 118 71| 28] 209] 56} ı15| 15] 90} 5.2 o| 54] 20 2.9 29:21 ,.12710.0:01,7.3 102% 

June 1956, D 606 June 1956, L 606 
U QI} 205 83} 85] 108| 46] 20] 13| 112| 4.4 00239, 782.742 || OV ees enO.8 No ne 2.0) 
SA | ror] 184] 165} 86] 89! og! 32 7| 90} 4.5 el. Sil ARB el Giga) a os | ate 
B 145| 239 66| 109| 133] 51| 18 6| 89] 4.4 or ur 199.00%2|0 0.4) 0.3) 7.2 
D 2| 89 39| 41 2,25 9 4l 69] 4.9 © 23) ro; HOME O|N2-41N 052 les 127. 
DB 57| 129! 75} 50] 58] 55| 13 6| 90] 4.6 ol 3815 LS.) 9:90 2:00. 0.7 et 2) aed 

July 1956, D 607 July 1956, L 607 
U 79| 12 28] 49] 57| 30| Io] 1o| 108} 4.7 0l0.237]18.8:2 | 2.911.0:5 1774|. 1:2120:9]0252 
SA 46] 71 HO ZA dls ind) > aA ell Zea 00 Aer SSI 0) 51.31.07 Lose 20:7, 
B 102| 169 Fi 381391387 2810277 34:7. 08.26) 21:6) 2-3) 2:7120:9| :0.5|.:0.2|..%0 
D 47 36 ZG" wales} | Ko] | ZT 6 5] 521 5.4 fe) TON Hoey| 8232 OM GE 
IDS \er24\242). 185] 60) 59, 131|- 33| 225) 93] 4-4 Ol 938] 73:917 4:9] 7-0) e250) no) 77] 21 

Brussels in October, 1956 A. Mertens. J. Bouquiaux 


Laboratoire de Chimie et de Physique, Institut d’Hygiéne et d’Epidémiologie, 14, rue Juliette Wytsman, Brussels 


Coordinates of New Sampling Stations for Air and Precipitation (cf Tellus 8, p. 825) 


ee ee Cade | ee | LE Le pou 

an nn ee | ee re ee een 
Baldenwegerhof (replaces Güntersthal Gü) .. Ba 370 47.8 7.9 
Sehiles wa Guar. ga Sara clos) sehe pete «ee Se, 48 54-5 9.6 

| Braunschweig— Völkenrode ................ BV 83 52.3 10.5 
Hoberpeissenbeusg se ee. er. Ho 900 47.8 I1.o 
Dons coton ee oe te ais Bn 62 50.7 Te 
Skurup (replaces Hammenhög Ha) ........ Sk 40 55.5 13.5 
LÉGER CR EC OP COR EE U 100 50.8 4.4 
BHO MOMENT Se Re SA 14 Bee 32 
BOURSE LES 0e B 694 50.5 6.1 
ID'OUHDES a chorea eGR ee ee eioe OncronCr cece D 224 50.1 4.6 
DIS ABE Yo oo o en Ci RO CTC rn eee DB Io aa 5.1 


JU RENE So 0 OOS) AE SCHEN EB Ele RER SRE) ig, en RS RE ES 


CO,-values May—July 1956 in Scandinavia. 
(Cf. Fonserius, KOROLEFF: Tellus 7, p. 258—265) 


May June July 
em | | 
Sta- 3 : a 8 = = 3 a 
-— a © an 
tion | 2] © 3 = [Sel2|o = aioe] © E £108 
a = 8 ö=|ä|° = 5 Jozıa| ° = QU DE 
ie = El > E = E 
Ti 4.01 WSW 6 TOES ee 1.0] W 6 Io] 3-18 || I 
Ab | 10 2.0| SE 3 013734 || TON ier Ans) E 2 o| 3.31 || 10 13-8] N I 2| 2.98 
20 2.5| N 3 IO} 3.18 |20| 10.5| ESE 2 |=° 7| 3.09 || 20 7.8) NW 3 IO] 3.02 
I 7OINSE $\\Set2i|e |) hor) Calm. TONNES Fi 15.5| E 2 8| 3.14 
Ojo 9.2| NE 2 ı|3.35 | ro | 20.3] E i o| 3-36 || 10 
20 7.6| NW » 8|3.41||20| 16.8] E 3>= 2313,72) 200 270.5, K@alım IO] 3.02 
ie 8.0] SW 4 5} 3.09 || ı | 10.0] NNW ı |@? 10] 3.12/]| ı 16.5| NW I 5| 3-08 
Bry) ro TO} 22-0]. Siw I 2| 3-39 || 10 21.0| NNW 4 2| 2.98 
20 5.0] NNW 1 9| 3.30 120 | r0.0] N 5 10| 3.02 || 20 13.0| N 4 10| 2.90 
I 10.3| ESE 2 6| 3.27 || 1 | 21.4| W 4 10| 3-07 
Sal Fo 15.0] ESE 2 81 3.41 ||} 10] 23-0] N a0 233272 
20 9.6| WNW 4] @ 10) 3.29 |20 | 16.8] SSW 2 10| 3.24 ; 
I 4.9| N 2 ala or DE 3 = 9|3,51|| I 
Er | tro 14.21, ESE I o| 3.21 || 10 | 21.8] N 2 O| 3.15 || 10 20.5| NE 3 I] 3.38 
20 DEREN! S| 3522111202], POLS 3 8] 3.18 || 20 24.2| W I 3| 3,06 
I 6.4| NE 3 ©! Ses DN ck | 24-2) 5S Wi I O153:A.0 ME 20.81 SW I 0|-3:27 
Fl | 10 10.0] S I] @ 10|3.30||10 | 16.8] N 2 5| 3.32 || 10 23.2| NE fo) o| 3-08 
20 12.4| N I 21:23 |) 20: | 73231. Oy) I |@ 10} 3.33 || 20 23.2] NE o o| 3.02 
it 6.0] NE 3 ro| 3.52.72 | 27:01 Ww 6 ESC TAUX 21.0| NE 3 2| 2.90 
Pl | ro 9.0] W 2 9| 3.33 || 10 |10 20.0] NW o| 3.02 
20 14.0]. W I 213.32 120 |20 20.0| W 2 I| 3.00 
I 7.6| NE I O| 3.301 || SL |, 2S-0] SW: 2. |@ 193.431 7 13.3) NE I| @ 10] 3.38 
Ka | 10 11.4| NE 4 0| 3.39 |Io| 20.1| E 5 0} 3.51 || 10 16.2) N 2 10| 3.25 
20 9.5| WNW 2 10] 3.20 120 | 12.8] W I |@ Io} 3.48 | 20 18.1] N 3 6| 3.20 
I 8.0| E 2 ©} 3247 BE | 27:08 I Io 3.36 | I 14.3| NE 2] @ 10) 3.40 
Ri | 10 10.8] NW 2 IO] 3.38 || 10 | 19.9] E 3 O} 3.60 | 10 12.8] N 3] @ 1013.50 
20 8.4) NW 2 10| 3.37 |20 | 20.4] SE 3 0| 3.70 || 20 12.6| N 2| @ I0} 3.44 
| 
I TIME 2 0| 3.26 || I | 19.4) S 3 IO] 3.20 | I 12.0| N 3| @ 10| 3.40 
Lu |ıo 7320 NINE 23 0! 3.41 || 10} 17.1] ENE Ss 5| 3.39 || 10 14.2] NNE 2 I0| 3.27 
20 9.6] WNW 3 10| 3.27||20] 18.8] ESE 2 10| 3.30 || 20 16.1] ENE 3 Io] 3.14 
I 6.0] E 3 SH 2635 IE OURS [6) o| 3 2| x 14.0] N 2] @ 10} 3.30 
Tv. |: to 12.0] S I 513.30 || 10 | 20.0] E 6 o| 3.42 || 10 21.0] WNW 2 o| 3.26 
20 6.0] SW 2 71321120 || 17.0 N fo) 10| 3.39 |20 17.01 SW 2| @ 10|3.32 
I 6.9] NNE s| =°10|3.26|| 1 9.31 WSW 4 |[—°10| 3.38] ı 9.3] W I O° x6 3.08 
Bo | 10 10,0] N 5 113.47 110 | 19-7] NNW 2 0| 3.08 ||10| 23.8] SE 21 @_ 3| 3-23 
20 4.9] NNE 3] =° 3] 3.06 ||20 | ı1.5| B 2 |=° 5] 3.11 || 20 10.0] NW 2) oO) 3-80 
I 8.8) NNE 2 313-30| I] 18.4) SW I 0] 3-16)| I 13-3] W 3 OY lessee ace! 
Va | 10 6.9| NE 2|=°10| 3.43 ||1o | 23.0] ENE 3 2| 3.28 ||10 
20| 10.0] W 4 5| 3-14 || 20 |20| 20.3] NE 2 o| 3.00 
| 
R ® 
I 11.6] SW 2), 513-76| z| 9-4] W s |v 8|3.06|| ı | 15.7| Calm 10| 3.30 
Be | 10 9.4| SSE 2|v 7[3.44|10| 22.2) NNW 2 3|3-17 |1o| 15.4| WNW 2 SAS 
20 8.1] WNW 1| @’10| 3.35 ||20 | 10.1] NW 5 10| 3.43 || 20 21.4| SSE I O| 3.17 
4 I 6.0| NNE 3 7\ 3-22 || 1| 16.4) W 4 TI a02 0 7 19.5| SE 6 O| 3.01 
Od | 10 13.5| W 4 SI 31110 22:6) IN 3 > 8| 3.15 || 10 22.218 4 5| 3.07 
20 12.6| WNW 3 713 30:11 poses 20 | 
I 5.9| NE 4 813.32 || I | 13.1] NW 3 5[2.99 | I 22, 20SE 3 10| 3.06 
As |ıo 13.0] W 4 5| 3-21 || 10 e 10 22.0| SE 3 10| 3.04 
20 0.2] W 4 6| 3:33 [20 | 14.0] W 5 [V1 201 20 PS AIISE 3 313.10 


Symbols used for description of weather conditions, are those, internationally used. The numbers in the weather 
column give the cloudiness and the wind force i Beauforts. The Station formerly listed at Ultuna, was move 
1 May 1956 to Erken lat. 59.9 N, long. 18.7 E. Both Stations were op erated from ı May 1956 to 30th June 1956. 
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